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1.  INTRODUCTION 


It  is  generally  taken  for  granted  that  the  spectrum  of  light  docs  not 
change  on  propagation.  This  belief  is  implicit  in  all  spectroscopy.  In  the 

last  few  years  it  was  discovered  that  while  such  spectral  invariance  holds  for 
radiation  from  commonly  used  sources  it  does  not  hold,  in  general,  because 
the  spectrum  of  emitted  radiation  depends  not  only  on  the  source  spectrum 
but  also  on  the  coherence  properties  of  the  source.  It  was  found  that  in  some 
cases  the  induced  changes  may  be  appreciable,  possibly  resulting  in 

frequency  shifts  of  spectral  lines.  It  has  also  been  demonstrated  that  similar 
spectral  changes  may  be  produced  by  scattering  on  random  media  whose 

dielectric  susceptibility  fluctuations  are  appropriately  correlated. 

Very  recently  it  was  shown  that  under  certain  circumstances  the 
frequency  shifts  generated  by  dynamic  scattering  may  imitate  the  Doppler 

effect,  even  though  the  source,  the  scallerer  and  the  observer  arc  all  at  rest 
relative  to  each  other. 

The  results  appear  to  have  considerable  potential  for  both  scientific  and 
technological  applications,  for  example  in  the  following  areas: 

(1)  Development  of  new  communication  techniques,  some  of  which 
could  be  employed  for  covert  operations. 

(2)  Improvements  in  the  accuracy  in.  standards  determination  for 
spectro-radiomelric  scales. 

(3)  Improvements  in  the  accuracy  in  tracing  of  satellites  by  laser 


beams. 


(4)  Testing  of  theories  of  turbulence. 

(5)  Astronomical  applications,  especially  in  connection  with  inter¬ 
pretation  of  some  puzzling  features  observed  in  the  spectra  of 
various  astronomical  sources,  particularly  quasars. 

The  latest  paper  entitled  “Correlation-induced  Doppler-like  frequency 
shifts  of  spectral  lines”  by  E.  Wolf,  published  in  the  November  13,  1989  issue 
of  Physical  Review  Letters  (Vol.  63,  220-223)  has  attracted  considerable 
attention,  as  is  evident  from  reports  that  have  since  then  appeared  in 
various  magazines  in  both  the  USA  and  abroad.  A  selection  of  them  follows. 

In  this  booklet  we  also  include  a  list  of  references  to  all  the  scientific 
papers  which  have  been  published  on  this  subject  so  far,  as  well  as  copies  of 
papers  which  describe  experimental  verification  of  the  theory.  Most  of  the 
theoretical  research  was  carried  out  under  the  sponsorship  of  the  ARO-UR1 
program  at  the  University  of  Rochester,  by  Professor  Emil  Wolf  and  his 
students  and  other  collaborators. 
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Physicist  Emil  Wolf 
extends  theory  of  spectral  shifts 


A  University  professor 
has  extended  his  theo¬ 
ry  challenging  a  tradi¬ 
tional  tenet  of  modern  physics 
—the  cause  of  spectral  shifts. 
The  work  by  Emil  Wolf,  Wilson 
Professor  of  Optical  Physics  in 
the  College  of  Arts  and  Sci¬ 
ence,  could  have  repercussions 
in  many  areas  of  science  and 
suggests  our  view  of  the  uni¬ 
verse  may  need  to  change.  Al¬ 
ready,  other  scientists  across 
the  globe  are  applying  his  theo¬ 
ry  to  long-standing  problems  of 
physics. 

Wolf’s  latest  paper,  “Correla¬ 
tion-Induced  Doppler-Like  Fre¬ 
quency  Shifts  of  Spectral 
Lines,”  appears  in  today’s  issue 
of  Physical  Review  Letters.  It 
represents  an  expansion  of  a 
principle  he  described  in  earlier 
work.  In  the  new  paper.  Wolf 
shows  how  something  other 
than  gravity  or  the  Doppler 
Effect  can  produce  dramatic 
changes  in  the  behavior  of  light 
or  sound;  in  fact,  the  mecha¬ 
nism  Wolf  has  discovered  can 
completely  imitate  the  Doppler 
Effect  to  any  magnitude. 

Scientists  have  known  for 
decades  that  light  from  differ¬ 
ent  sources  has  different  spec¬ 
tral  patterns  and  that  these  pat¬ 
terns  can  help  identify  the 
nature  of  a  source.  By  studying 
the  spectral  patterns  emanating 
from  a  star,  for  example,  as¬ 
tronomers  can  determine  its 
chemical  composition. 

Scientists  have  also  known 
that  light  spectra  can  be  shifted: 
The  spectral  pattern  that  is 
characteristic  of  a  certain  chem¬ 
ical  at  a  certain  wavelength,  for 
example,  can  be  shifted  to 
shorter  (bluer)  or  longer  (red¬ 
der)  wavelengths.  Scientists 
have  long  assumed  that  only 


Three  years  ago  Wolf  theo¬ 
rized  a  new  cause  for  spectral 
shifts.  What  has  come  to  be 
known  as  the  “Wolf  Effect”  or 
"Wolf  Shift"  has  since  been 
confirmed  several  times  ex¬ 
perimentally.  In  his  latest  paper, 
Wolf  shows  that  the  Wolf  Effect 
can  shift  spectral  lines  as  com¬ 
pletely  as  the  Doppler  Shift.  In 
addition,  the  shift  is  consistent 
—each  line  within  a  spectrum  is 
shifted  in  the  same  proportion. 

“This  mechanism  can  exactly 
imitate  the  Doppler  Shift,  using 
known  principles  of  physics 
and  optics,"  Wolf  says.  “The 
redshift  can  be  as  large  or  as 
small  as  you  want  -  in  principle, 
it  can  have  any  magnitude." 

Wolf’s  theory  could  have 
consequences  in  a  wide  range  of 
areas,  including  astronomy, 
standards  specification,  coding 
of  signals,  and  satellite  tracking 

Previously,  Wolf  showed  that 
the  manner  in  which  atoms  in  a 
light  source  are  ordered  (source 
correlation)  affects  the  way 
these  atoms  emit  light  and  the 
way  this  light  travels  through 
space.  When  the  fluctuations, 
or  light  emissions,  from  these 
atoms  are  neither  fully  ordered, 
nor  coherent  (as  in  a  laser),  nor 
fully  incoherent  (as  in  a  light 
bulb  or  candle  flame),  shifts  can 
occur.  Spectral  shifts  occur  in 
the  realm  of  partial  coherence. 

To  explain  partial  coherence. 
Wolf  uses  this  analogy:  “Picture 
a  contingent  of  1,000  soldiers 
walking  across  a  bridge,”  Wolf 
says.  “If  all  the  soldiers  are  walk¬ 
ing  in  step,  thar’s  an  example  of 
coherence.  If  they  are  all  drunk 
and  wander  randomly  across 
the  bridge,  that’s  an  example  of 
incoherence.  But  if  some 
walked  in  step,  and  some  wan¬ 
dered  across,  that's  an  example 
of  partial  coherence." 


Despite  confirmations  in  the 
laboratory,  it  is  not  known 
whether  and  under  what  con¬ 
ditions  the  Wolf  Effect  occurs 
in  nature.  Wolf  and  several 
associates  are  continuing  their 
investigations.  A  graduate  stu¬ 
dent,  Daniel  James,  has  devel¬ 
oped  a  model  of  a  medium  giv¬ 
ing  rise  to  such  shifts.  Whether 
such  media  exist  in  nature  is 
now  being  studied  by  Wolf, 
James,  and  University  as¬ 
tronomer  Malcolm  Savedoff. 
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Newly  Discovered  Property  of  Light  Spectra 
May  Resolve  Major  Dispute  in  Astronomy 


■y  um  a.  McDonald 

Aa  unusual  effect  of  light  discovered  by 
•  University  of  Rochester  physicist  may 
help  astronomers  resolve  (he  punk  of 
why  some  objects  in  the  universe  that  are 
connected  appear  to  be  (raveling  at  differ¬ 
ent  speeds. 

In  a  paper  published  in  this  week's  issue 
of  (he  journal  Physical  Review  Letters. 
Emil  Wolf,  a  professor  of  optical  physics  at 
Rochester,  contends  that  (he  newly  dis¬ 
covered  effect,  which  he  predicted  three 
years  ago,  may  exist  on  a  much  larger  teak 
than  has  been  observed  in  laboratory  ex¬ 
periments. 

If  so.  it  may  explain  a  major  puzzk  that 
has  plagued  modern  astronomers  far  (wo 
decades. 

Measuring  Light's  Characteristics 

The  problem  concerns  quasars .  unusual¬ 
ly  bright,  compact  sources  of  rutiaiion  lo¬ 
cated  in  the  centers  of  ancient  galaxies.  ' 
(hat  are  thought  to  be  among  the  most  dis¬ 
tant  objects  in  the  universe.  By  measuring 
the  spectral  characteristics  of  the  light 
emitted  by  quasars,  astronomers  have  de¬ 
duced  that  many  quasars  are  near  the  edge 
of  the  expanding  universe,  traveling  out¬ 
ward  at  incredible  velocities. 

Those  deductions  have  been  baaed  on 
(he  fact  that  the  spectra  from  quasars  are 
typically  shifted  far  to  the  red  end  of  the 
light  spectrum.  Scientists  know  that  the 

■"■■w**  “frt  •M'-;  ft* 


EmR  Wolf:  Hit  prediction  would  support 
astronomers  who  believe  quasars  art 
receding  neither  aa  raptdty  nor  ea  far 
from  Earth  aa  their  red  ktffts  Mlcate. 

lance  and  speed  of  the  object  receding  from 
Earth.  This  is  because  of  the  Doppler  ef¬ 
fect,  a  property  of  sound  and  light  demon¬ 
strated  by  the  sound  of  a  train  whistk  as  it 
travels  past  an  observer— becoming  high¬ 
er-pitched  as  it  moves  toward  the  observer 
and  lower  pitched  (or  shifting  toward  the 
Conrimir^yvivffeariAk  > 


Continued  from  Page  Ad 
other  end  of  the  frequency  spec¬ 
trum  I  as  it  moves  away. 

Physicists  have  long  known  that 
the  Doppler  effect  and.  to  a  lesser 
extent,  gravity,  can  produce  red 
shifts  m  light  spectra  In  recent 
years,  researchers  in  the  United 
Slates.  India,  and  Italy  have  con¬ 
firmed  experimentally  that  the  ef¬ 
fect  predicted  by  Mr  Wolf,  now 
known  as  the  "Wolf  effect. "  also 
can  affect  red  shifts  But  because 
the  Wolf  effect,  like  gravity,  is  so 
smalt.  astronomers  have  generally 
discounted  it  in  their  calculations 
In  his  paper,  however.  Mr.  Wolf 
shows  how  the  Wolf  effect  could 
produce  a  red  shift  as  large  as  the 
Doppler  effect  for  objects  such  as 
quasars  He  said  in  an  interview 
that  some  physicists  had  already 
begun  testing  his  theory  see  if  il 
holds  up  experimentally. 

If  true.  Mr.  Wolfs  newest  pre¬ 
diction  would  bolster  the  claims  of 
astronomers  who  believe  quasars 
are  receding  neither  as  rapidly  nor 
as  far  from  Earth  as  their  red  shifts 
appear  to  indicate. 

it  would  also  suggest  that  astron- 
omen  may  have  miscalculated  the 
age  of  quasars  and  theii  distance 
from  Earth — long  considered  an 
important  yardstick  for  measuring 
(he  size  of  the  universe 
The  question  of  whether  red 
'  'Muffs  reflect  real  coamofagicil  dis¬ 


tances  has  been  one  of  the  central 
issues  of  debate  among  astrono¬ 
mers  for  the  pasi  20  years 
The  debate  has  become  particu¬ 
larly  heated  since  the  discovery  of 
quasars  connected  to  galaxies,  in 
which  the  quasars  appear  to  be 
traveling  away  from  Earth  at  a 
much  faster  rate  than  the  galaxies. 


The  quetfon  of  whether 
red  shifts  reflect  real 
coemojogcal  dtoteoc— 
hee  been  debated  by 
eefronoroen  for  the 
peet20  yeera. _ 

an  incongruity  that  astronomers 
have  been  unable  to  explain. 

"That's  the  conflict  that  has 
brought  us  to  this  point."  said 
Geoffrey  R.  Burtxdge.  a  professor 
of  physics  at  the  University  of  Cali¬ 
fornia  at  San  Diego  who  is  one  of 
the  main  proponents  of  the  non 
cosmological  explanation  of  red 
shifts  "The  data  keep  piling  up 
and  people  keep  arguing  shout  the 
dilemma  " 

Mr  Wolf  said  the  Wolf  effect 
might  esptain  the  dilemma,  since 
receot  calculations  led  bun  to  con- 
‘bm  Wr  uriic  Ur  U  IDs  than 


have  been  observed  experimental¬ 
ly  could  be  produced  when  a  cer¬ 
tain  type  of  light  passes  through  a 
secondary  medium,  such  as  the  at¬ 
mosphere  that  surrounds  a  quasar 
In  the  laboratory,  the  Wolf  effect 
produces  a  slight  red  shift  in  light 
that  is  partially  coherent  Such 
light  is  produced  when  the  fluctua¬ 
tions.  or  light  emissions,  from  the 
atoms  in  a  light  source  are  neither 
fully  ordered— or  coherent,  as  in  a 
laser — nor  fully  incoherent,  as  in  a 
light  bulb,  candle  flame,  or  star. 

Source  of  Radiation  la  Unclear 
Although  Mr.  Wolf's  model  as¬ 
sumes  the  light  from  quasars  it  not 
stellar  and  probably  partially  co¬ 
herent  .  astronomers  say  the  source 
of  radiation  from  quasars  is  not 
well  understood.  But  if  the  light  is 
partially  coherent.  Mr  Wolf  laid, 
as  the  fluctuations  from  the  aec- 
ondary  source  become  greater,  so 
loo  could  the  red  shift,  until  it  ap¬ 
proximated  the  Doppler  effect 
“If  one  provides  a  good -enough 
model  (or  the  atmosphere  around 
quasars,  one  could  quantify  it."  he 
Hid 

Mr  Burbidge.  who  has  exam¬ 
ined  Mr  Wolf's  paper,  said  he  be¬ 
lieved  the  concept  was  worthy  of 
further  examination  and  testing 
"My  feeling  is  that  he  may  be  on 
,  (he  pray  to  hetpipg  us  understand 
these  rod  shifts."  fee  (aid. 
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Stellar 

‘yardstick’ 

disputed 

By  David  L.  Chandler 

GLOBE  STAFF 

New  work  by  a  respected  physicist 
could  make  the  almost-sacred 
yardstick  by  which  the  universe 
has  been  measured  seem  rubbery 
and  unreliable,  forcing  a  reexa¬ 
mination  of  some  central  beliefs  of  modern 
astronomy. 

In  a  paper  being  published  today  in  the 
journal  Physical  Review  Letters,  Emil  Wo£T 
of  the  University  of  Rochester  in  New  York 
casts  doubt  on  the  accepted  view  that 
changes  in  the  color  of  light  are  an  absolute 
and  accurate  indicator  of  astronomical  dis¬ 
tances. 

The  new  "Wolf  effect"  could  help  to  ex¬ 
plain  some  strange  observations  involving 
distant,  highly  energetic  objects  called  qua¬ 
sars,  Wolf  said.  But  some  other  physicists 
and  astronomers  say  the  effect,  while  prob¬ 
ably  real  in  principle,  is  likely  to  be  irrelevant 
in  the  real  world  of  astronomy. 

The  paper  strikes  to  the  heart  of  one  of 
astronomy’s  most  cherished  laws,  a  principle 
laid  down  by  astronomer  Edwin  Hubble  in 
1929.  Hubble  discovered  the  expansion  of  the 
universe  -  the  fact  that  all  the  distant  objects 
we  see  in  the  sky  seem  to  be  racing  away 
from  us. 

Hubble's  law  says  that  the  farther  away  a 
galaxy  is,  the  faster  it  is  going.  Therefore,  its 
speed  can  be  used  as  a  reliable  yardstick  to 
measure  its  distance  from  Earth. 

And  speed  is  a  relatively  easy  thing  t» 
measure  for  distant,  luminous  objects  like 
galaxies  -  which  Hubble  proved  to  be  island 
universes  containing  billions  of  stars  each.  To 
measure  their  speed,  astronomers  need  only 
spread  out  their  light  through  a  prism  or  oth¬ 
er  device  to  separate  out  the  rainbow  colors 
that  make  up  white  light 

This  produces  a  spectrum  containing 
characteristic  lines  of  much  brighter  light. 
These  lines,  produced  by  gases  in  the  stars 
that  emit  light  of  particular  colors,  have  dis¬ 
tinctive  apadngs  bke  the  barcodes  on  super¬ 
market  products  -  unique  “fingerprints"  that 
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The  effect,  Wolf  Mid  In  an  inter-  It  is  farther  ewmy  than  its  it-  ka,"  but  “whether  or  not  it  will  really  tnfl  as  a  demonstration  of  nature’s 

slew  last  week,  means  that  if  tight  tached  galaxy,  and  that  the  objects  explain  all  of  these  things  we  are  capacity  to  surprise,  he  said, 

peases  through  certain  kinds  of  mat-  realJy  are  connected.  But  most  as-  finding,  I  don’t  know.  It’s  certainly  a  “Amongst  all  those  billions  of  stars, 

ter  it  could  end  up  with  a  redshift  tronomers  believe  such  pairs  are  step  in  the  right  direction.  To  me,  It  to  me  it  would  be  a  miracle  if  there 

Indktingukhable  from  tisft  produced  '  aimply  chance  alignments  efid  the  looks  promising."  wen  none  that  shewed  this  effect" 
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Expanding  a  theory  for  shifting  starlight 


If  Oscar  Peterson  were  running  along¬ 
side  a  speeding  Ella  Fitzgerald  as  s!ie 
belted  out  an  A,  he  would  hear  an  A.  But  if 
the  pianist  took  a  breather  while  singing 
Ella  raced  on.  the  Doppler  effect  would 
make  the  note  sound  lower  to  him. 

For  close  to  a  century,  scientists  have 
known  that  the  wavelengths  of  light  from 
a  rapidly  receding  source  in  outer  space 
also  get  Doppler-shifted  to  lower  frequen¬ 
cies  and  so  appear  redder  than  they 
would  if  their  source  remained  station¬ 
ary.  Moreover,  researchers  have  assumed 
that  this  so-called  redshift  can  result  only 
from  Doppler  motion  or  from  a  gravity- 
based  mechanism  also  uncovered  early 
this  century. 

Not  so,  contends  physicist  Emil  Wolf  of 
the  University  of  Rochester  (N.Y).  In  1986, 
he  published  a  theoretical  sketch  of  a 
third  mechanism  that  could  account  for 
part  of  the  redshift  of  light  from  exotic 
cosmic  objects  such  as  quasars.  Several 
researchers  have  since  reported  labora¬ 
tory  confirmations  of  the  process,  which 
some  of  them  refer  to  as  the  “Wolf  shift" 
(SN:  9/13/86.  p.166;  7/11/87,  p.22). 

Now  Wolf  reports  an  important  exten¬ 
sion  of  his  theory,  suggesting  the  process 
could  account  for  arbitrarily  large  shifts 
toward  either  the  red  or  blue  end  of  the 
electromagnetic  spectrum. 

In  its  earlier  form,  the  theory  proposed 
a  mechanism  that  could  produce  small 
redshifts  in  the  spectra  of  light  emitted 
from  certain  exotic  sources.  Wolf  sug¬ 
gested  at  the  time  that  the  shifting  mech¬ 
anism  could  emerge  physically  from 
partially  synchronized,  or  coherent,  fluc¬ 
tuations  in  the  wavelengths  of  light 
emitted  from  the  countless  individual 
atomic  and  molecular  "microlamps’  that 
make  up  such  a  source.  As  these  emis¬ 
sions  travel  through  space,  their  original 
spectrum  would  appear  to  shift  in  the 
same  way  as  Doppler-shifted  light. 

Though  Its  astronomical  conse¬ 
quences  remain  unknown,  the  theory 
could  change  estimates  of  the  size  of  the 
universe  and  help  explain  some  anoma¬ 
lous  astronomical  observations,  Wolf 
says.  However,  astronomers  have  not 
rushed  to  adopt  it.  No  known  light  source 
has  components  that  display  the  required 
correlated  fluctuations,  Wolf  notes.  He 
also  blames  the  theory’s  unorthodoxy 
and  its  arcane  mathematical  formulation 
for  its  limited  consideration  by  astrono¬ 
mers. 

In  the  updated  theory,  described  in  the 
Nov.  13  Physical  Review  Letters,  Wolf 
outlines  a  more  general  -  and  perhaps 
physically  more  plausible  —  mechanism 
that  could  Imitate  Doppler  shifts  of  any 
magnitude.  Instead  of  requiring  the  mi¬ 
crolamps  In  the  source  to  fluctuate  In 
some  correlated  fashion,  he  now  pro¬ 
poses  that  a  complex  “scattering  me¬ 
dium*  such  as  the  electrically  charged 


and  frenetic  atmosphere  thought  to  sur¬ 
round  quasars,  might  serve  as  an  unusual 
lens  that  restructures  Incoming  light  to 
have  redshifting  or  blueshifting  correla¬ 
tions  upon  leaving  the  medium.  ‘A  scat¬ 
tering  medium  of  the  right  type  between 
the  source  and  an  observer  should  pro¬ 
duce  these  effects,"  Wolf  told  Science 
News. 

Wolf  concedes  that  astronomers  have 
never  reported  such  a  scattering  medium 
and  notes  that  he  used  simplifying  as¬ 


sumptions  In  both  the  original  and  up 
dated  theories.  Nonetheless,  the  ex 
pension  ol  the  theory  strengthens  thi 
case  for  a  third  physical  mechanisn 
underlying  spectral  shifts  even  in  light 
from  stationary  sources.  University  of 
Rochester  astronomer  Malcolm  P  Save- 
doff  says  the  soundness  of  Wolf's  theory 
demands  that  scientists  take  it  seriously 
In  a  paper  submitted  to  Astrophysical 
Journal,  Savedoff,  Wolf  and  graduate 
student  Daniel  FV  James  have  mathe¬ 
matically  modeled  scattering  media  that 
are  consistent  with  typical  models  of  the 
environment  near  quasars.  -  /.  Amato 
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Novel  red  shifts 

UNUSUAL  effects  in  which  light  is 
scattered  might  be  causing  large  shifts 
in  the  wavelengths  of  light  in  the  optica! 
spectra  of  distant  galaxies  and  quasars.  This 
could  mimic  the  effects  of  the  expansion  of 
the  Universe,  which  astronomers  believe 
causes  the  light  in  distant  objects  to  be 
greatly  red  shifted.  The  phenomenon  might 
explain  some  puzxling  associations  of  qua¬ 
sars  with  galaxies  which  have  much  smaller 
red  shifts  ( Physical  Review  Letters,  vol  63, 

p  2220). 

Emil  Wolf,  an  optical  theorist  at  the 
University  of  Rochester  in  New  Y ork  State , 
has  carried  out  calculations  together  with 
graduate  student  Daniel  James.  They  find 
that  large  shifts  in  wavelength  are  possible 
when  light  is  scattered  from  media  with 
optical  properties  that  change  randomly  in 
space  and  time.  Such  media  may  occur  in 
tne  vicinity  of  the  “engines”  that  drive 
quasars,  violent  galaxies  that  astronomers 
do  not  understand  well  at  all.  The  shift 
would  be  similar  at  all  wavelengths,  thus 
mimicking  the  shift  caused  by  the  expansion 
of  the  Universe. 

Wolf  first  showed,  a  few  years  ago,  that 
small  shifts  could  occur  in  light  that  is 
partially  coherent,  when  viewed  from  dif¬ 
ferent  angles  ( Physical  Review  Letters,  vol 
58,  p  2646).  Researchers  have  since  seen 
such  changes  in  the  laboratory.  Until  his 
discovery,  physicists  thought  that  such  shifts 
in  wavelength  could  be  caused  only  by  the 
motion  of  the  light  source,  or  by  strong 
gravitational  fields. 

The  large  shifts,  which  Wolf  and  James 
predict,  have  yet  to  be  demonstrated.  So 
tar,  astronomers  are  unconvinced.  “Even  if 
they  can  be  produced  in  the  laboratory ,  they 
may  not  exist  in  nature,”  says  Malcolm 
Savedeff,  an  astronomer  at  the  University 
of  Rochester.  □ 
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Cosmic 
theory 
gets  red 
light 

By  Robert  Matthews 
Science  Correspondent 


RESEARCH  by  a  distinguished 
professor  of  physics  is  likely  to 
start  a  major  controversy  about 
the  origin  of  the  universe. 

The  work,  published  in  an 
Influential  journal  of  physics, 
challenges  one  of  *hp  funda¬ 
mental  tenets  of  cosmology: 
that  light  from  distant  galaxies 
reveals  how  fast  the  universe  is 
expanding. 

According  to  current  think¬ 
ing,  the  universe  began  in  a  cat¬ 
aclysmic  explosion  —  the  Big 
Bang  —  about  15  billion  years 
ago  and  has  been  expanding 
ever  since. 

The  proof  for  this  assertion 
dates  from  the  1920s,  when  the 
American  astronomer  Edwin 
Hubble  discovered  that  galax¬ 
ies  were  racing  away  from  one 
another,  with  their  speed 
increasing  as  ritatanrp 
increased.  By  measuring  this 
rate  of  expansion  of  the  uni¬ 
verse,  it  is  possible  to  work  out 
the  date  of  the  Big  Bang. 

But  Hubble's  observations 
rely  crucially  on  the  so-called 
“red  shift”  effect.  This  is  simi¬ 
lar  to  the  change  in  pitch  heard 
from  a  police  car  siren  as  it 
passes  by.  Light  coming  from  a 
receding  galaxy  is  shifted 
towards  longer,  redder,  wave¬ 
lengths  of  light.  By  measuring 
the  "red  shift”  of  a  galaxy  It  is 
possible  to  estimate  its  speed. 

The  fundamental  assumption 
is  that  only  movement  away 
from  an  observer  can  cause  a 
red  shift.  But  now  Emil  Wolf, 
professor  of  optical  physics  at 
the  university  of  Rochester, 
New  York  State,  and  graduate 
student  Daniel  James  have 
found  another  cause  of  red 
shifts,  which  has  nothing  to  do 
with  movement. 


It  casts  doubt  on  one  of  the 
few  things  cosmologists 
thought  they  knew  for  certain: 
that  red  shifts  can  be  relied  on 
to  work  out  the  age  and  fate  of 
the  universe. 

Professor  Wolf  has  discov¬ 
ered  that  it  is  possible  to  mimic 
a  red  shift  when  light  coming 
from  a  source  bounces  off  mat¬ 
erial  with  certain  optical  prop¬ 
erties  that  change  with  both 
Hm>  and  location  relative  to  the 
source. 

Details  of  the  research  are 
published  in  the  November  13 
Issue  of  Physical  Revine  Letters. 

The  process  can  account  for  a 
red  shift  of  any  size.  Including 
those  associated  with  quasars, 
the  incredibly  bright  centres  of 
galaxies,  thought  to  be  the  most 
distant  objects  in  the  universe. 

Backing  for  Professor  Wolfs 
work  may  have  already  come 
from  a  long-standing  contro¬ 
versy  concerning  these  objects. 

Photographs  taken  with  the 
world's  largest  telescopes  by  Dr 
Halton  Arp,  an  American 
astronomer,  appear  to  show 
physical  links  between  some 
quasars  and  galaxies.  Yet  the 
red  shifts  of  the  two  objects  are 
vastly  different,  suggesting 
they  are  billions  of  light-years 
apart. 

The  standard  explanation  is 
that  the  two  objects  were  lying 
in  the  same  part  of  the  sky,  and 
were  superimposed  on  each 
other. 

Professor  Wolf  believes  his 
new  research  may  mean  that 
Arp’s  contention  that  the  pho¬ 
tographs  prove  that  there  is 
something  wrong  with  the  red 
shift  idea  will  have  to  be  re¬ 
examined.  Others  think  not. 

Professor  Martin  Rees,  a  lead¬ 
ing  astrophysicist  at  the  Insti¬ 
tute  of  Astronomy  in 
Cambridge,  said  there  was 
plenty  of  evidence  supporting 
the  conventional  interpretation 
of  red  shifts. 


England) 


Scrnmnc  American  January  1990 


The  Redshift  Blues 

New  redshift  theory  challenges 
both  physicists  and  cosmologtsts 

Redshifts  are  to  astronomy  what 
tape  measures  are  to  carpentry: 
a  well-understood  tool  whose 
validity  hardly  seems  open  to  ques¬ 
tion.  Modem  cosmology  Is  rooted  in 
the  belief  that  essentially  all  observed 
redshifts  are  caused  by  the  Doppler 
effect,  whereby  light  emitted  by  a  re¬ 
ceding  object  is  “stretched  out”  and 
so  shifted  toward  the  red  end  of  the 
spectrum.  Emil  Wolf  of  the  University 
of  Rochester,  writing  in  Physical  Re¬ 
view  Letters,  now  questions  this  tradi¬ 
tional  Interpretation  by  proposing  a 
previously  unrecognized  mechanism 
that,  he  claims,  can  completely  mimic 
the  Doppler  effect 
Wolf  suggests  that  light  can  be  red- 
shifted  if  It  passes  through  a  scat¬ 
tering  medium  in  which  the  Index  of 
refraction  varies  randomly  over  both 
space  and  time.  If  suitable  correlations 
exist  within  the  medium  the  light  will 
change  frequencies  even  though  the 
overall  source  is  at  rest  His  calcula¬ 
tions  show  that  the  result  can  In  prin¬ 
ciple  be  a  Doppler-like  shift  across  the 
entire  spectrum 

Wolf’s  ideas  have  generated  con¬ 
siderable  interest  and  controversy 
among  his  colleagues,  in  part  because 
the  big-bang  theory  rests  on  obser¬ 
vations  of  redshifts  that  are  Interpret¬ 
ed  as  evidence  that  the  universe  is 
expanding.  Wolf  hesitates  to  suggest 
that  the  blg-bang  concept  might  be 
wrong.  “It  is  something  that  should  be 
looked  at,"  he  says  cautiously. 

More  plausible.  Wolf  believes,  is  the 
possibility  that  his  mechanism  may 
explain  long-disputed  quasar-galaxy 
pairings.  In  some  instances,  quasars 
and  galaxies  that  appear  to  be  physi¬ 
cally  associated  have  vastly  different 
redshifts.  which  indicates,  by  conven¬ 
tional  reckoning,  that  they  are  billions 
of  light-years  apart  Many  astrono¬ 
mers— Including  Christopher  L.  Carilll 
of  Harvard  University,  who  with  two 
colleagues  found  the  most  recent 
such  association— dismiss  these  asso¬ 
ciations  as  chance  line-of-slght  pair¬ 
ings.  A  few  workers,  among  them  Jack 
W.  Sulentic  of  the  University  of  Ala¬ 
bama  at  Tuscaloosa,  maintain  that 
they  reveal  either  anomalies  unex¬ 
plained  by  the  big-bang  theory  or  else 
the  existence  of  some  kind  of  non- 
Dopplerian  redshift 
Naturally,  Wolf  favors  the  latter  ln- 
terpretatloa  Matter  surrounding  qua¬ 
sars  Is  probably  anisotropic  (not  the 


tame  In  all  directions)  because  of  tur¬ 
bulence  or  the  powerful  Jets  of  matter 
that  quasars  often  emit  Such  anisot¬ 
ropy  could  produce  the  sort  of  corre¬ 
lation  mechanism  that  Wolf  has  stud¬ 
ied,  generating  redshifts  that  have 
hitherto  been  attributed  to  the  qua¬ 
sars'  great  distance  from  the  earth.  If 
confirmed,  this  finding  could  solve  an¬ 
other  quasar  mystery:  some  quasar 
jets  appear  to  be  expanding  faster 
than  the  speed  of  light  If  quasars  were 
significantly  less  distant  than  general¬ 
ly  believed,  the  jets  would  be  much 
smaller  and  their  rate  of  expansion 
within  the  cosmic  speed  limit 
Wolf  is  searching  for  more  down-to- 
earth  applications,  such  as  methods 
for  correcting  satellite-tracking  sig¬ 
nals  and  for  Improving  reference  stan¬ 
dards.  He  also  reports  Interest  from 
the  Department  of  Defense,  which  is 
intrigued  by  the  possibility  that  artifi¬ 
cially  Induced  spectra]  modulation 
could  provide  the  ultimate  in  coded 
communications.  — Corey  S.  Powell 
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Invariance  of  the  Spectrum  of  Light  on  Propagation 
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(Received  27  January  1986) 

The  question  is  raised  as  to  whether  the  normalized  spectrum  of  light  remains  unchanged  on 
propagation  through  free  space.  It  is  shown  that  for  sources  of  a  certain  class  that  includes  the  usu¬ 
al  thermal  sources,  the  normalized  spectrum  will,  in  general,  depend  on  the  location  of  the  observa¬ 
tion  point  unless  the  degree  of  spectral  coherence  of  the  light  across  the  source  obeys  a  certain  scal¬ 
ing  law  Possible  implications  of  the  analysis  for  astrophysics  are  mentioned. 


PACS  numbers  42.10  Mg.  07.6S.-b,  42  68. Hf 

Measurements  of  the  spectrum  of  light  are  generally 
made  some  distance  away  from  its  sources  and  in 
many  cases,  as  for  example  in  astronomy,  they  are 
made  exceedingly  far  away.  It  is  taken  for  granted  that 
the  normalized  spectral  distribution  of  the  light  in¬ 
cident  on  a  detector  after  propagation  from  the  source 
through  free  space  is  the  same  as  that  of  the  light  in 
the  source  region.  I  will  refer  to  this  assumption  as 
the  assumption  of  invariance  of  the  spectrum  on  prop¬ 
agation.  This  assumption,  which  is  implicit  in  all  of 
spectroscopy,  does  not  appear  to  have  been  previously 
questioned,  probably  because  with  light  from  tradition¬ 
al  sources  one  has  never  encountered  any  problems 
with  it  However,  with  the  gradual  development  of 
rather  unconventional  light  sources  and  with  the  rela¬ 
tively  frequent  discoveries  of  stellar  objects  of  an  un¬ 
familiar  kind,  it  is  obviously  desirable  to  understand 
whether  all  such  sources  generate  light  whose  spec¬ 
trum  is  invariant  on  propagation,  and  if  so,  what  the 
reasons  for  it  are.  Actually  it  is  not  difficult  to  con¬ 
ceive  of  sources  that  generate  light  whose  spectrum  is 
not  invariant  on  propagation.  In  this  note  I  will  show 
what  are  the  characteristics  of  a  certain  class  of  sources 
that  generate  light  whose  spectrum  is  invariant,  at  least 
in  the  far  zone 

From  the  standpoint  of  optical  coherence  theory,  in¬ 
variance  of  the  spectrum  of  light  on  propagation  from 
conventional  sources  is  a  rather  remarkable  fact,  as  can 
be  seen  from  the  following  simple  argument.  Consid¬ 
er  an  optical  field  generated  by  a  stationary  source  in 
free  space.  The  basic  field  variable,  say  the  electric 
field  strength  at  the  space-time  point  (r,r),  may  be 
represented  by  its  complex  analytic  signal1,2  £(r,f). 
According  to  the  Wiener-Khintchine  theorem3  the 
spectral  density  of  the  light  at  the  point  r  is  then 
represented  by  the  Fourier  transform, 

S(r.w)-  J"_r(r-T  >*""**•  (,) 

of  the  autocorrelation  function  (known  in  the  optica! 
context  as  the  self-coherence  function)  of  the  field 
variable.  It  is  defined  as 

f(r. r)  -  (£*(r,/)£(r,t  +  r)),  (2) 


where  the  angular  brackets  denote  the  ensemble  aver¬ 
age.  Now  the  spectral  density  and  the  self-coherence 
function  are  the  “diagonal  elements”  (r2  —  r,  —  r )  of 
two  basic  optical  correlation  functions,  viz.,  the  cross- 
spectral  density 

(^(ri,  r2,w)  -  J*~^r(r1,r2.  r)e'“T  dr.  (3) 

and  the  mutual  coherence  function 

r(r,,r2,r)-  (£*(r,./)£(r2,r  +  t)).  (4) 

It  is  well  known  that  both  the  mutual  coherence  func¬ 
tion  and  the  cross-spectral  density  obey  precise  prop¬ 
agation  laws.  For  example,  in  free  space4 

(V,2  +  *2)H'(r1,r2.a))-0  0-1.2).  (5) 

where 

k  -  w/c,  (6) 

with  c  being  the  speed  of  light  in  vacuo  and  V/  being 
the  Laplacian  operator  acting  with  respect  to  the  vari¬ 
able  r j.  Consequently,  both  the  mutual  coherence 
function  and  the  cross-spectral  density  and,  in  fact, 
also  their  normalized  values  change  appreciably  on 
propagation.  For  example,  for  a  spatially  incoherent 
planar  source  W''(r1.r2,ci>)  and  r(r|,r2,  r)  will  be 
essentially  6  correlated  with  respect  to  rt  and  r2  at  the 
source  plane  but  will  have  nonzero  values  for  widely 
separated  pairs  of  points  which  are  sufficiently  far  away 
from  the  source.  This  is  the  essence  of  the  well 
known  vsn  Cittert-Zernike  theorem  (Ref.  1,  Sect. 
10.4.2).  In  physical  terms,  the  correlation  in  the  field 
generated  by  a  spatially  incoherent  source  may  be 
shown  to  have  its  origin  in  the  process  of  superposi¬ 
tion.  We  thus  have  the  following  rather  strange  situa¬ 
tion:  The  correlations  of  the  light  may  change  drasti¬ 
cally  on  propagation;  yet,  under  commonly  occurring 
circumstances,  their  (suitably  normalized)  diagonal 
elements,  which  represent  the  spectrum  of  the  light  or 
its  Fourier  transform,  remain  unchanged. 

To  obtain  some  insight  into  this  problem  we  con¬ 
sider  light  generated  by  a  very  simple  model  source; 
namely,  a  planar  source  occupying  a  finite  domain  D  of 
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a  plane  z  -  0  and  radiating  into  the  half  space  z  >  0, 
which  has  the  same  spectral  distribution  S<0’(«)  at 
each  source  point  P(p)  and  whose  degree  of  spectral 
coherences  pt0)(pi.p2, w)  is  statistically  homogene¬ 
ous,  i.e.,  has  the  functional  form  pm(pi~  p\,u>). 
The  cross-spectral  density  of  the  light  across  the 
source  plane  is  then  given  by 

M/'0,(pi,p2.  w) 

-<(pi)«(p2)5(0,(w)^(0,(pj-p|.  w),  (7) 

where  « (p )  —  1  or  0  according  to  whether  the  point 
P(p)  is  located  within  or  outside  the  source  area  D  in 
the  plane  z  -  0. 

We  will  also  assume  that  at  each  effective  frequency 
u  present  in  the  source  spectrum,  the  linear  dimen¬ 
sions  of  the  source  are  much  larger  than  the  spectral 
correlation  length  (the  effective  width  A  of  |p<0)(p\ 
<»)|].  Sources  of  this  kind  belong  to  the  class  of  so- 
called  quasihomogeneous  sources ,6  which  have  been  ex¬ 
tensively  studied  in  coherence  theory  in  recent  years. 
Most  of  the  usual  thermal  sources  are  of  this  kind. 

The  radiant  intensity  y„(u),  i.e.,  the  rate  at  which 
energy  is  radiated  at  frequency  u  per  unit  solid  angle 
around  a  direction  specified  by  a  unit  vector  u,  is  given 
by  the  expression  (cf.  Ref.  6,  Eq.  (4.8)] 

^.(u)  -fc2y4S(0,(cu)p<0,(/cui.w)cosJe.  (8) 

In  this  formula,  A  is  the  area  of  the  source, 

A<0,(f.‘*')-TTT7  JV0,(p'.w>c-'f  V<*V  W 

is  the  two-dimensional  spatial  Fourier  transform  of  the 
degree  of  spectral  coherence,  ux  is  the  transverse  part 
of  the  unit  vector  u,  i.e.,  the  component  of  u  (con¬ 
sidered  as  a  two-dimensional  vector)  perpendicular  to 
the  z  axis,  and  0  is  the  angle  between  the  u  and  the  r 
directions  (see  Fig.  1).  Evidently  the  normalized  spec¬ 
tral  density  S<">(u,w)  at  a  point  in  the  far  zone,  in 
the  direction  specified  by  the  unit  vector  u,  is  given  by 

S(m)(u.w)-Jm(u)/fjm(u)du>.  (10) 


z 


FIG.  I.  Illustration  of  the  notation. 


On  substituting  Eq.  (8)  into  Eq.  (10)  we  obtain  for  the 
normalized  spectrum  in  the  far  zone  the  expression 

*2S<#,U.)M(0'Uui.a») 


^'-'(u.w)- 


fA2S(0>(w)p<0,(kui.w)dw 


(11) 


It  is  clear  from  Eq.  (11)  that  the  normalized  spec¬ 
trum  of  the  light  depends  on  the  direction  u;  i.e.,  it  is 
in  general  not  invariant  throughout  the  far  zone. 
However,  it  is  seen  at  once  from  Eq.  (11)  that  it  will 
be  invariant  throughout  the  far  zone  if  the  Fourier 
transform  of  the  degree  of  spectral  coherence  of  the 
light  in  the  source  plane  is  the  product  of  a  function  of 
frequency  and  a  function  of  direction,  i.e.,  it  is  of  the 
form 


MM,,(KuA.u*)-F(a,)//(u1). 
In  this  case  Eq.  (11)  reduces  to 


Stm)(  u,w) 


*2.S(0'(a,)fU) 

JVs,0,(a>)FU)rfa/ 


(12) 


(13) 


and  the  expression  on  the  right  is  independent  of  the 
direction  u. 

I  will  now  show  that  the  condition  (12)  has  some  in¬ 
teresting  implications,  which  follow  from  the  fact  that 
ft*0’  is  a  correlation  coefficient.  Before  doing  this  we 
note  that  since  u  is  a  unit  vector,  luj  <  1.  However, 
we  will  now  assume  that  the  factorization  condition 
(12)  holds  for  all  two-dimensional  vectors  (0 
«  |ux|  <  oo).  This  assumption  will  be  trivially  satis¬ 
fied  if  the  degree  of  spectral  coherence  p,0,(p\  w)  is. 
at  each  effective  temporal  frequency  u>,  band  limited 
in  the  spatial  frequency  plane  to  a  circle  of  radius  k 
about  the  origin;  in  more  physical  terms  this  condition 
means  that  p <0,(p\«)  does  not  vary  appreciably  over 
distances  of  the  order  of  the  wavelength  k^lnc/ut. 
With  this  being  understood  let  us  take  the  Fourier 
transform  of  Eq.  (12).  We  then  find  at  once  that 

M<0,(p\«> 

-  F  (to)J*//(ux)exp(ifcux  •  p')  d2(k\ii).  (14) 

i.e., 

p<0)(p.to)-k2F(<o)f/(kp').  (15) 

where  H  is,  of  course,  the  two-dimensional  Fourier 
transform  of  H.  Since  p,0>(p‘.  «)  is  a  correlation  co¬ 
efficient  it  has  the  value  unity  when  p'-0,  i.e., 

p<0)(0,  w )  —  I,  for  all  <u.  (16) 

and  hence  Eq.  (15)  implies  that 

k2F(u)  —  I  //  ( 0)  1  ~ 1 .  (17) 

Since  the  left-hand  side  of  Eq.  (17)  depends  on  the 
frequency  but  the  right-hand  side  is  independent  of  it. 
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each  side  must  be  a  constant  (a  say)  and  consequently 


F(w)~alk\ 


(18) 


Two  important  conclusions  follow  at  once  from 
these  results.  If  we  substitute  Eq.  (18)  into  Eq.  (13) 
we  obtain  the  following  expression  for  the  normalized 
spectrum  of  light  in  the  far  zone: 


S'-'tu.oJ-S'-’fw) 


S(0i(a») 

J"s,0,(«u)do) 


(19) 


This  formula  shows  that  not  only  is  the  normalized 
spectrum  of  the  light  now  the  same  throughout  the  far 
zone,  but  it  is  also  equal  to  the  normalized  spectrum  of 
the  light  at  each  source  point. 

Next  we  substitute  Eq.  (18)  into  Eq.  (15)  and  set 
<*//«/;,  p ** p2 —  p i ■  We  then  obtain  for  p(0>  the  ex¬ 
pression 

p<0>(p2~  pi.w)  -  h  [k  (p2  —  Pi )) 

U««/c);  (20) 

i.e.,  the  complex  degree  of  spectral  coherence  is  a 
function  of  the  variable  A,(py-  pt)  only.  We  will 
refer  to  Eq.  (20)  as  the  scaling  lave.  Obviously  for  a 
source  that  satisfies  this  law,  the  knowledge  of  the  de¬ 
gree  of  spectral  coherence  of  the  light  in  the  source 
plane  at  any  particular  frequency  w  specifies  it  for  all 
frequencies. 

The  scaling  law  (20),  which  ensures  that  for  sources 
of  the  class  that  we  are  considering  the  normalized 
spectrum  of  the  light  is  the  same  throughout  the  far 
zone  and  is  equal  to  the  normalized  spectrum  of  the 
light  at  each  source  point  (Eq.  (19)),  is  the  main  result 
of  this  note. 

It  is  natural  to  inquire  whether  sources  are  known 
that  obey  this  scaling  law.  The  answer  is  affirmative. 
Many  of  the  commonly  occurring  sources,  including 
blackbody  sources,  obey  Lambert's  radiation  law  (Ref. 
1,  Sect.  4.8.1).  It  is  known7  that  all  quasi-homo- 
geneous  Lambertian  sources  have  the  same  degree  of 
spectral  coherence,  viz. 


<z<0,(p2-pi.  w)  -sin(k  |p2-pil)/k  lp2-pil,  (21) 


which  is  seen  to  satisfy  the  scaling  law  (20).  Accord¬ 
ing  to  the  preceding  analysis  such  sources  will  generate 
light  whose  normalized  spectrum  is  the  same 
throughout  the  far  zone  and  is  equal  to  the  normalized 
spectrum  at  each  source  point.  This  fact  is  undoubted¬ 


ly  largely  responsible  for  the  commonly  held,  but 
nevertheless  incorrect,  belief  that  spectral  invariance  is 
a  general  property  of  light. 

This  Letter  has  dealt  with  what  is  probably  the  sim¬ 
plest  problem  regarding  spectral  invariance  on  prop¬ 
agation.  It  would  seem  that  some  significant  questions 
in  this  area  might  be  profitably  studied.  Among  them 
are  the  elucidation  of  the  physical  origin  of  the  scaling 
law,  spectral  properties  of  light  from  a  broader  class  of 
sources  than  considered  here,  the  relation  between  the 
scaling  law  and  Mandel's  results  regarding  cross- 
spectrally  pure  light,*- 9  and  relativistic  effects.  Appli¬ 
cations  of  the  results  to  problems  of  astrophysics 
might  be  of  particular  interest;  at  this  stage  one  might 
only  speculate  whether  source  correlations  may 
perhaps  not  give  rise  to  differences  between  the  spec¬ 
trum  of  the  emitted  light  and  the  spectrum  of  the 
detected  light  that  originates  in  some  stellar  sources. 

It  is  a  pleasure  to  acknowledge  stimulating  discus¬ 
sions  with  Professor  Leonard  Mandel  about  the  subject 
matter  of  this  note.  This  research  was  supported  by 
the  National  Science  Foundation  and  by  the  Air  Force 
Geophysics  Laboratory  under  Air  Force  Office  of 
Scientific  Research  Task  No.  2310G1. 


‘•’Also  at  the  Institute  of  Optics,  University  of  Rochester, 
Rochester,  N.  Y.  14627. 

■M.  Born  and  E.  Wolf,  Principles  of  Optics  (Pergamon,  Ox¬ 
ford  and  New  York,  1980),  6th  ed..  Sect  10  2. 

JL.  Mandel  and  E.  Wolf.  Rev.  Mod.  Phys.  37,  231  (1965). 

3C  Kiltel,  Elementary  Statistical  Fhysics  (Wiley,  New  York, 
1958),  Sect.  28. 

4E.  Wolf,  J.  Opt.  Soc.  Am.  68.  6  (1978),  Eqs.  (5.3). 

5The  degree  of  spectral  coherence  is  defined  by  the  for¬ 
mula  Icf.  L.  Mandel  and  E.  Wolf,  J.  Opt.  Soc.  Am  66,  529 
(1976)) 

,o,r  j  _ W'°'(p>.p:.u>) _ 

M  lH''0,(p).p1.<v)),/I(W,,0,(p3.p2.a,))1/J 

*W.  H.  Carter  and  E.  Wolf,  J.  Opt.  Soc.  Am.  67,  785 
(1977). 

7W.  H.  Carter  and  E.  Wolf,  J.  Opt.  Soc.  Am.  65.  1067 
(1975). 

*L.  Mandel,  J.  Opt  Soc  Am.  51,  1342  (1961). 

’See,  Mandel  and  Wolf,  Ref.  5. 


1372 


Reprinted  from  Nature,  Vol.  326,  No.  6111,  pp.  363-365,  26  March  1987 
■c  Macmillan  Journals  Lid.,  1987 


4.2 


Non-cosmological  redshifts 
of  spectral  lines 

Emil  Wolf 
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We  allowed  in  a  recent  report'  (see  also  refs  2-4)  that  the  normal¬ 
ized  spectrum  of  light  will,  in  general,  change  on  propagation  in 
free  space.  We  also  showed  that  the  normalized  spectrum  of  light 
emitted  by  a  source  of  a  well-defined  class  will,  however,  be  the 
same  throughout  the  far  zone  if  the  degree  of  spectral  coherence 
of  the  source  satisfies  a  certain  scaling  law.  The  usual  thermal 
sources  appear  to  be  of  this  kind.  These  theoretical  predictions 
were  subsequently  verified  by  experiments’.  Here,  we  demonstrate 
that  under  certain  circumstances  the  modification  of  the  normal¬ 
ized  spectrum  of  the  emitted  light  caused  by  the  correlations 
between  the  source  fluctuations  within  the  source  region  can  pro¬ 
duce  redshifts  of  spectral  lines  in  the  emitted  light.  Our  results 
suggest  a  possible  explanation  of  various  puzzling  features  of  the 
spectra  of  some  stellar  objects,  particularly  quasars. 

To  explain  why  source  correlations  influence  the  spectrum  of 
the  emitted  light  consider  a  very  simple  example.  Suppose  that 
two  point  sources  P,  and  P2  have  identical  spectra  SQ(w)  and 
that  measurements  on  the  emitted  field  are  made  at  some  point 
P.  The  sources  are  assumed  at  rest  relative  to  an  observer  at  P. 
Assuming  that  the  source  fluctuations  can  be  described  by  a 
stationary  ensemble,  the  field  at  P  may  be  characterized  by  an 
ensemble  { V( P,  w)}  of  frequency-dependent  realizations6,  each 
of  the  form 

ikR,  -.kR, 

V(P,«u)=Q(P„a,)~— K?(P,,<u)  —  (1) 

K,  R2 

where  {QiP^u)),  (;'=  1,2),  characterize  the  strengths  of  the 
two  fluctuating  point  sources,  R ,  and  R2  are  the  distances  from 
P,  to  P  and  from  P2  to  P  respectively  (see  Fig  1)  and  It  =  w/c, 
c  being  the  speed  of  light  in  vacuo.  For  simplicity  polarization 
effects  are  ignored  and  hence  V  and  Q  are  taken  to  be  scalars. 
The  spectrum  of  the  light  at  P  is  then  given  by 

S\  (P,  to)  =  <  V*(P,  w)V(P,  w))  (2) 

where  the  asterisk  denotes  the  complex  conjugate  and  the 
angular  brackets  denote  the  ensemble  average.  On  substituting 
from  equation  (1)  into  equation  (2)  and  using  the  fact  that 

(Q0(P„u>)Q(P„w))  =  (Qt(P2.w)Q(P2.u>))  =  SoM  (3) 

the  following  expression  is  obtained  for  the  spectrum  of  the 
emitted  light  at  P: 

SV(P,  ->)  =  + 

[  .UR.-R,)  -) 

W0(P„PI,m)-  —  +C.C.J  (4) 

Here 

Wq(P„P2,w)  =  (Q*(Puu,)Q(P2,w))  (5) 

is  the  so-called  cTOss-spectral  density  of  the  source  fluctuations 
and  c.c.  denotes  the  complex  conjugate. 

The  formula  (4)  shows  that  the  spectrum  SV(P,  <v)  is,  in 
general,  not  just  proportional  to  S0{u>)  but  is  modified  by  the 
correlation,  characterized  by  W0(  P,,  P2l  <d),  between  the  fluctu¬ 
ations  of  the  two  source  strengths  Q(P,,  w)  and  Q(  P2,  o»).  Only 
in  some  very  special  cases,  for  example,  when  the  source  fluctu¬ 
ations  are  uncorrelated  [  W0(P,,  P2,  w)  =  0]  will  Sy(P,u>)  be 
proportional  to  S9(<u).  Hence,  in  general,  the  spectrum  of  the 
light  generated  by  two  point  sources  depends  not  only  on  their 
spectra  but  also  on  the  correlation  between  the  fluctuations  of 


their  strengths. 

A  generalization  of  the  elementary  formula  (4)  for  radiation 
from  three-dimensional  steady-state  (that  is,  statistically  station¬ 
ary)  sources  of  any  state  of  coherence  is  known’.  Of  special 
interest  in  the  present  context  is  the  form  that  the  formula  takes 
when  the  source  has  the  same  normalized  spectrum  SqIu), 
(ir  *<?(")  dw  =  1)  at  each  point  in  the  source  region  and  has  a 
degree  of  spectral  coherence5  (appropriately  normalized  cross- 
spectral  density)  ^Q(t„  tj,  *>)  that  depends  on  the  position  vec¬ 
tors  r,  and  r2  of  any  source  points  P,  and  P2  only  through  their 
difference  r2-r,.  If,  in  addition,  for  each  frequency  that  sig¬ 
nificantly  contributes  to  the  source  spectrum,  the  spectral  corre¬ 
lation  length  [the  effective  spatial  width  |Arj  of  |*i(r',  w)!]  is 
small  compared  to  the  linear  dimensions  of  the  source,  the 
normalized  spectrum  s'^^u,  <o)  of  the  emitted  light  in  the  far 
zone,  in  a  direction  specified  by  a  unit  vector  u,  becomes  (see 
equation  (3.11)  of  ref.  8) 


s'f'(  u,w)« 


S0(ttl)/IQ(fcu,  cd) 

\  s0(u))/I0(ku,  w)  dw 


where  ft0(K,u)  is  the  three-dimensional  spatial  Fourier  trans¬ 
form  of  the  degree  of  spectral  coherence  /i0(r',ai)^ 
Mo(r2-r,,w). 

Let  us  now  choose  as  the  normalized  source  spectrum  sqUu  i 
a  spectral  line  with  a  gaussian  profile 


s©(w)  =-£-^==exp  [-(<u  —  <u°)I/26‘]  (6«  wc) 


and  suppose  that  at  each  effective  frequency  w,  the  source 
correlation  decreases  with  the  separation  |rj  =  |r2  —  r,|  of  any  two 
source  points  in  a  gaussian  manner,  that  is 

/i0(r’,  <i>)  =  exp[-r'5/2<r5(w)]  (81 

On  taking  the  Fourier  transform  of  equation  (8)  and  substituting 
the  resulting  expression  into  equation  (6)  we  obtain  the  follow¬ 
ing  expression  for  the  normalized  spectrum  of  the  emitted  light 
in  the  far  zone  (see  equation  (3.21)  of  ref.  8) 


sn«v)= 


z<?(iu)o-5(«u)exp{-j[lc(r).(iu)]:} 

Jo  s0(u»)cr5  (<u)  exp  {-4[kow(«u)]:}  du> 


Here,  s'v‘(<v)  is  written  in  place  of  s‘'’(u,  w),  because  the 
spectrum  of  the  far  field  is  now  independent  of  u,  as  a  con¬ 
sequence  of  the  assumed  isotropy  of  (see  equation  (8)1. 

The  formula  (9)  shows  that  the  spectrum  of  the  emitted  light 
in  the  far  zone  depends  both  on  the  spectrum  of  the  source 
fluctuations  and  on  the  manner  in  which  the  effective  source 
correlation  length  crh(w)  depends  on  the  frequency  w. 

Let  us  consider  two  particular  cases.  (1)  Suppose  hist  that 
ct„(w)  is  independent  of  w.  Letting  {  denote  the  (now  constant) 
value  of  o-„  and  with  ^(oi)  given  by  equation  (7),  one  can 
readily  evaluate  the  integral  in  the  denominator  on  the  right  of 
equation  (9)  and  one  then  finds  that 


'=1+7 


When  the  source  is  effectively  spatially  incoherent,  {  -» 0  Then 
according  to  equation  (11)  A-*oc  and  o-*l  and  it  follows  from 
equations  (10)  and  (7)  that  in  this  case 

jV^'Iw)  -•  Jp(«t»)  (1-1 

Hence,  in  the  limiting  case  of  a  completely  incoherent  source 
of  the  class  that  is  considered  here,  the  normalized  spectrum  of 
the  emitted  light  in  the  far  zone  is  identical  with  the  normalized 


p 


1.0 


-C*o 

(incoherent  limit) 


Fig.  1  Illustrating  the  notation  relating  to  derivation  of  the 
formula  (4). 
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spectrum  of  the  source  fluctuations. 

However,  when  the  source  fluctuations  are  correlated  over  an 
effective  distance  (  >0,  equation  (10)  shows  that  the  spectrum 
St* '(«*),  although  it  is  also  a  line  with  a  gaussian  profile,  is 
centred  at  a  lower  frequency  u'a^u0/a2<u>0-  Hence  the  source 
correlations  give  rise  to  a  spectral  line  s'* ’(or)  that  is  redshifted 
with  respect  to  the  spectral  line  produced  by  the  completely 
spatially  incoherent  source  with  the  source  spectrum  s0(w).  The 
shifted  line  is  narrower,  having  root-mean-square  width  6'  = 
6/a  <6  and  has  a-times  greater  height.  Examples  of  spectra 
of  light  in  the  far  zone,  produced  by  several  sources  which  emit 
the  same  spectral  line  but  which  have  different  correlation 
lengths  are  shown  in  Fig.  2.  From  the  formula  (10)  one  can 
readily  deduce  that  the  relative  shift  of  the  line,  namely. 


A0-  A  0  _  <Uq—Wq 

A0  <Uo 

( A0  =  2 Trc/u0,  A o  =  2 7 tc/u'0)  is  given  by 


(13) 


(14) 


which  shows  that  in  this  case  the  redshift  increases  quadratically 
with  the  spectral  source-correlation  length  (.  (2)  Next  consider 
the  situation  when  trM(<u)  =  a/w  where  a  is  a  positive  constant. 
The  expression  (9)  for  the  normalized  spectrum  of  the  emitted 
light  in  the  far  zone  now  reduces  to 


s'v"’(«) 


Sq(u)/u3 

io  [ *Q< «< )/<*»’]  do) 


(15) 


Note  that  this  expression  is  independent  of  the  value  of  the 
constant  a. 

When  sQ( w  )  is  a  line  with  a  gaussian  profile,  given  by  equation 
(7),  the  spectrum  Sy'c,(w),  given  by  equation  (15)  is  no  longer 
strictly  gaussian  but  it  can  be  closely  approximated  by  a  gaussian 
and  can  be  shown  to  be  redshifted  with  respect  to  sQ(w)  by  the 
relative  amount 


2  •» 


(16) 


An  example  of  this  situation  is  illustrated  in  Fig.  3. 

This  case  [o„(<i/)  *  o/w]  is  of  special  interest  because,  accord¬ 
ing  to  equation  (8),  the  degree  of  spectral  coherence  is  now 
given  by 


M<?(r', «) «  exp  [ — ( kr’)2/  2(  o/  c)*],  (17) 

that  is,  it  has  the  functional  form 


w)-/(kr)  (k  -  w/c«2ir/A)  (18) 

Thus  the  degree  of  spectral  coherence  of  the  source  distribution 
now  satisfies  the  three-dimensional  analogue  of  a  requirement 
(called  the  scaling  law)  derived  in  ref.  1,  as  a  sufficient  condition 


Fig.  2  Spectra  '(<*>)  of  the  far  field  from  sources  with  spectrum 
iQ(a>)“(6v/5»r1exp[-(<D-«o)V2AJ]  and  degree  of  spectral 
coherence  s»c>(r’, «)«exp(-r'J/2fJ),  with  <v0“ 3-887 x  10” s'1 
(A0»4,86I  A)  and  8«9.57x  101Js~\  for  several  selected  values 
of  the  effective  source-correlation  length  {.  The  solid  curve  ( { -♦  0 1 
also  represents  the  source  spectrum  j0(u»). 


for  the  spectrum  of  the  light  emitted  by  a  planar  secondary 
source  of  a  well-defined  class  to  have  certain  invariance 
properties  on  propagation.  It  will  be  shown  in  another  publica¬ 
tion  (J.  T.  Foley  and  E.  Wolf,  in  preparation)  that  for  three- 
dimensional  primary  sources  of  an  analogous  class,  whose 
degree  of  coherence  satisfies  this  law,  the  spectrum  of  the  emitted 
light  has  similar  invariance  properties.  We  conjecture  that  the 
usual  thermal  sources  obey  such  a  scaling  law. 

Now  briefly  consider  the  question  of  a  physical  mechanism 
for  producing  source  correlations.  Such  correlations  must  clearly 
be  manifestations  of  some  cooperative  phenomena.  At  the 
atomic  level  possible  candidates  may  perhaps  be  superradiance 
and  superfluorescence9.  An  effect  of  this  kind  was  first  predicted 
by  Dicke  in  1954  when  he  showed10  that  under  certain  circum¬ 
stances  energy  from  excited  atoms  may  be  released  cooperatively 
in  a  much  shorter  time  than  the  natural  lifetime  of  the  excited 
states  of  the  atoms  and  with  much  larger  emission  intensity  than 
would  be  obtained  were  the  atoms  radiating  independently. 


Fig.  3  The  spectrum  s,*,(i»)  of  the  far  field  from  a  source  with 
source  spectrum  Splw) « (fivTrr)"1  exp  [-(w  -  and 

degree  of  spectral  coherence  >*o(r ’•  •*) “ **P (-(Ar "V 72(a/r) •) 
(a«an  arbitrary  constant),  with  •>„»  3.887  x  10”  s"'  (A0  = 
4,861  A)  and  8  «9.57x  io”  s*1.  The  source  spectrum  Sq<w)  is 
shown  for  comparison.  Note  that  •>)  now  obeys  the  scaling 
law. 
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However  not  enough  is  known  at  present  about  the  coherence 
properties  of  large  three-dimensional  systems  of  this  kind  to 
make  it  possible  to  determine  whether  superradiance  and  super¬ 
fluorescence  might  involve  correlations  that  could  give  rise  to 
spectral  line  shifts. 

There  is,  however,  quite  a  different  mechanism,  which  can  be 
described  at  the  macroscopic  level,  and  which  can  imitate  effects 
of  source  correlations;  namely  effects  of  correlations  between 
the  refractive  index  at  pairs  of  points  in  a  spatially  random  but 
statistically  homogeneous,  time-invariant  medium.  If  a  wave 
illuminates  such  a  medium,  say  a  dilute  gas,  then,  as  is  well 
known,  the  medium  acts  as  a  secondary  source,  namely  as  a  set 
of  oscillating  charges  set  in  motion  by  the  incident  wave.  The 
secondary  waves  produced  by  the  oscillating  charges  then  com¬ 
bine  with  eact-  -  and  with  the  incident  wave  and  generate 
the  scattered  if  the  gas  is  not  too  dilute  the  collective 
response  of  the  microscopic  charges  to  the  incident  field  can  be 
described  by  macroscopic  parameters  such  as  the  dielectric 
susceptibility  or  the  refractive  index.  Now  within  the  accuracy 
of  the  first  Bom  approximation  the  basic  equation  for  scattering 
is  of  the  same  form  as  the  basic  equation  for  radiation  from 
primary  sources,  the  ‘equivalent  source'  for  scattering  being  the 
product  of  the  scattering  potential  (which  is  a  simple  function 
of  the  refractive  index )  and  of  the  incident  wave.  This  correspon¬ 
dence  clearly  implies  that  our  results  regarding  the  effects  of 
source  correlations  on  the  spectrum  of  the  emitted  light  must 
have  analogues  regarding  the  effects  of  a  spatially  random 
medium  with  correlated  refractive  index  distribution  on  the 
spectrum  of  the  light  that  is  scattered  by  it.  This  topic  will  be 
discussed  elsewhere. 

Let  us  now  consider  some  implications  of  this  analysis.  Using 
equation  (14),  the  spectral  line  in  Fig  2,  produced  by  the  source 
whose  correlation  length  (  =  A0  is  readily  found  to  have  a  redshift 
given  by  z  =  0.0241  with  respect  to  the  source  spectrum.  It  is  of 
interest  to  note  that  if  an  observer  detected  such  a  redshift 
unaware  of  its  true  origin  and  interpreted  it  on  the  basis  of  the 
Doppler  shift  formula  t/c  =  AA/A0  =  z  he  would  incorrectly 
conclude  that  the  source  was  receeding  from  him  with  a  speed 
v  =  0.024^  =  7,230100  s'1. 
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It  seems  worthwhile  to  note  that  there  is  a  maximum  line  shift 
that  can  be  produced  by  source  correlations  This  can  be  seen 
from  the  basic  formula  (6)  which  indicates  that  s'v  '(u,  w)  =  0 
when  Sq(o»)  =  0,  implying  that  the  spectrum  of  the  far  field  can 
only  contain  those  frequencies  that  are  already  present  in  the 
source  spectrum.  Consequently  the  maximum  attainable 
frequency  shift  of  the  line  cannot  exceed  its  effective  frequency 
range.  However,  any  frequency  contribution  from  the  source 
spectrum  to  the  normalized  spectrum  of  the  far  field  can  be 
greatly  magnified  or  greatly  reduced,  as  is  evident  from  equation 
(6)  and  from  Fig.  2. 

We  have  mainly  considered  effects  of  source  correlations 
under  circumstances  when  the  source  spectrum  consists  of  a 
single  line  and  when  the  degree  of  spectral  coherence  n0  that 
characterizes  the  source  correlations  depends  on  a  single  para¬ 
meter.  Preliminary  calculations  show  that  with  a  suitable  chosen 
Hq  which  depends  on  a  larger  number  of  parameters,  redshifts 
of  several  lines  may  be  produced,  all  of  which  will  have  approxi¬ 
mately  the  same  z-values. 

In  this  article  we  have  considered  redshifts  of  spectral  lines. 
However,  it  is  not  difficult  to  specify  source  correlations  which 
will  produce  blueshifts.  Examples  of  this  kind  are  given  in  a 
forthcoming  publication". 

It  seems  plausible  that  the  mechanism  discussed  in  this  article 
may  be  responsible  for  some  of  the  so  far  unexplained  features 
of  quasar  spectra,  including  line  asymmetries  and  small  differ¬ 
ences  in  the  observed  redshifts  of  different  lines.  In  this  connec¬ 
tion  it  is  of  interest  to  recall  that  the  role  of  coherence  in  the 
emission  of  radiation  from  quasars  was  stressed  by  Hoyle. 
Burbidge  and  Sargent  in  a  well-known  article11. 

I  thank  Mr  A.  Gamliel  and  Mr  K.  Kim  for  carrying  out 
computations  relating  to  the  analysis  presented  in  this  article 
The  fact  that  scattering  can  also  produce  shifts  of  spectral  lines 
was  noted  independently  by  Professor  Franco  Gori,  who  in¬ 
formed  me  of  this  result  when  commenting  on  an  early  version 
of  the  manuscript  of  this  article.  This  investigation  was 
supported  by  the  NSF  and  by  the  US  Air  Force  Geophysical 
Laboratory. 


I  Cantr.  W  H  a  Wolf.  E  Optica  Ana  a,  295-589  (19811 

9  Scb'junnam.  M  F  H  ,  Vrthen.  0  H  F,  Rolder.  D  A  Gibb«.  H  M  in  Adioncts  m  4mmi. 
mud  Moiemlmr  Physics  Vol  17  <eds  lam.  D  A  Bederwn.  B  i  16*- 228  (Academic.  Sew 
York.  1981) 

10  Dicke,  R  H  Phyt  An.  M.  99-1 10  <  1954 1 

11  Wolf.  E  Optics  Commun  (in  the  pmu 

12  Hoyle,  F..  Burbidge.  G  R.  A  Sargem,  W  L  W  Natvrr  lit.  751 .75?  ( 1966 


Volume  58,  Number  25 


PHYSICAL  REVIEW  LETTERS 


22  June  1987 


Red  Shifts  and  Blue  Shifts  of  Spectral  Lines  Emitted  by  Two  Correlated  Sources 

Emil  Wolf(,) 

Department  of  Physics  and  Astronomy,  l/niversity  of  Rochester.  Rochester,  New  York  14627 

(Received  24  February  1987) 

It  has  recently  been  shown  theoretically  that  correlations  between  fluctuations  of  the  source  distribu¬ 
tion  at  different  source  points  can  produce  red  shifts  or  blue  shifts  of  emitted  spectral  lines.  To  facilitate 
experimental  demonstration  of  this  effect  a  simple  example  is  analyzed.  It  involves  only  two  small  ap¬ 
propriately  correlated  sources  and  brings  out  the  essential  physical  features  of  this  new  phenomenon. 

PACS  numbers  42.68.Hf,  07.6S.-b,  42.IO.Mg 

I  showed  not  long  ago  that  the  spectrum  of  light  pro¬ 
duced  by  a  fluctuating  source  depends  not  only  on  the 
source  spectrum  but  also  on  the  correlation  that  may  ex¬ 
ist  between  the  source  fluctuations  at  different  points 
within  the  domain  occupied  by  the  source. 1  This  result 
was  recently  confirmed  experimentally.2  I  also  showed 
that  under  certain  circumstances  source  correlations  may 
produce  red  shifts  or  blue  shifts  of  spectral  lines  in  the 
emitted  radiation.3,4  This  prediction  has  obviously  im¬ 
portant  implications,  particularly  for  astronomy,  and  it  is 
therefore  desirable  to  verify  it  also  by  experiment. 

In  this  Letter  I  analyze  theoretically  one  of  the  sim¬ 
plest  systems  that  will  generate  spectral  shifts  by  this 
mechanism;  namely,  two  small  correlated  sources,  with 
identical  spectra  consisting  of  a  single  line  of  Gaussian 
profile.  I  show  that  with  an  appropriate  choice  of  the 
correlation,  the  spectrum  of  the  emitted  radiation  will 
also  consist  of  a  single  line  with  a  Gaussian  profile,  iicw- 
ever,  this  emitted  line  will  be  red  shifted  or  blue  shifted 
with  respect  to  the  spectral  line  that  would  be  produced 
if  the  sources  were  uncorrelated,  the  nature  of  the  shift 
depending  on  the  choice  of  one  of  the  parameters  that 
specifies  the  exact  form  of  the  correlation  coefficient. 


(/>,  co )  -  ( 1 //?  f  +  1 //?^  )S0  (w )  +  [IVC  (/> ,  ,/>j,  a.  )e'* 

Here 

SQ{a))-(Q*(Ph(i>)Q(P,,(o))-{Q*(P2,(o)Q(P2>(o))  (4) 

is  the  spectrum  (assumed  to  be  the  same)  of  each  of  the  two  source  distributions, 

WQ{PuPhw)-{Q*(Puw)Q{P2,<»))  (5) 

is  the  cross-spectral  density  of  the  source  fluctuations  [first  paper  of  Ref.  6,  Eqs.  (3.3)  and  (5.9)1,  and  c.c.  denotes  the 
complex  conjugate. 

The  degree  of  spectral  coherence  at  frequency  c o,  which  is  a  measure  of  correlation  that  may  exist  between  the  two 
fluctuating  sources,  is  given  by  the  formula® 

PQ(Pt,P2,a>)-WQ(P{%P2,<o)/SQ(o>).  (6) 

The  normalization  in  Eq.  (6)  ensures  that  0<  \pQ(.P\.Pi,to)  i  S  1.  The  extreme  value  |pc|  -1  characterizes  com¬ 
plete  correlation  (complete  spaticl  coherence)  at  frequency  <u.  The  other  extreme  value,  p  "0,  characterizes  complete 
absence  of  correlations  (complete  spatial  incoherence). 

On  substituting  for  Wq  from  Eq.  (6)  into  Eq.  (3),  we  find  that 

Si'(P,w)  “Sp (to Ml //??-♦■  l/7?f  +  [pp(a>)ert<*’*'V/?i/?j+c.c.l),  (7) 


The  main  features  of  this  theoretical  prediction  have 
been  confirmed  by  Bocko,  Douglass,  and  Knox,  using 
acoustical  rather  than  optical  sources.  An  account  of 
their  experiments  is  given  in  the  accompanying  Letter.5 

Let  us  consider  two  small  fluctuating  sources  located 
at  points  P\  and  Pi-  I  assume  that  the  fluctuations  are 
statistically  stationary.  Let  [Q(Pi,to)}  and  {Q(P2,(o)} 
be  the  ensembles  that  represent  the  source  fluctuations6 
at  frequency  <u.  Furthermore,  let  (l/(/>,<o)l  be  the  en¬ 
semble  that  represents  the  field  at  point  P  generated  by 
the  two  sources  (Fig.  1).  Each  realization  U(P, to)  may 
then  be  expressed  in  the  form7 

IkR,  ikR , 

U(P,to)  ^Q(.P\,to)— —  +  g(/>2,&))-£— — ,  (1) 

R)  R : 

where  R i  and  Ri  are  the  distances  from  P\  to  P  and 
from  Pi  to  P,  respectively,  and  k  «tu/c,  c  being  the  speed 
of  hgnt  in  free  space.  The  spectrum  of  the  field  at  the 
point  P  is  given  by 

Su(P,to)  miU*  (P,to)U(P,o)),  (2) 

where  the  angular  brackets  denote  ensemble  average. 
On  substitution  from  Eq.  (1)  into  Eq.  (2),  we  find  that 

*,,/R,R2+c.cl  (3) 
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FIG.  1.  Geometry  and  notation  relating  to  the  determina¬ 
tion  of  the  spectrum  Su(P,u)  of  the  field  at  P  produced  by  two 
small  sources  with  identical  spectra  Sq(to)  located  at  P\  and 

Pi 


where  I  have  omitted  the  arguments  P\  and  P3  in  ji.  For 
the  sake  of  simplicity,  let  us  choose  the  field  point  P  to 
lie  on  the  perpendicular  bisector  of  the  line  joining  /*, 
and  Pi.  Then  Ry^-Ri  (■■/?,  say)  and  formula  (7) 
reduces  to 

Su(P,to)  m,l2/R2)Sg(io)l\  +Re*c(a>)],  (8) 

where  Re  denotes  the  real  part. 

We  note  in  passing  that  when  either  ftQ(m)sO  (mutu¬ 
ally  completely  uncorrelated  sources)  or  when  nq  (<y)s  1 
(mutually  completely  correlated  sources),  the  spectrum 
Su(.P,a> )  of  the  field  at  the  point  P  will  be  proportional 
to  the  spectrum  Sq(.<o)  of  the  source  fluctuations.  How¬ 
ever,  in  general  this  will  not  be  the  case.  In  fact,  it  is 
clear  from  formula  (8)  that  the  field  spectrum  may  differ 
drastically  from  the  source  spectrum,  the  difference  de¬ 
pending  on  the  behavior  of  the  correlation  coefficient 
l iq(o))  as  a  function  of  frequency. 

Suppose  now  that  the  spectrum  of  each  of  the  two 
sources  consists  of  a  single  line  of  the  same  Gaussian 
profile, 


s 


(9) 


where  A,  too,  and  $o  («wo)  are  positive  constants.  Sup¬ 
pose  further  that  the  correlation  between  the  two  sources 
is  characterized  by  the  degree  of  spectral  coherence 

Mq(<o)-af-(m-*),/Uf- 1.  (10) 

where  a,  at,  and  S  ( «  a>i )  are  also  positive  constants.  In 
order  that  expression  (10)  is  a  degree  of  spectral  coher¬ 
ence,  I  must  also  demand  that  o  <  2.  On  substituting 
from  Eqs.  (9)  and  (10)  into  Eq.  (8),  I  obtain  the  follow¬ 
ing  expression  for  the  spectrum  of  the  field  at  the  point 
P: 


Sv(P,a>) 


2 Aa  -(w-mt)l/2ti[ 

R 2 


(11) 


FIG.  2.  Red  shifts  and  blue  shifts  of  spectral  lines  as  pre¬ 
dicted  by  formula  (12).  The  spectrum  Sq(u)  of  each  of  the 
two  source  distributions  is  a  line  with  a  Gaussian  profile  given 
by  Eq.  (9)  with  A- 1,  «o-4.32201  *  10,s  sec-'  (Hg  line 
X  "4358.33  A),  *-5*10*  sec-1,  (a)  The  field  spectrum 
SyiP.ot)  at  P  when  the  two  sources  are  uncorrelated  (jigSO). 
(b),(c)  The  field  spectra  at  P  when  the  two  sources  are  corre¬ 
lated  in  accordance  with  Eq.  (10),  with  o"l.8,  di“7.5*10’ 
sec-1,  and  (b)  #>i“«o-25o  (red-shifted  line),  (c)  a»i—  oto 
+  2So  (blue-shifted  line). 


By  straightforward  calculation  one  can 
expression  may  be  rewritten  in  the  form 

show  that  this 

St/(F.cu)-^'f't-'-i,,/24?, 

(12) 

where 

A,-{2Ao/R2)'~'m'~-mi‘i+ti\ 

(13) 

«u6*(£?a>o+^f>i  )/(6tf+$?). 

(14) 

and 

l/$oJ-l/5£+l/5?. 

(15) 

On  the  other  hand,  were  the  two  sources  uncorrelated, 
the  correlation  coefficient  hq  would  have  zero  value  and 
we  would  then  have,  according  to  Eqs.  (8)  and  (9), 

lSu</\o>)]..wr,  -  (2A/R2)e  ( 1 6) 

Comparison  of  Eq.  (12)  with  Eq.  (16)  shows  that  al¬ 
though  both  the  spectral  lines  have  Gaussian  profiles, 
they  differ  from  each  other.  Since  according  to  Eq.  (15) 
6i><  6o,  the  spectral  line  from  the  correlated  sources  is 
narrower  than  the  spectral  line  from  the  uncorrelated 
sources.  Further,  we  can  readily  deduce  from  Eq.  (14) 
that 

o>ii£  a>o 
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according  as 
a>i$oo- 

Hence  if  ®i  <  <uo  the  spectral  line  (12)  produced  by  the 
correlated  sources  is  centered  on  a  lower  frequency  than 
the  spectral  line  ( 1 6)  from  two  uncorrelated  sources,  i.e., 
it  is  red  shifted  with  respect  to  it;  and  if  <u,  >  o>o  the 
spectral  line  (12)  is  blue  shifted  with  respect  to  the  spec¬ 
tral  line  (16).  Figure  2  illustrates  these  results  by  simple 
examples. 

The  preceding  considerations  show  clearly  the  possibil¬ 
ity  of  generating,  by  means  of  correlations  between 
source  fluctuations,  either  red  shifts  or  blue  shifts  of  lines 
in  the  spectrum  of  radiation  emitted  by  sources  that  are 
stationary  with  respect  to  an  observer. 

I  am  obliged  to  Mr.  A.  Gamliel  for  carrying  out  the 
computations  relating  to  Fig.  2.  This  research  was  sup¬ 
ported  by  the  U.S.  National  Science  Foundation  and  by 
the  U.S.  Air  Force  Geophysics  Laboratory  under  Air 
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li  h»s  recently  been  demonstrated  both  theoretically  and  experimentally  that  the  spectrum  of  light  emmrd  b>  a  source  can  be 
modified  b>  chanting  the  correlation  properties  of  the  source  We  examine  the  possibility  of  using  this  mechanism  to  modulate 
spectra  in  a  desired  manner  A  simple  configuration  is  considered,  consisting  of  two  small  sources,  emitting  identical  spectra!  lines 
and  »e  show  that  by  suitably  correlating  the  sources  one  can  control  the  width,  the  amplitude  and  the  frequency  of  the  emitted 
light  The  possibility  of  generating  set  era!  lines  from  a  single  one  is  also  considered 


1.  Introduction 

Ji  has  been  shown  theoretically  not  long  ago  J 1 ] 
that  the  spectrum  of  light  emitted  b>  a  source  de¬ 
pends  not  only  on  the  source  spectrum  but  also  on 
the  correlation  properties  of  the  source.  This  pre¬ 
diction  has  been  verified  by  experiments  [2a].  The 
spectral  changes  that  may  arise  from  source  corre¬ 
lations  may  be  such  as  to  produce  frequency  shifts  of 
spectral  lines  [  3-5]  either  towards  lower  frequencies 
(red  shifts)  or  towards  higher  frequencies  (blue 
shifts).  Such  shifts  have  recently  been  demonstrated 
by  experiments  with  optical  [2b]  and  also  with 
acoustical  [6]  sources. 

In  the  present  note  we  show  that  source-correla¬ 
tions  can  give  rise  to  other  interesting  modifications 
of  spectra.  We  consider  only  a  simple  arrangement, 
consisting  of  two  small  sources  which  generate  ra¬ 
diation  of  identical  spectra,  and  we  analyze  the  ef¬ 
fects  of  correlation  between  the  two  sources  on  the 
spectrum  of  the  emitted  radiation.  We  show  that 

*  Research  supported  b>  the  National  Science  Foundation  »nd 

b>  the  U  S  Arm>  Research  Office 
1  Alsou  ihe  Institute  of  Optics.  University  of  Rochester 


spectral  lines  can  be  frequency  shifted,  made  nar¬ 
rower  or  broader  and  that  several  lines  may  be  gen¬ 
erated  from  a  single  line  by  this  mechanism.  These 
results  suggest  that  it  might  be  possible  to  develop  a 
new  technique  for  modulating  spectra  in  a  desired 
manner  by  control  of  source  correlations 


2.  Formulation  of  the  problem 

Consider  light  generated  by  two  small  fluctuating 
sources  located  at  points  F,  and  P2.  We  assume  that 
the  fluctuations  are  statistically  stationary  Let 
w)}  and  {Q(P},  to)}  be  the  ensembles  that  represent 
the  source  fluctuations  at  frequency  w.  Further  let 
{V(P,  u))}  be  the  ensemble  that  represents  the  field 
at  P,  generated  by  the  two  sources  (fig  I ).  We  as- 
aume  that  the  spectra  of  the  two  source-distributions 
•re  identical  and  we  denote  them  by  Sq{o>).  More 
specifically 

5<?(cu)«<0*(/,„w)0(/>„w)> 

•(PIP:. to)  Q(P2 . to)), 
where  the  angular  brackets  denote  ensemble  average 
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Fig  I .  Illustrating  the  configuration  and  notation 


It  has  been  shown  in  ref.  [5]  that  the  spectrum  of 
the  field  at  a  point  P  is  given  by 

SL(P,a))sSe(io) 


:{*y+*r24 


PeUa) 


exp[iA(/?:  -R  )) 
RtR: 


]}■ 

(2.1) 


Here  R,  and  R:  denote  the  distances  between  P,  and 
P  and  P2  and  P  respectively:  3  denotes  the  real  pan 
and  Hy(aj).  known  as  the  degree  of  spectral  coher¬ 
ence.  characterizes  the  correlation  between  the  two 
fluctuating  sources.  Explicitly. 

H0{u)  =  [<Q'(Pl'U)Q‘(P:<io)y)lS0{u>).  (2.2) 

The  degree  of  spectral  coherence  satisfies  the 
inequality 

l/ipfoi)!  <1  (2.3) 

for  all  frequencies. 

For  simplicity  we  will  only  consider  the  spectrum 
of  the  emitted  radiation  at  points  located  on  the  per¬ 
pendicular  bisector,  which  we  will  refer  to  as  the  axis, 
of  the  line  joining  the  two  sources.  In  this  case 
uy)  and  eq.  (2.1 )  reduces  to 

Sl(/>,cu)«(2//?:)S<?(cu)  [1  +  .  (2.4) 

We  note  that  when  3tpQ{o))  is  independent  of  to,  the 
spectrum  of  the  field  at  all  axial  points  will  be  pro¬ 
portional  to  the  source  spectrum  S^fo;).  This  in¬ 
cludes  the  case  where  the  two  sources  are  mutually 
incoherent  at  each  frequency  [/^(co)  ■0).  It  also  in¬ 
cludes  the  case  when  StpQ(<o)m  1;  the  two  sources 
are  then  mutually  fully  coherent  at  each  frequency  . 
These  are,  however,  exceptional  cases.  In  general, 
Stp9{<o)  wil  be  frequency-dependent  and  eq.  (2.4) 
shows  that  the  field  spectrum  5(  (oi)  will  then  be  no 
longer  proportional  to  the  source  spectrum  SQ(to). 


Hence,  in  general,  not  only  the  source  spectrum  but 
also  the  correlation  between  the  two  sources  will  de¬ 
termine  the  spectrum  of  the  emitted  light 
Not  long  ago  it  was  shown  theoretically  [  3-5]  that 
in  some  cases  sources-correlations  can  produce  shifts 
of  spectra]  lines  either  towards  lower  frequencies  ( red 
shifts),  or  towards  higher  frequencies  (blue  shifts). 
This  prediction  has  recently  been  verified  experi¬ 
mentally  with  optical  [2b]  as  well  as  acoustical  [6] 
sources.  In  this  note  we  will  show  that  with  suitable 
choices  of  the  correlation  other  types  of  modulations 
of  the  spectrum  of  the  emitted  radiation  may  be 
produced. 


3.  A  realizability  condition 

It  will  be  convenient  to  set 

SL(a))m(R:/2)  SL(P,  to)  .  (3.1) 

We  will  refer  to  St{tu)  as  the  reduced  field  spectrum 
For  the  Mke  of  simplicity  we  will  consider  source- 
correlations  that  are  characterized  by  a  real  degree  of 
coherence.  Then  eq.  (2.4)  becomes 

St(eo)=S0<(o)  [1  +uQ(w)}  .  (3.2) 

It  follows  at  once  from  this  formula,  that  in  terms  of 
Si  and  Sq. 

/i0(w)*[Sl(ai)/50(ai)]-l  .  (3.3) 

From  the  inequality  (2.3)  and  from  eq  (3  3)  it  fol¬ 
lows  that  only  reduced  field  spectra  5(  (ui)  can  be 
generated  for  which 

St(cuK25e(a>) .  (3.4) 

Si  (to)  and  S0I. o> )  are,  of  course,  necessarily  non¬ 
negative. 

Conversely,  when  the  inequality  (3.4)  is  Mtisfied 
one  finds  at  once  from  eq.  (3.3)  that 

-»I</jc(q>)<  1  .  (3.5) 

Since  the  inequality  (3.5)  is  the  only  constraint  that 
the  degree  of  spectral  coherence  p0(w)  must  satisfy, 
we  see  from  eq.  ( 3.4 )  that  any  reduced  field  spectrum 
Stiw)  which  does  not  exceed  twice  the  magnitude  of 
the  source  spectrum  Sq(uj)  can,  in  principle .  be  gen¬ 
erated  by  this  mechanism. 

We  will  now  consider  a  few  examples 
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4.  Change  of  linewidth 

4.1.  Lorentzian  line  shapes 

Let  us  assume  first  that  the  source  spectrum  con¬ 
sists  of  a  spectral  line  of  lorentzian  profile,  viz., 

Se(oj)  =  .40t[(io  —  4-ri)  (4.1) 

(cu&.  r0.  A(.  are  positive  constants  and  rc« cu0). 
Suppose  that  we  wish  to  produce  a  reduced  field 
spectrum  that  consists  also  of  a  line  of  lorentzian 
profile  centered  on  the  same  frequency  a>c.  but  is  of 
a  different  width  and  of  different  strength,  say 


i 


i  it  »  pm  vrii  urn 

|«ftatt  JC  ’  •  ■) 


S((o»)=.«,/[(cu-cuo):  +  r-;]  (4.2) 

(/",.  A,  are  positive  constants  and  r,«:ci>o).  On 
substituting  from  eqs  (4.1)  and  (4.2)  into  the  in¬ 
equality  (3.4)  we  find  that  we  must  have 

M,//*o)K(w))m,.«2.  (4.3) 

where  l^(cu)]m„  is  the  maximum  value,  in  the  range 
Occucoc.  of  the  function 

/0(<u)=((£j-w0):+rs]/l(tu-(u0):+rt]-  (4.4) 

Straightforward  calculation  shows  that  for  all  (po¬ 
sitive)  frequencies  w. 

(/"o/r,)'</0(cu)<l,  iff",  >r0,  (4.5a) 


IIT  I  PM  s  P71  in 
initi  10"  .  >) 


1  </0(ai)€(r0/ri):,  iff,  <f0-  (4.5b) 

Using  these  inequalities  we  deduce  at  once  from  (4.3) 
that  we  must  have 

A,/A0<2,  when >/"0  ,  (4.6a) 

and 


Fig  2  Lint  nirrowing  tnd  lint  broadening  b'  sourer  correla¬ 
tions  The  source  spectrum  and  lhe  reduced  field  spectrum  art 
lines  of  loremrian  profiles  (eqs  <4  1 1  and  (4  2 1]  with  w .  * 
3.g'7xl0  's-'.r,  »  6xl0  s-'.r,  .  2x10  |-‘..4,,.4l« 
0.22  (a). To  ■  Jxl0'0s-',r,  «  5x10"  «  2  lb) 

The  curves  art  normalized  so  that  the  source  spectrum  has  ihe 
value  units  at  the  center  frequencs  w,. 


Al/Ao*.2(ri/roy,  When  /",<ro.  (4.6b) 

In  the  first  case  ir,>r0)  the  emitted  (reduced) 
spectral  line  is  broader  than  the  spectral  line  of  the 
source;  in  the  second  case  (f,  <f0)  it  is  narrower. 

With  the  conditions  (4.6)  assumed  to  be  satisfied, 
the  degree  of  spectral  coherence  that  will  produce  the 
reduced  field  spectrum  (4.2)  from  the  source  spec¬ 
trum  (4.1 )  is  obtained  at  once  on  substituting  from 
these  equations  into  the  formula  (3.3).  One  then 
finds  that  the  required  degree  of  spectra!  coherence 
is  given  by 


a,  (tu— tu0):4 ri 

*A0  (to-a>o )J+fl 


These  results  are  illustrated  in  figs  2. 


(4.7) 


4.2.  Gaussian  line  shapes 


Next  let  us  consider  a  source  spectrum  that  con¬ 
sists  of  a  line  of  gaussian  profile  viz., 

Sc(a»)«40  exp[ -((a-Uo)2l2dl]  (4  8) 

(a>o.  6q.  Ad  are  positive  constants  and  <5C«  a>0)  We 
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examine  the  possibility  of  producing  a  reduced  field 
spectrum  that  also  consists  of  a  line  of  gaussian  pro¬ 
file  ceniered  at  the  same  frequenc>  <j  bui  of  dif¬ 
ferent  width  and  different  strength,  say 

S,  (v)*.4,  expl-(w-ioJ'/ 26]]  (4.9) 

(d ,,  A,  are  positive  constants  and  <5,  «a>0).  On  sub¬ 
stituting  from  eqs  (4.8)  and  (4  9)  into  the  ine¬ 
quality  (3.4)  we  deduce  at  once  that  the  following 
condition  must  be  satisfied 

(AJAb)  .  (4.10) 

where  (£(a>)]m,v  is  the  maximum  value,  in  the  range 
0<a><oc,  of  the  function 

£(u)»exp[- (fa>-w..)-V2J)  .  (4.1J) 

with 

l/J=(l/<5-;)-(]/<5j  .  (4.12) 

Consider  first  the  case  when  (line  narrow¬ 

ing)  Then  J>0  and  evidently  (g(a>)]m„eg(iut)) 
*e  1.  Hence  the  realizability  condition  (4.10)  becomes 

A>!A<$2.  (4.13) 

On  the  other  hand  when  <5,  ><5„.  J  becomes  negative 
and  g(u)  has  then  no  upper  bound  in  the  range 
0<u><x.  Hence  a  broader  line  of  gaussian  profile 
centered  at  the  same  frequency  cu,  cannot  be  pro¬ 
duced  by  this  mechanism  However,  in  practice  one 
is  unlikely  to  be  interested  m  situations  where  the 
spectra  Sp(o>)  and  5,  (cu)  have  gaussian  forms  for 
all  frequencies  If  one  requires  that  the  reduced  field 
spectrum  has  gaussian  shape  only  over  a  finite  range 
around  c uc„  say 

tuc,-o€tu€cuc+^,  (4  14) 

where  a  and  /?  are  positive  constants,  the  inequality 
(4.10)  needs  only  to  be  satisfied  when  the  maximum 
of  g(io)  is  taken  over  the  restricted  range  (4.14).  In¬ 
stead  of  the  inequality  (4.10)  we  then  have  the 
constraint 

(A,/A(,)np(y:/2\J\)<2,  (4.15) 

where  >■  is  the  largest  of  the  constants  a,  ft. 

Let  us  return  to  the  first  case  (d,  <  d0)  The  degree 
of  spectral  coherence  needed  to  achieve  this  modi¬ 
fication  of  the  spectral  line  is  according  to  eqs  ( 3.3 ), 
(4.8)  and  (4.9)  given  by 


lUMtt  10  »  j 

Fig  3  Line  narrowing  and  line  broadening  b>  source  cortela- 
lions  The  source  spectrum  and  the  reduced  field  specuw  are 
lines  with  gaussian  profiles  |eqs  (4  8)  and  (4  9 1)  o\er  ihc  fre- 
quenc>  range  shown  in  the  figures  with  u,  c  3  S"'x  1C  '  s"  . 
4,  -  6x10"  -  2x10"  s‘'./4,/.4,  «  2  (a),  and  i ,  « 

18  x  10"  s*  ‘.A  JA  m  0.5  (bl  The  normalization  is  the  same  as 
■n  figs  2 


Vo(co)  =  (A ,/Ao)  expl-  (a/-a0):/2J]  - 1 ,  (4  16) 

where  6  is  defined  by  eq.  (4.12).  In  the  second  case 
( 6 1  >  <50)  the  degree  of  spectral  coherence  is  given  by 
eq  (4.16)  only  for  frequencies  that  are  within  the 
range  (4.14);  for  frequencies  outside  this  range  the 
degree  of  spectral  coherence  can  take  on  arbitrary 
values,  subject  to  the  constraint  expressed  b'  eq 
(3.5). 

Some  computed  curves  illustrating  these  results  are 
shown  in  fig  3. 
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S.  Change  of  width  of  a  lorentzian  line  accompanied 
by  shift  of  the  center  frequency 

Next  let  us  consider  the  possibility  of  changing  not 
only  the  width  of  the  line  but  also  its  center  fre¬ 
quency,  from  to  to,  say.  Suppose  that  the  source 
spectrum  Sv(to)  is  again  the  single  spectral  line  (4.1 ) 
of  lorentzian  profile,  but  that  the  reduced  field  spec¬ 
trum.  w  hilst  also  a  line  of  lorentzian  profile,  is  cen¬ 
tered  at  a  different  frequency  to  i.e.  that 

Si(w)=-4,/[(cu-a)l):+rf]  (5.1) 

(to,,  r,.  A,  are  positive  constants,  r,«cti|).  On 
substituting  from  eqs.  (4.1 )  and  (5.1 )  intoeq.  (3.4) 
we  deduce  at  once  that  the  following  inequality  must 
now  be  satisfied: 

(AM  Lf,(w)]m,v^2.  (5.2) 

Here  [/,(cu)]m„.  is  the  maximum  value  in  the  range 
0  <  oj  <x  of  the  function 

/,(ai)«  ((tt-ct>u):+re,]/((w-C£>0):  +  r;].  (5.3) 


Fig  4.  Frequcnc)  shift  b>  source  correlations  The  source  spec¬ 
trum  and  the  field  spectrum  arc  lines  with  lorenmar.  profiles 
with  uL.  ■  3.877x  10*‘  s"  '.r,.  *  5x  10  s'  .T  *1x10  s"  . 
A,!A.  m  1.58Sx  I0'\  The  normalization  is  the  same  as  in  figs 


spectra  given  by  eq  (4.1)  we  wish  to  generate  a  field 
whose  reduced  spectrum  has  the  form 


Unlike  in  the  case  considered  in  sec.  4.1  (when 
cu,  =  cu0)  an  explicit  expression  for  the  maximum 
value  of  this  function  cannot  readily  be  obtained 
However  it  seems  that  with  suitable  choices  of  the 
constants  that  specify  the  reduced  field  spectrum 
(5.1 )  the  inequality  (5.2)  can  be  satisfied  for  all  non¬ 
negative  frequencies  cu.  The  degree  of  spectra!  co¬ 
herence  which  gives  rise  to  the  reduced  field  spec¬ 
trum  (5.1)  is  then  obtained  on  substituting  from  eqs 
(51 )  and  (4.1 )  intoeq  (3.3).  The  resulting  expres¬ 
sion  is 


»v(<o) 


Aj_  (to-co0):  +rj 
au  (u)-oj,y  +r:, 


1 . 


(5.4) 


An  example  of  spectral  changes  in  which  both  the 
line  width  and  the  center  frequency  are  modified  by 
source  correlations  of  this  type  is  presented  in  fig.4 


6.  Generation  of  several  lines  from  a  single  line 

As  a  last  example  we  consider  the  possibility  of 
generating  from  a  source  spectrum  that  consists  of  a 
single  line  of  lorentzian  profile  a  field  spectrum  that 
consists  of  several  lines  of  lorentzian  form  More 
specifically,  with  two  sources  that  have  identical 


A. 


5,  (to)-  - TrTTT 

„Ti(w-tu„)-  +r- 


(6.1) 


where  A',  A„.  tu„.  and  are  positive  constants  and 
r„«cu„  (1 «  n$.V).  For  this  to  be  possible,  the  fol¬ 
lowing  condition  obtained  on  substituting  eqs  (6.1) 
and  (4.1 )  into  eq.  (3.4).  must  be  satisfied 

t  (AJA0)  LUto))  m»\  <2.  (6.2) 


Here  |/,.(£u))n,„  is  the  maximum  value,  in  the  range 
0«  to  <  oc,  of  the  function 

f„(oi)~  l ( u> -u>0):  +/■£]/[ (cj-oO:  +  T;]  (6  3) 

Assuming  that  the  constraint  (6.2)  is  satisfied,  the 
degree  of  spectral  coherence  needed  for  generating 
the  reduced  field  spectrum  (6.1 )  is  obtained  at  once 
on  substituting  from  eqs  (6.1 )  and  (4.1 )  into  eq 
(3.3).  One  then  finds  that 


(fcl-CUei^-f  rl\ 

(w-t on):+r:J 


(6.4) 


An  example  of  the  generation  of  a  field  spectrum 
consisting  of  three  lines  from  a  source  spectrum  con¬ 
sisting  of  a  single  line  is  illustrated  m  fig  5. 
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Fig  5.  Oeneration  of  three  spectral  lines  from  a  single  spectra) 
line,  all  oflorentzian  profiles,  with  «u,.  ■  3.877  x  10” s*  ’.To  * 
5x10"  s  ' 1 .  The  parameters  that  speci  f>  the  three  lines  a  re 
w,  m  3.8769x  10  ‘s-'.r,  «  8x  10'  I"1.  A,IAV  m  3  4x10-*; 
iu;  >  3.8'>70xlO"s-',  r:  m  5  x  10‘  $" A;/A,  m  6  1x10*  >. 
tu,  *  3.8 — I  x  10  ‘  s*1.  T,  ■  8x  10'  x-'.AxlA^  m  3  4x10-’. 
The  normalization  is  the  time  as  in  figs  2 

7.  Discussion 

We  have  considered  in  this  note  a  very  simple  ra¬ 
diating  system  consisting  of  two  small  sources  with 
identical  spectra  and  we  showed  that  by  appropri- 
atel>  correlating  the  two  sources,  the  spectrum  of  the 
emitted  radiation  can  take  on  may  different  forms. 
With  systems  consisting  of  a  larger  number  of  ra¬ 
diating  sources  one  could,  of  course,  produce  more 
diverse  spectral  changes  than  those  considered  in  this 
note. 

The  basic  question  that  we  have  not  addressed  is 
how  to  produce  the  prescribed  correlations.  Several 
methods  for  generating  and  modifying  source  cor¬ 
relations,  at  least  for  secondary  sources,  have  been 
developed  in  recent  years.  They  include  the  use  of 
scattering  by  liquid  crystals  under  the  influence  of  an 
external  d.c.  field  (7),  the  use  of  rotating  ground  glass 


plates  [8]  and  of  holographic  filters  (9).  interaction 
of  light  with  ultrasonic  waves  (10).  imaging  and 
lensless  feedback  systems  (11)  and  use  of  achro¬ 
matic  Fourier  transform  lenses  [2].  Unfortunately 
none  of  these  techniques  appear  to  be  flexible  enough 
to  produce  the  kind  of  correlations  that  are  needed 
to  modulate  spectra  in  the  manner  that  we  discussed 
in  this  note.  However,  the  fact  that  controlled  cor¬ 
relations  have  been  produced,  at  least  to  a  limited 
extent,  seems  to  be  promising.  It  might  be  possible 
to  control  correlations  of  sources  of  electromagnetic 
radiation  more  easily  at  lower  frequencies,  e  g  in  the 
microwave  region.  Correlations  between  acoustical 
sources  could  probably  also  be  controled  more  easily 
(6}.  Such  techniques,  when  utilized  to  generate 
spectral  modulation,  might  perhaps  find  potential 
applications  in  communications  systems 
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ABSTRACT 

A  new  method  for  modulating  spectra  by  correlating  source  fluctuations  is  discussed.  Potential  advantages  and  limitations  of 
this  method  are  considered  for  both  scalar  and  electromagnetic  sources. 

1.  INTRODUCTION 

Current  methods  for  modulating  spectra  may  be  divided  into  two  main  types.  The  first  is  active  modulation  in 
which  the  spectrum  is  modified  by  controlling  the  gain  or  feedback  of  the  system.  The  second  type  is  passive  modulation 
where  the  source  emits  a  very  broad  spectrum  and  the  modulation  is  accomplished  by  filtering  out  or  attenuating  specific 
spectral  intervals. 

In  this  paper  we  discuss  a  new  method  for  spectral  modulation  that  is  based  on  controlling  source  correlations.  We 
consider  the  effects  of  correlating  the  fluctuations  of  two  small  scalar  sources.  We  review  some  of  the  possible  spectral 
modifications  and  examine  the  advantages  and  limitations  of  this  approach.  Next  we  present  an  example  in  which  the  sources 
are  two  partially-correlated  linear  dipoles;  it  demonstrates  correlation  effects  on  the  spectrum  and  on  the  radiation  pattern. 

2.  SCALAR  SOURCES 

Consider  two  small  scalar  sources  located  at  points  Pj  and  P2  (see  Fig.  1).  We  assume  that  the  source  fluctuations 
are  characterized  by  stationary  ensembles  (Q(Pj,a)))  and  (Q(P2,(o))  that  give  rise  to  field  fluctuations  represented  by  an 
ensemble  { V(P,co)} .  We  further  assume  that  both  sources  have  identical  spectra,  i Jt.  that 

Sq(g»  =  (Q-(P1,cd)Q(P1,(o))  *  (Q*(Pro»Q(P2,G»)  ,  (1) 

where  the  asterisk  denotes  complex  conjugation  and  the  angular  brackets  denote  ensemble  averaging. 


Fig.  1.  Notation  relating  to  the  calculation  of  the 
field  spectrum  produced  by  two  scalar  sources. 

In  this  notation  the  spectrum  of  the  field  at  a  point  P  is  expressed  by  a  similar  formula 

Sv(P,<0)  =  (v\P,(D)V(P,£0))  , 

and  the  degree  of  correlation  Pq(w)  is  given  by1 

(Q*(Pr<o)Q(P2,o>)) 


UQ(w) 


Sq(co) 


(2) 


(3) 


It  has  been  shown  that  under  the  conditions  stated,  the  spectrum  of  (he  field  at  any  point  P  is  given  by 


2.3 


SPie  Vol  976  Sievsvcet  Optics  (f$88i  /  137 


Sy(CO)  *  Sq((D) 


kCRj-R.n 

e 

*1*2  . 


(4) 


where  5t  denotes  the  real  part  and  k*(ofc  is  the  wavenumber  ytyitH  with  frequency  to. 

In  order  to  simplify  the  equations  and  illustrate  their  physical  implications  it  is  instructive  to  consider  first  the 
spectrum  at  observation  points  located  on  the  axis  (a  perpendicular  bisector  to  the  line  adjoining  the  two  sources).  At  such 
points  R^R  «R  and  the  formula  (4)  for  the  field  spectrum  reduces  to 


Sv(P,co)  =  ^SQ(G>)[l  +  9^Q(co)]  .  (5) 

R 

Since  -1  S  91Pq(io)  S  1,  it  is  apparent  from  Eq.  (5)  that  changes  of  the  source  correlations,  represented  by  Mq(w).  give  rise 
to  modulation  of  the  field  spectrum  Sy(co),  subject  to  the  constraint 


0  £  Sv(C0)  £  “2Sq(C0).  (6) 

lx 

It  follows  from  this  inequality  that  the  field  spectrum  Sy(ta)  must  vanish  at  every  frequency  for  which  SQ(to)  vanishes.  In 
other  words,  this  mechanism  does  not  produce  new  frequency  components. 

Both  limits  in  Eq.  (6)  occur  when  the  sources  are  fully  correlated.  Other  choices  of  source  correlations  can  produce 
field  spectra  that  are  significantly  different  from  their  corresponding  source  spectra.  To  illustrate  the  effect  of  source 
correlations  let  Sq(oj)  consist  of  a  line  of  Loren tzian  shape,  centered  at  frequency  a>0  and  of  width  8Q,  Le., 


V'») ! 


1 


1  +  (co-cd0)2/8o 


and  with  degree  of  correlation  p^(to)  given  by 


HQ«0)  = 


1  +  (0)  -  ©J  )2/By  1 +  (CD  -  ci>2)2/8* 


-1 


(7) 


(8) 


In  Eqs.  (7)  and  (8)  Aj,  are  positive  constants  chosen  so  that  ! p^fti))  |  Z  1  throughout  the  frequency  interval  of  interest. 

«V  *°r  <02‘  ®2  ***  Positive  constants  of  dimension  tec~  .  For  the  sake  of  simplicity  the  dimensions  are  not 

indicated  in  Figs.  2. 4-6.  The  resulting  field  spectrum  Sy(w)  at  points  on  axis  is  obtained  on  substituting  from  Eqs.  (7)  and 
(8)  into  Eq.  (5). 

Examples  of  field  spectra  computed  from  these  formulas  are  shown  in  Fig.  2.  We  see  that  while  for  fully  correlated 
sources  (I  Hq(io)I=  1)  and  uncorrelaied  sources  (IPq(<i))|*  0),  the  spectrum  of  the  field  is  proportional  to  the  spectrum  of  each 

source,  appreciable  modifications  of  the  emitted  spectrum  are  achieved  with  appropriate  choice  of  the  degree  of  correlation. 

In  general,  given  a  source  spectrum  Sq(co)  and  a  desired  field  spectrum  Sy(o>)  the  appropriate  choice  of  the  degree  of 
correlation  is  obtained  from  Eq.  (5)  as 


R2  V<*» 

u  (&})  =  — - - - 

v  ]  2  S_(co) 


-1 


(9) 
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Of  course  only  those  field  spectn  Sy(<o)  can  be  generated  which  are  consistent  with  the  constraint  ||1q(s1.  A  limitation  of 

this  method  is  that  there  is  an  upper  bound,  for  each  frequency  tt,  for  the  strength  of  the  emitted  field  spectrum.  This  upper 
bound  is  given  by  the  spectrum  that  is  produced  when  the  two  sources  are  fully  correlated  [cf.  Eq.  (6)]. 


Frequency  (arbitrary  units) 

Fig.  2.  Field  spectra  on  axis  for  different  types  of  source  correlations,  calculated 
from  Eqs  (5),  (7)  and  (8).  The  constants  were  chosen  as  follows:  o>0*  4, 6Q*  0.25, 

Aj=  1.95,  (Dj*  3.5,  0.03,  A2*  2,  co2*  4.5,  and  S2«  0.03. 

Returning  to  Eq  (4)  we  note  that  the  field  spectrum  also  depends  on  the  direction  of  observation  because  of  the 
presence  of  Rj  and  R^  We  postpone  the  discussion  of  the  angular  dependence  and  of  complex-valued  degree  of  correlation  to 

the  electromagnetic  case  that  we  will  now  consider. 

3.  LINEAR  DIPOLES 

Consider  two  linear  dipoles  located  at  distances  y^  from  the  origin  and  oscillating  in  the  z  direction  as  shown  in  Fig. 
3.  Let  the  electric  polarization  vectors  of  the  dipoles  be  given  by 

PjCr.coJsp^tDjSCr-y^i,  P2(r,co)  =  p2(to)5(r  +  y^)4 .  (10) 


Here  y  and  z  are  unit  vectors,  6  is  the  Dirac  delta  function  and  p^co),  (i=l,2)  represent  the  polarization  fluctuations  as  a 
function  of  frequency.  If  we  assume  that  these  fluctuations  have  the  same  spectrum,  i.e.  that 


(lPj(<0)|2)  =  (lp2(G>)|2)  *  SQ  (0))  , 


(11) 


it  can  be  shown  that  the  spectrum  of  the  emitted  field  in  the  far  zone  in  the  direction  specified  by  a  unit  vector  u  is4 

S(Ru,co)  =~-k4  SQ(o))  sin28  Jl  +9^MQ(G>)exp(2iky0sine  sin$)jj  .  (12) 


Here  C  is  a  constant  depending  on  the  choice  of  units  and  (0,$)  art  the  spherical  polar  coordinates  of  the  unit  vector  u,  with 
the  polar  axis  taken  along  the  z  direction. 
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(13) 


To  simplify  formula  (12)  we  express  the  degree  of  correlation  in  the  form 

|Xq(co) ■  I j1q(g>) I  e**a)  , 

where  y((o)  is  a  real  valued  function*  of  to.  It  follows  that  when  the  two  dipoles  are  fully  correlated  I  Mq(<o)  I  *1  the 
spectrum  of  the  field  in  the  far-zone  is  given  by 

[SCRu.to)]^  «  2  k4  SQ(to)  sin2®  cos2  [kyQsine  sin<>  +  v(m)]  .  (14) 


Similarly  when  the  dipoles  are  uncorrelated  Pq(o))  *  0  and  the  far-field  spectrum,  obtained  from  Eq.  (12)  is 


[S(Ru,oj)] 


uncorr. 


-^•k4So(t0)sin20  . 
IT  W 


(15) 


Z 


Fig.  3.  Notation  relating  to  the  calculation  of  the  far-field 
spectrum  produced  by  two  linear  dipoles  Pj  and  Pj. 

For  every  degree  of  correlation  Pq(co)  *  0  the  frequency-dependence  of  the  far  field  spectrum  is  coupled  to  the  angular  variables 
through  the  product  ky0sin9sin$.  To  illustrate  this  remark  consider  observation  points  in  the  x,  y-plane.  The  spectrum  is 
then  a  function  of  frequency  to  and  of  the  azimuthal  angle  $  only,  viz., 

S(R,  6=|- ,  $;co)  =  -^-k4  Sq(£d)  jl  +  |nQ(co)|  cos [2ky0sin$  +  2y(0))]j  .  (16) 

R 

Let  us  consider  some  implications  of  Eq.  (16).  Suppose  the  spectrum  Sq(w)  of  each  dipole  consists  of  a  single 

Lorentzian  line  [Eq.(7)J.  We  now  compare  the  angular  distribution  of  the  far  field  spectrum  when  the  degree  of  correlation  is 
given  by  Eq.  (8)  and  when  the  dipoles  are  fully  correlated.  Both  cases  are  examined  for  fixed  frequencies  and  for  a  fixed  value 
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of  yQ.  Figures  4  and  S  show  the  angular  distributions  of  the  spectra  far  two  selected  frequencies.  It  is  to  be  noted  that  larger 
values  of  yQ  give  rise  to  multiple  lobes  in  the  radiation  pattern  (see  Fig.  6),  but  it  can  be  shown  that  such  systems  will 
radiate  relatively  little  power. 


•o* 


Fig.  4.  Angular  distribution  of  the  far  field  spectra 
S(R,0=n/2,$;  co)  of  fully  correlated  sources  (dashed 
line)  and  for  Pq  given  by  Eq.  (8)  (solid  line).  Both 
curves  are  for  a  fixed  frequency  to  *  4  and 


so1 


Fig.  5.  Angular  distribution  of  the  far  field  spectra 
S(R,0«*/2,$;  to)  of  fully  correlated  sources  (dashed  line) 
and  for  Pq  given  by  Eq.  (8)  (solid  line).  Both  curves  are 
for  a  fixed  frequency  to  «  3.515  and  k^*!. 
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Fig.  6  Angular  distribution  of  the  far  field  spectra 
S(R,0=r/2,$;  to)  of  fully  correlated  sources  (dashed  line) 
and  and  for  Pq  given  by  Eq.  (8)  (solid  line).  Both  curves 
are  for  a  fixed  frequency  to*  4  and  k^y0«3. 


Another  consequence  of  Eq.  (16)  is  that  when  the  observation  points  are  located  along  the  x-axis  (sin$=0)  the  far 
zone  spectrum  of  the  electromagnetic  field  exhibits  the  same  frequency  dependence  as  we  found  before  for  the  scalar  case  [Eq. 
(5)].  For  other  directions  of  observation  the  far  zone  spectrum  has,  in  general,  a  more  complicated  form  due  to  the  presence  of 
the  fKtor  cos[2kyQsin^42y((0)].  With  the  same  choices  of  the  source  spectrum  and  of  the  degree  of  correlation  as  before,  we 

show  in  Figs.  7  and  8  the  far  field  spectrum  along  the  x-and  the  y  axes. 
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On  comparing  Fig  7  with  the  conupooding  figure  2  for  the  malar  case,  we  note  that  the  curves  in  Fig.  7  show 
asymmetry  not  present  in  Fig.  2.  This  difference  is  due  to  the  presence  of  the  factor  k4  in  Eq.  (16)  relating  to  the 
electromagnetic  case. 


t  t 


t 

TJw 


/Q<“)  - 1 


\  n0(«ifllv»nbyEq.(8) 

\i 

JjJb . T— 

S  4  S  •  7 


Frequency  to  (eibirary  unite) 

Fig.  7  Far  zone  field  spectrum  for  points  located  on 
the  x-axis. 


Frequency  c»  (erbttary  units) 

Fig.  8.  Far  zone  field  spectrum  for  points  located  on 
they-axis. 


DISCUSSION 

We  have  shown  in  this  note  how  spectral  modulation  can  be  achieved  by  controlling  correlations  between  source 
fluctuations.  Both  scalar  sources  and  electromagnetic  sources  were  considered  and  they  illustrate  typical  spectral  features  that 
may  be  produced  by  relatively  simple  choices  of  the  functional  forms  of  the  correlation  coefficient.  The  main  limitation  of 
this  method  is  the  requirement  that  all  spectral  components  of  the  desired  spectrum  must  be  present  in  the  source  spectrum. 

The  practical  difficulty  of  inducing  the  desired  correlation  in  scalar  or  electromagnetic  sources  has  been  overcome  in 
some  special  cases.6-*  At  present  it  seems  that  more  complicated  functional  forms  of  the  degree  of  correlation  are  easier  to 
achieve  with  microwaves  than  in  the  optical  region  of  the  electromagnetic  spectrum. 
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We  examine  the  spatial  directivity  and  the  temporal  spectrum  of  the  field  radiated  by  two  partially  correlated  linear 
dipoles.  The  two  dipoles  are  parallel  and  have  identical  spectra  in  their  source  excitation.  Expressions  are  derived 
for  the  radiated  spectrum,  for  the  total  emitted  power,  and  for  the  directivity.  W'e  illustrate  these  results  by 
calculating  the  radiated  spectrum  for  a  particular  choice  of  the  frequency-dependent  correlation  function.  Results 
for  the  limiting  cases  of  fully  correlated  dipoles  and  for  uncorrelated  dipoles  are  also  derived. 


INTRODUCTION 

Recent  advances  in  coherence  theory  include  investigations 
of  the  effects  of  source  correlations  on  the  spectrum  of  the 
radiated  field.1-7  The  type  of  source  correlation  that  gives 
rise  to  line  shifts  in  the  observed  spectrum  was  considered  in 
a  number  of  those  studies.  Some  of  the  theoretical  predic¬ 
tions  have  since  been  verified  experimentally  by  several 
groups.5-7  One  of  the  possible  applications  of  the  present 
theory  is  a  method  for  modulating  field  spectra  that  is  based 
on  the  controlling  of  source  correlations.  The  effect  of  spec¬ 
tral  modulation  was  illustrated  in  an  earlier  publication  in 
which  a  simple  configuration  of  two  scalar  sources  was  ana¬ 
lyzed  for  observation  points  that  are  located  on  the  axis.4 

In  the  present  paper  we  are  interested  in  spectral  effects 
that  are  also  associated  with  the  direction  of  observation.  In 
our  treatment  we  consider  a  system  of  two  partially  correlat¬ 
ed  linear  dipoles.  First  we  derive  expressions  for  the  ob¬ 
served  spectrum  at  an  arbitrary  point  in  the  far  zone;  then 
we  examine  how  the  degree  of  correlation  affects  the  spec¬ 
trum  of  the  field  observed  in  particular  directions  and  its 
effects  on  the  angular  distribution  of  the  radiant  intensity 
for  fixed  frequencies.  In  Sections  3  and  4  expressions  are 
derived  for  the  total  radiated  power  and  the  directivity  of 
the  system  as  a  function  of  the  degree  of  correlation  and  the 
spatial  separation  between  the  two  dipoles. 

Throughout  the  analysis,  we  compare  our  results  for  the 
partially  correlated  dipoles  with  the  results  in  the  well- 
known  limiting  cases  of  fully  correlated  and  uncorrelated 
dipoles.  This  comparison  gives  a  valuable  measure  for  the 
range  of  possible  spatial  and  spectra]  modulation  effects 
that  can  be  produced  by  controlling  source  correlations. 

1.  FAR-ZONE  SPECTRUM  OF  PARTIALLY 
CORRELATED  DIPOLES 

Consider  two  linear  dipoles  situated  at  points  ±yo  and  vi¬ 
brating  in  the  z  direction  as  shown  in  Fig.  1.  Let 

P,(r,  1)  «  p,(t)i(r  -  y$)i,  (1.1a) 


the  dipoles  as  a  function  of  time,  and  S'  and  i  are  unit  vectors 
in  the  positive  y  and  z  directions,  respectively.  W  e  assume 
that  Pj(t)  are  random  functions  of  time,  characterized  by 
stationary  ensembles. 

In  the  space-frequency  representation,8  the  electric  Hertz 
potential  of  the  field  produced  by  the  dipoles  is  given  by 


n,(r,w) 


i  />,( w) 


npUkRj) 

*i 


+  P2(tt.') 


exp(ikR:)' 

*2  V 


(1.2) 


in  which  R i  *  Ir  -  yoS'l.  r2  *  |r  +  yo>l  and 

Pj(w)  *  j"  Pj(t)e‘*'dt  (J  «=  1,  2).  (1.3) 


In  order  to  calculate  the  radiant  intensity  in  the  far  zone,  it 
is  sufficient  to  evaluate  the  magnetic  field  only;  this  is  given 
by  the  expression 


B(r,  u>)  -  -ikV  x  U,(r,  u>).  (1.4) 


Since  Il,(r,  u>)  is  a  vector  along  the  z  direction,  it  follows  that 
vxn((r,«)*^i-^jn,(r,«).  (1.5) 


On  substituting  Eq.  (1.2)  into  Eq.  (1.5),  we  obtain  the  formu¬ 
la 

,xn.<r 

For  field  points  r  *  rti  in  the  far  zone  kR,  »  1  (j  *  1,  2),  we 
have 


R\~r  —  y<>G  •>,  R2  -  r  +  y^Ci  •  $, 


(1.7) 


in  which  Ci  is  a  unit  vector  for  the  polar  radius  r. 

Using  Eq.  (1.6)  and  expressions  (1.7),  we  may  express  the 
far  zone  B  field  in  the  form 


P2(r,  t)  -  p2(t)«(r  +  y^)i  (1.1b) 

be  the  electric  polarization  vectors  of  the  two  dipoles.  Here 
Pj(t)  0  “  1.  2)  characterize  the  polarization  fluctuations  of 


B(rCt, «) - k2  sin  0[p,(o>)  expi-iky^  •  5) 

+  p2(w)exp(ik>0Ci  •  $)!$  as  kr -*  ®.  (1.8) 
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Fig.  1.  Illustration  of  the  configuration  and  notation  used  in  this 
study. 


Here,  we  have  used  the  spherical  polar  coordinates  (r,  6,  <t>) 
with  the  polar  axis  along  the  z  direction  and  with  4>  *  sin  9 
(cos  <t>y  -  sin  4>x). 

The  radiant  intensity  J(&,  v),  i.e.,  the  power  per  unit 
frequency  at  frequency  u>,  per  unit  solid  angle  centered  about 
direction  specified  by  the  unit  vector  u,  is  given  by 

J(u,  u>)=  lim  r^u  •  (S{‘\ru,w)),  (1.9) 

kr 

where  u)  is  the  Poynting  vector  in  the  far  zone  and 

the  angle  brackets  denote  the  ensemble  average.  In  terms  of 
the  B  field,  we  have 

J(l H,  u)  «  lim  r2  •  [B(ru,  w)  X  ti\*  X  B(ru,  w)}, 

*r—  Sir 

(1.10) 

where  %  denotes  the  real  part.  Using  the  identity 

6  •  (B  X  £)•  X  B  ■  IB  X  Ct\2,  (1.11) 

we  simplify  Eq.  (1.10)  for  the  radiant  intensity  to  the  form 
«/(£,«)«  Urn  r2^-(iB(rii,«)|2).  (1.12) 

*r—«  Or 

Let  us  now  assume  that  both  dipoles  have  the  same  spec¬ 
trum,  i.e.,  that 

<lp,(w)!2)  -  <|/>2(u0l2  *  Sp(u>).  (1.13) 


We  also  introduce  the  complex  degree  of  spatial  coherence  at 
frequency  w  that  characterizes  the  correlation  between  the 
two  dipoles  by  the  formula9 


SPM 


(1.14) 


On  substituting  from  expression  (1.8)  and  Eqs.  (1.13)  and 
(1.14)  into  Eq.  (1.12),  we  obtain  the  following  expression  for 


the  radiant  intensity  produced  by  the  two  partially  correlat¬ 
ed  dipoles: 

J(a,u)-^SJu)  sin2  8 

4*  v 

X  |1  +  [Mp(w)exp(2iky0  sin  9  sin  *)]).  (1.15) 

If  we  express  the  degree  of  correlation  in  the  form 

Mp(a))  ■  Up(«)l  exp[i2iM <*>)],  (1-16) 

in  which  the  phase  2iMw)  corresponds  to  the  steering  angle  in 
the  fully  coherent  case,10  we  find  that  the  radiant  intensity  of 
the  two  partially  correlated  linear  dipoles  is  given  by 

J(d,  w)  -  ^  Sp(u)  sin2  911  +  lup(u>)l 

X  coe[2Ay0sin  9  sin  4>  +  2^(w)]|.  (1.17) 


2.  EFFECTS  OF  SOURCE  CORRELATION  ON 
THE  RADIATED  SPECTRUM 

We  now  examine  some  special  cases  that  help  to  illustrate 
Eq.  (1.17).  First  let  us  consider  two  uncorrelated  dipoles. 
In  this  case  mp(<*>)  =  0,  and  Eq.  (1.17)  reduces  to 

Sp(u>)  sin2  9.  (2.1) 

As  one  might  expect,  the  same  expression  is  obtained  if  the 
radiating  spectrum  originated  from  a  single  dipole  located  at 
the  origin,11  whose  strength  equals  the  sum  of  the  strength  of 
the  two  uncorrelated  dipoles.  Similarly,  when  the  dipoles 
are  fully  correlated,  i.e.,  when  UP(«)I  *  1,  Eq.  (1.17)  for  the 
radiant  intensity  reduces  to 

[«7(fi,  <*>)]„,„  -  —  Sp(«)  sin2  9  cos2[fey0  sin  9  sin  <t>  +  ^(w)]. 

(2.2) 

We  observe  that  for  every  frequency  w,  the  phase  angle  Hu) 
and  the  product  ky0  completely  determine  the  angular  dis¬ 
tribution  of  the  radiated  power. 

Returning  to  the  general  case  in  which  the  two  dipoles  are 
partially  correlated,  we  note  that  when  the  point  of  observa¬ 
tion  is  on  the  x  axis  (9  “  *72,  <t>  *  0),  the  radiant  intensity 
according  to  Eq.  (1.17)  is  given  by 

J{St,  w)  -  Sp(«)|l  +  *Im>)1|.  (2.3) 

The  result  for  this  special  case  is  in  the  same  form  as  the 
corresponding  expression  for  the  radiant  intensity  from  two 
small  partially  correlated  scalar  sources  [cf.  Ref.  4,  Eq.  (2.4)]. 
The  only  significant  difference  is  in  the  factor  k*  appearing 
in  Eq.  (2.3),  in  which  the  product  k*Sp(w)  is  shifted  to  higher 
frequencies  relative  to  Sp{w). 

The  radiation  pattern  produced  by  the  partially  correlat¬ 
ed  dipoles  differs  from  those  produced  in  the  two  limiting 
cases  of  fully  correlated  and  completely  uncorrelated  dipoles 
in  several  ways.  If  we  denote  the  direction  for  which  the 
radiant  intensity  is  a  maximum  by  fim,  it  can  readily  be 
shown  that  the  maximum  possible  radiant  intensity  [</(um, 
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w)]mu  from  the  two  dipoles  is  obtained  when  they  are  fully 
correlated,  and  it  is  given  by 


ck* 

£sP(»). 


(2.4) 


On  the  other  hand,  when  the  dipoles  are  partially  correlated 
the  maximum  radiant  intensity  is  smaller  by  a  factor  Vail  + 

|jip(a>)|]. 

In  the  two  limiting  cases  of  fully  correlated  dipoles  and 
uncorrelated  dipoles,  the  nulls  of  the  radiant  intensity  dis¬ 
tribution  are  determined  by  the  factor  sin2  8  in  Eq.  (2.1)  and 
by  the  factor  sin2  6  cos2[feyo  sin  8  sin  <t>  +  ^(oj)]  in  Eq.  (2.2). 
By  contrast,  it  follows  from  Eq.  (1.17/  that  when  the  dipoles 
are  partially  correlated  (0  <  Up(<u)l  <  1)  there  are  no  nulls  of 
the  radiant  intensity  outside  the  plane  9*0.  This  fact  is 
significant  in  connection  with  the  theorems  regarding  the 
approximation  of  desired  radiation  patterns  by  arrays  of 
such  sources.12 

For  a  numerical  illustration,  let  the  dipole  spectrum  be  a 
Loren tzian  line  of  width  So  centered  at  frequency  wq,  i.e., 


_ 1 _ 

1  +  (a;  —  u0)2/S02 


(2.5) 


We  choose  a  real -valued  degree  of  correlation  given  by 


Mp<u)  ■  A,  exp 


(u>  —  u,,)2 
25 12 


(<*>  -  o)2)2"] 

“srj"1* 

(2.6) 


Here  u\,  u>2,  5 1,  and  S2  are  positive  constants,  and  we  select 
real-valued  At  and  A2  so  that  mp(u)  satisfies  the  constraint 

Up(w)l  £  1  (2.7) 

throughout  the  frequency  range  of  interest.  We  note  that,  if 
we  choose  a  real-valued  degree  of  correlation,  the  steering 
angle  is  set  to  zero  [i.e.,  2^(«)  s  0]  in  the  following  examples. 

Figures  2  and  3  show  the  angular  distributions  of  the 
radiant  intensity  in  the  (x,  y)  plane  (8  *  *72)  for  two  fully 
correlated  dipoles  and  for  partially  correlated  dipoles  whose 
degree  of  correlation  is  given  by  Eq.  (2.6).  In  Fig.  2  the 
angular  distribution  is  calculated  for  the  center  frequency  u 
~  w0  of  the  polarization  spectrum;  Fig.  3  shows  the  angular 
distribution  at  a  different  frequency.  When  the  separation 
between  the  dipoles  increases,  the  lobe  structure  becomes 
more  complicated.  In  Figs.  4  and  5  we  show  the  angular 
distributions  of  the  radiant  intensity  for  dipole  separation 
*oyo  *  3.  Comparing  the  angular  distributions  for  dipole 
separation  *oy0  *  1  (Figs.  2  and  3)  and  k^yo  ■  3  (Figs.  4  and 
5),  we  note  that  at  each  frequency  the  number  of  lobes  of  the 
two  choices  of  correlations  is  identical,  although  their  angu¬ 
lar  distributions  are  somewhat  different. 

At  a  fixed  direction  in  space,  we  can  compare  the  spectrum 
produced  by  the  partially  correlated  dipoles  with  the  spec¬ 
trum  produced  by  fully  correlated  dipoles.  In  Fig.  6  we  show 
the  spectra  for  two  choices  of  correlation  when  the  observa¬ 
tion  is  along  the  x-direction.  Although  it  is  not  shown  in  the 
figure,  we  remark  that  the  spectrum  produced  by  the  fully 
correlated  dipoles  differs  from  that  of  the  uncorrelated  di¬ 
poles  by  a  factor  of  2,  whereas  the  spectrum  produced  by  the 
partially  correlated  dipoles  displays  features  that  are  not 


present  in  the  polarization  spectra.  We  also  point  out  that 
the  spectrum  of  the  field  differs  from  the  Lorentzian  line 
shape  of  the  polarization  spectrum  because  of  the  factor  of  k 4 
in  Eq.  (1.17).  A  similar  comparison  of  spectra  for  observa¬ 
tion  points  along  the  y  direction  is  shown  in  Fig.  7.  Here  we 
note  that  the  spectrum  produced  by  fully  correlated  dipoles 
has  a  lower  magntiude  than  the  spectrum  produced  by  the 


90° 


Fig.  2.  Angular  distributions  of  the  radiant  intensity  at  an  offset 
frequency  w  «  4  sec*1  and  a  dipole  separation  koy0  •  1  when  the  two 
dipoles  are  partially  correlated  (solid  curves)  and  fully  correlated 
(dashed  curves).  The  constants  used  in  the  calculations  are  (in 
reciprocal  seconds)  wo  *  4,  *  3.5,  w2  m  4.5, 4o  “  0.25,  and  6i  «  5:  » 

0.03.  These  simple  numerical  constants  can  be  scaled  to  any  fre¬ 
quency  band  that  is  being  considered. 


90° 


Fig.  3.  Angular  distributions  of  the  radiant  intensity  at  an  offset 
frequency  u  m  3.52  sec*1  and  a  dipole  separation  *o>o  “  1  when  the 
two  dipoles  are  partially  correlated  (solid  curves)  and  fully  correlat¬ 
ed  (dashed  curves).  The  constants  used  in  the  calculation  are  the 
same  as  in  Fig. 2. 
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Pig.  4.  Angular  distribution  of  the  radiant  intensity  at  an  offset 
frequency  u  «  4  sec'1  and  a  dipole  separation  kayo  “  3. 


90° 


Fig.  5.  Angular  distribution  of  the  radiant  intensity  at  an  offset 
frequency  u>  ■  3.52  sec'1  and  a  dipole  separation  kayo  ■  3. 


partially  correlated  dipoles.  This  fact  may  be  explained  by 
Fig.  2,  for  example,  in  which  the  radiant  intensity  in  the  y 
direction  is  higher  for  the  partially  correlated  dipoles.  Fig¬ 
ures  6  and  7  demonstrate  that  the  spectral  features  of  the 
partially  correlated  dipoles  vary  considerably,  depending  on 
the  direction  of  observation,  whereas  the  spectrum  produced 
by  the  fully  correlated  dipoles  remains  essentially  un¬ 
changed. 

The  two  limiting  cases  of  fully  correlated  and  uncorrelated 
dipoles  constitute  the  boundaries  of  possible  modification  of 
spectra.  In  Fig.  8  we  illustrate  the  range  of  modulation  that 
can  be  achieved  at  every  frequency  w,  by  the  variation  of  the 
magnitude  of  the  degree  of  correlation  Up(o>)l  and  the  phase 
2^(u>).  The  figure  shows  four  concentric  circles,  (i)-(iv), 


representing  the  relative  ranges  of  the  angular  distributions 
of  the  radiant  intensity  in  the  (x,  y)  plane.  The  external 
circle  (i)  corresponds  to  the  limit  of  fully  correlated  dipoles. 
As  one  may  observe  from  Eq.  (2.2),  the  radiant  intensity  of 
two  correlated  dipoles  may  attain  any  value  inside  the  circle 
of  radius  ck*Sp(u)/2r.  The  location  of  the  peak  is  controlled 
by  the  choice  of  the  steering  angle  2^(u>).  Similarly,  it  fol¬ 
lows  from  Eq.  (2.1)  that,  when  the  dipoles  are  uncorrelated, 
the  radiant  intensity  is  constrained  to  values  on  a  circle  (iii) 
of  radius  ck*Sp(w)/4r.  When  the  dipoles  are  partially  corre¬ 
lated,  we  see  from  Eq.  (1.17)  that  the  radiant  intensity  can 
have  any  value  in  an  annular  domain  bounded  by  circles  (ii) 
and  (iv).  It  is  also  apparent  from  this  representation  that 


Fig.  6.  Comparison  of  spectra  along  the  x  direction  when  the  two 
dipoles  are  partially  correlated  (solid  curves)  and  fully  correlated 
(dashed  curves) .  The  constants  used  in  the  calculation  are  the  same 
as  in  Fig.2. 


Frequency  u  [sec-1) 

Fig.  7.  Comparison  of  spectra  along  the  y  direction  when  the  two 
dipoles  are  partially  correlated  (solid  curves)  and  fully  correlated 
(dashed  curves). 
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Fig.  8.  Range  of  control  over  angular  distribution  of  the  radiant 
intensity.  The  shaded  area  indicates  the  region  in  which  the  maxi¬ 
ma  and  minima  of  the  radiant  intensity  for  partially  correlated 
dipoles  are  found. 
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Fig.  9.  The  behavior  of  the  function  /  in  Eq.  (3.3)  as  a  function  of 
the  product  2Ayo- 


unless  the  two  dipoles  are  fully  correlated  there  are  no  nulls 
of  the  radiant  intensity  outside  the  plane  6  m  0. 

3.  TOTAL  EMITTED  POWER 

The  total  power  Pu*(v)  radiated  by  the  system  at  frequency 
u>  is  defined  by 

P^M  -  [  J(Ci,  a>)dn,  (3.1) 

in  which  the  integration  extends  over  the  whole  4r  solid 
angle.  On  substituting  Eq.  (1.17)  into  Eq.  (3.1)  and  per¬ 
forming  the  integration,  we  find  that 

PmM  "  V*4Vw)|l  +  f(2*yomMpM]l.  (3.2) 

in  which 


and  j0(z)  and  jt(z)  are  spherical  Bessel  functions.  The  maxi¬ 
mum  of  the  function  f(z)  can  be  shown  to  occur  when  z  ■  0, 
and  it  has  the  value  f(0)  ■  1  (see  Fig.  9).  It  follows  that  the 
upper  bound  for  the  total  radiated  power  is 

[PM].*,  -  V*4SpM-  (3.4) 

When  the  two  dipoles  are  uncorrelated,  the  total  radiated 
power  is  given  by 

[P(«)]u»eorr  «  V*4ty«>-  (3.5) 

On  comparing  the  total  powers  o"  two  partially  correlated 
dipoles  PM  and  uncorrelated  dipoles  (PMluncorr.  we  see 
from  Eqs.  (3.2)  and  (3.5)  that 

rp/y~-  -  1  +  /<2*y0)tf  [„>)].  (3.6) 

{ty&flixncoTT 

Since  the  function  f(z)  given  by  Eq.  (3.3)  decreases  rapidly 
with  increasing  z,  it  is  evident  that  the  total  radiated  power 
from  two  dipoles  that  are  separated  by  a  distance  much 
larger  than  a  wavelength  is  equivalent  to  the  total  radiated 
power  from  two  uncorrelated  dipoles. 


4.  DIRECTIVITY 

The  directivity  D(u,  w)  of  a  radiating  system  is  defined  by 
the  ratio13 


D(u,  u) 


4t J(u,  ii>) 
Pu>.M  ’ 


(4.1) 


where  Ct  is  the  direction  of  observation.  Here  we  are  particu¬ 
larly  interested  in  the  maximum  directivity  DM  -  D(um,  u) 
that  occurs  in  a  particular  direction  um.  It  follows  from  Eqs. 
(1.17)  and  (3.2)  that  the  directivity  of  the  two  partially  corre¬ 
lated  dipoles  is  given  by 


0  10  20 
2ky0 


Fig.  10.  Comparison  of  the  directivity  D(u>)  for  fully  correlated 
dipoles  (dashed  curve)  and  uncorrelated  dipoles  (solid  curve)  as  a 
function  of  the  dipole  separation  2Jkyo-  In  both  cases  the  value  of 
the  steering  angle  was  2^<w)  ■  r. 


f(z )  ■  %[/0(*)  -  )\(z)/z\ 


(3.3) 
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^3  1  +  lp,(w)l 

"  2  |1  +  /(2*yo>*|M,<«Mr 


(4.2) 


It  is  straightforward  to  show  that,  when  the  dipoles  are 
uncorrelated,  (D(u>)]uiworT  ■  3/2.  Similarly,  when  the  dipoles 
are  fully  correlated  the  directivity  is  given  by 

DM  "  |1  +  /(2fcy0)  cos(2^(u))]j  ’  (4,3) 


tion  between  the  dipoles  on  the  total  radiated  power  and  the 
directivity  of  the  two  dipole  system.  In  both  cases  we 
showed  that  the  degree  of  the  correlation  was  coupled  with  a 
function  of  the  separation  constant  feoy o.  With  increasing 
separation  between  the  two  dipoles,  the  total  radiated  power 
becomes  independent  of  the  correlation,  and  as  a  result  the 
directivity  becomes  bounded  by  a  constant  value. 
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in  which  2 is  the  phase  of  pp(a>). 

Equations  (4.2)  and  (4.3)  also  indicate  that,  when  the  two 
dipoles  are  separated  by  a  distance  that  is  much  larger  than  a 
wavelength,  the  directivity  of  partially  correlated  dipoles 
approaches  the  value  D{u)  —  3[1  +  Up(o>)l]/2,  whereas  the 
directivity  of  fully  correlated  dipoles  approaches  the  value 
[DMjnm  ~  3  (see  Fig.  10).  This  result  is  consistent  with  our 
earlier  remarks  regarding  the  total  power  radiated  from  such 
dipoles.  In  addition,  if  we  use  Eq.  (1.16)  in  Eq.  (4.2)  and 
differentiate  with  respect  to  |pp(«)|,  it  is  straightforward  to 
show  that  the  directivity  is  maximized  for  a  fixed  value  of 
/(2A>0)cos[2^(<»)]  when  l^p(ai)l  «  1,  i.e.  when  the  dipoles  are 
fully  correlated. 

5.  DISCUSSION 

Although  the  properties  of  the  Fields  produced  by  linear 
dipoles  have  been  known  for  many  years,  most  of  the  previ¬ 
ous  studies  were  confined  to  the  limits  of  fully  correlated  and 
uncorrelated  dipoles.  In  this  paper  we  described  the  spec¬ 
trum  of  the  field  and  the  angular  distribution  of  the  radiant 
intensity  produced  by  two  linear  dipoles  with  polarization 
fluctuations  that  are  subject  to  an  arbitrary  choice  of  fre¬ 
quency-dependent  correlations. 

In  Section  2  we  used  a  simple  choice  of  frequency-depen¬ 
dent  correlation  function  to  illustrate  its  effects  on  the  angu¬ 
lar  distribution  of  the  radiant  intensity  for  selected  values  of 
the  frequency.  We  showed  that  the  angular  distribution  for 
different  frequencies  changes  significantly  despite  the  fact 
that  we  used  a  fixed  steering  angle  2^(u>)  a  0  in  all  cases.  We 
also  showed  that  the  spectrum  observed  in  one  direction  may 
have  distinct  features  that  are  not  present  in  other  direc¬ 
tions. 

In  all  the  examples  that  we  considered,  the  precise  spec¬ 
tral  and  angular  details  depend  on  the  choice  of  the  correla¬ 
tion  function.  We  believe  that  the  essential  facets  of  the 
theory  are  illustrated  well  by  our  simple  choice  of  correlation 
function  (Eq.  (2.6)].  Additional  choices  of  the  correlation 
function  are  discussed  in  Ref.  4. 

In  Sections  3  and  4  we  examined  the  effects  of  the  correla¬ 
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Scattering  of  polychromatic  light  by  a  medium  whose  dielectric  susceptibility  is  a  random  function  of  position  is 
considered  within  the  accuracy  of  the  first  Born  approximation.  It  is  shown,  in  particular,  that  if  the  two-point 
spatial  correlation  function  of  the  dielectric  susceptibility  has  Gaussian  form  and  the  spectrum  of  the  incident  light 
has  a  Gaussian  profile,  the  spectrum  of  the  scattered  light  may  be  shifted  toward  the  shorter  or  the  longer 
wavelengths,  depending  on  the  angle  of  scattering  The  results  are  analogous  to  those  derived  recently  in  connec¬ 
tion  with  radiation  from  partially  coherent  sources  [Nature  (London)  326, 363  (1987)]. 


I.  INTRODUCTION 

It  was  predicted  theoretically  not  long  ago,  that,  in  general, 
correlations  between  source  fluctuations  produce  changes  in 
the  spectrum  of  the  emitted  light.12  For  radiation  from 
planar  secondary  sources  this  prediction  was  subsequently 
verified  by  experiment.3  It  was  also  shown  theoretically 
that  the  changes  may  be  such  as  to  produce  red  shifts  or  blue 
shifts  of  spectral  lines, 4-6  and  this  prediction  too  has  been 
verified.7-9  A  similar  effect  can  also  be  expected  to  arise 
with  acoustical  waves  and  was,  in  fact,  observed  not  long 
ago.10  Some  other  modifications  of  spectra  due  to  source 
correlations  have  also  been  discussed.11 

Because  of  the  well-known  analogy  that  exists  between 
the  processes  of  radiation  and  scattering,  one  might  expect 
that  similar  phenomena  will  arise  when  a  polychromatic 
wave  is  scattered  by  a  medium  whose  dielectric  susceptibil¬ 
ity  is  a  random  function  of  position.  We  show,  in  the 
present  paper,  that  this  indeed  is  the  case.  First  we  derive 
an  expression,  valid  within  the  accuracy  of  the  first  Born 
approximation,  for  the  spectrum  of  the  scattered  light  in 
terms  of  the  spectrum  of  the  incident  light  and  the  two-point 
spatial  correlation  function  of  the  dielectric  susceptibility  of 
the  medium.  We  then  show  that  if  the  spectrum  of  the 
incident  light  consists  of  a  single  line  of  Gaussian  profile  and 
the  correlation  function  of  the  dielectric  susceptibility  is  also 
a  Gaussian  function,  the  spectrum  of  the  scattered  field  will 
consist  of  a  line  that  has  approximately  a  Gaussian  profile. 
However,  this  line  is,  in  general,  red  shifted  or  blue  shifted 
with  respect  to  that  of  the  incident  light,  depending  on  the 
angle  of  scattering.  This  result  may  appear,  at  first  sight,  to 
contradict  the  well-known  fact  that  there  is  no  frequency 
change  in  linear  scattering  on  a  time-invariant  medium  and 


that  different  frequency  components  of  the  incident  and  also 
of  the  scattered  light  are  uncorrelated.  The  resolution  of 
this  apparent  paradox  is  briefly  discussed  in  the  concluding 
section. 

2.  EXPRESSION  FOR  THE  SPECTRUM  OF 
THE  SCATTERED  FIELD 

Let  us  consider  a  field  incident  upon  a  scattering  medium 
occupying  a  finite  volume  V.  Suppose  that  the  incident 
field  propagates  in  a  direction  specified  by  a  real  unit  vector 
lo-  We  do  not  assume,  however,  that  the  field  is  monochro¬ 
matic;  rather  we  consider  it  to  fluctuate  at  each  point,  gener¬ 
ally  in  a  random  manner,  characterized  by  an  ensemble  that 
is  statistically  stationary.  The  cross-spectral  density 
W^fri,  r2,  <*>)  of  the  incident  field  at  points  whose  location  is 
specified  by  position  vectors  rt  and  rz  may  be  expressed  in 
the  form12 

W<"(r„  r2,  w)  -  (U<,**(r„ «)£/"( r* «)>.  (2.1) 

where  |  u>)|  represents  a  statistical  ensemble  of  random 

functions,  all  of  the  form 

(/‘‘(r,  w)  •  a(«)exp(ifcj0  •  r),  (2.2) 

with 

(2.3) 

c 

c  being  the  speed  of  light  in  vacuo.  In  Eq.  (2.2)  a(«)  are 
(generally  complex)  frequency-dependent  random  vari¬ 
ables,  and  the  angle  brackets  in  Eq.  (2.1)  denote  the  expecta¬ 
tion  value,  taken  over  the  ensemble  of  the  incident  field. 
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The  spectrum 

^‘’(r.w)*  ^‘’(r.r.u)  -  <l/,',(r,«)I/,’(r.w)>  (2.4) 

of  the  incident  field  is  then,  according  to  Eq-  (2.2),  given  by 

#■'(«)  *  ^''(r,  w)  -  (o*(w)a(«)>  (2.5) 

and  is  seen  to  be  independent  of  position. 

We  assume  that  the  scatterer  is  weak  in  the  sense  that  the 
amplitude  of  the  scattered  field  (/,((r,  u)  is  small  compared 
with  the  amplitude  of  the  incident  field  lll/ul| «  ll/'1!).  We 
may  then  calculate  U{,)  on  the  basis  of  the  first  Born  approx¬ 
imation,  which  yields13 

u>)  ■  o(u>)  j  F(r', u>)G(r,  r',  w)exp(iAs0  •  r')dV, 


where 


F(r,  u>)  *  (w/c)2i?(r,  uj), 


i|(r,  <o)  *  [n2(r,  w)  -  l]/4*  being  the  dielectric  susceptibility 
of  the  scattering  medium,  n(r,  u >)  its  refractive  index,  and 


C(r,r',  u,)--CX-PP(r-~r')] 


the  outgoing  free-space  Green’s  function. 

We  will  consider  only  the  scattered  field  in  the  far  zone.  If 
we  set  r  ■  rs,  (I2  ■  1),  the  Green's  function  may  then  be 
approximated  by  its  asymptotic  form  (see  Fig.  1) 

G(r,  r',  oj)  ~ exp(-iks  ■  r')  (*r— «).  (2.9) 

On  substituting  from  the  formula  (2.9)  into  Eq.  (2.6),  we  find 
at  once  that 


Fig.  1.  Illustrating  the  notation  relating  to  the  asymptotic  approxi¬ 
mation  (2.9). 

On  substituting  from  Eq.  (2.12a)  intoEq.  (2.13)  and  on  using 
Eq.  (2.5),  we  find  at  once  that 

S<“'(rs,  w) 

-  ~  ~  *0>.  ~  l0),  »)S"'(u:).  (2.14) 

Up  to  now  we  have  assumed  that  the  scattering  medium 
was  deterministic.  Let  us  now  assume,  instead,  that  its 
physical  properties  are  characterized  by  random  functions 
of  position  but  are  independent  of  time.  An  example  of  such 
a  medium  would  be  the  atmosphere  under  conditions  when 
its  temporal  fluctuations  are  slow  enough  to  be  ignored 
Then,  for  each  u,  the  dielectric  susceptibility  r;(r,  u)  will  be  a 
random  function  of  r,  and  consequently  r)(K,  u)  will  be  a 
random  function  of  K.  A  meaningful  measure  of  the  spec¬ 
trum  of  the  scattered  field  in  the  far  zone  is  then  the  expecta¬ 
tion  value  of  the  right-hand  side  of  Eq.  (2.14),  taken  over  an 
ensemble  of  different  realizations  of  the  scattering  medium; 
i.e.,  it  is  given  by  the  expression 

S"“'(rl,  w) 


(/*’(«,  ui)  *  a(w)  ~ j  F(r',  ui)exp[-ifc(s  -  |0)  •  r’]dV,  *<>)•  i0),  a)]),S"  (w).  (2.15) 


where  we  have  written  U{m)  rather  than  (/*’  to  stress  that  the 
expression  on  the  right-hand  side  of  Eq.  (2.10)  represents 
the  scattered  field  in  the  far  zone.  If  we  introduce  the 
Fourier  transform  of  the  scattering  potential  by  the  formula 

£(K,  w)  «  |  F(r\  «)exp(-iK  •  r')dV,  (2.11) 
Eq.  (2.10)  becomes 

l/*’(ri,  «)  -  a(w)F IMS  -  «0),  w]  ~ .  (2.12) 

or,  using  Eq.  (2.7), 

l/*’(rl,  u)  -  a(w)^Z  n|fe(l  -  l0),  w]  (2.12a) 

where  v  (K,  w)  is  the  Fourier  transform,  defined  by  a  formula 
of  the  form  (2.11),  of  the  dielectric  susceptibility  i?(r',  u>). 

Let  us  next  consider  the  spectrum  of  the  scattered  field  in 
the  far  zone.  It  can  be  determined  from  the  following  for¬ 
mula  analogous  to  Eq.  (2.4),  which  is 

S<”<rl, «)  «  <l/-‘*(ri,  w)l/*’(ri,  »)>.  (2.13) 


Here  the  angle  brackets,  with  the  subscript  77,  denote  the 
average  value  taken  over  this  ensemble.  It  is  to  be  noted 
that  the  expression  (2.15)  involves  two  averaging  proce¬ 
dures,  one  over  the  ensemble  of  the  incident  field  (cf.  Eq 
(2.5))  and  the  other  over  an  ensemble  of  scatterers.  We 
implicitly  assume  here  that  these  two  sources  of  randomness 
are  mutually  independent. 

Although  the  formula  (2.15)  for  the  far-field  spectrum 
involves  an  average,  obtained  from  experiments  with  macro- 
scopically  similar  but  microscopically  different  scatters,  one 
can  in  some  cases  deduce  the  value  of  this  average,  at  least  to 
a  good  approximation,  from  experiments  performed  with  a 
single  scatterer.  For  example,  often  the  necessarily  finite 
size  of  the  detector  aperture  will  provide  spatial  averaging, 
which  is  essentially  equivalent  to  ensemble  averaging  14 15 
The  value  of  the  ensemble  average  may  also  be  obtained 
sometimes  by  the  use  of  a  moving  aperture  in  front  of  the 
scatterer.16 

If  we  express  5j  in  terms  of  its  Fourier  inverse  n  and  inter¬ 
change  the  order  of  averaging  and  integrations,  we  find  that 


(n*(K,,i*>)n(K2, «)>„•  j  j  C,(r,\ r2\ w) 

X  exp[-i(K2  •  r2'  -  K!  •  r,'))d3r1'd3r2'.  (2.16) 


(2.13) 


1144  J.  Opt.  Soc.  Am.  A/Vol.  6,  No.  8/August  1989 


Wolf  et  at. 


where 


C,(r,',  r2\  uj)  -  <i?*(r1\  w)v(r2’,  <*>)),  (2.17) 

is  the  spatial  correlation  function  of  the  dielectric  suscepti¬ 
bility.  If  we  introduce  its  six-dimensional  Fourier  trans¬ 
form 

£,(KJfK2,  w)  *  J  j  C,( r,',  r2',  u) 

X  exp[-i(K,  .  r,'  +  K2  •  r2')]d3r,'dV.  <2-18) 
Eq.  (2.16)  becomes 

<q*(K„  w)n(K2.  «)>,  -  f,(-K„  K2.  u).  (2.19) 

Finally,  on  substituting  from  Eq.  (2.19)  into  Eq.  (2.15),  we 
obtain  the  following  expression  for  the  spectrum  of  the  scat¬ 
tered  field  in  the  far  zone: 

S'-'(rs,  u)  «=  ~  («/c)4C„|-fc<l  -  s0),  -  SoK 

r 

(2.20) 

The  formula  (2.20)  shows  that  the  spectrum  of  the  scat¬ 
tered  field  in  the  far  zone  differs,  in  general,  from  the  spec¬ 
trum  S^Hu.)  of  the  incident  field  by  the  effect  of  two  multi¬ 
plicative  factors,  namely,  a  factor  proportional  to  u4  and  the 
factor  C„.  The  a-4  factor  is  a  reflection  of  the  fact  that  on  the 
microscopic  level  the  medium  responds  to  the  incident  field 
as  a  set  of  dipole  oscillators  [cf.  Ref.  17,  Secs.  2.2.1  and  2.2.3]. 
The  other  factor,  £„  may  be  regarded  as  representing  the 
correlation  that  exists  between  them.  Because  C„  in  Eq. 
(2.20)  depends  on  the  momentum  transfer  vector  K  *  fed  - 
S0),  it  is  a  function  of  the  angle  of  scattering,  and  hence,  in 
general,  the  far-field  spectrum  will  be  different  in  different 
directions  of  observation. 

The  situation  is  somewhat  different  in  the  idealized  spe¬ 
cial  case  when  the  dielectric  susceptibility  of  the  scattering 
medium  is  completely  uncorrelated,  i.e.,  when 

C„< r,\  r2\  u)  *  /(ai)<5l3,(r2'  -  r,')  when  r,'  e  V,  r2'  e  V 

■  0  otherwise, 

(2.21) 

where  4,3)(r')  is  the  three-dimensional  Dirac  delta  function 
and  Hu)  is  a  nonnegative  function  of  frequency.  We  then 
have,  according  to  Eqs.  (2.18)  and  (2.21), 

C„[-fed  -  S0),  fed  -  S0),  u)  -  Hu) V,  (2.22) 

and  the  formula  (2.20)  reduces  to 

lS<-»(rS.  -  7  (7  Jnu)SMiu).  (2.23) 

We  see  that  in  this  special  case  the  spectrum  of  the  scattered 
field  is  the  same  for  all  angles  of  scattering  and  that  it  differs 
from  the  source  spectrum  Hu)  only  by  a  multiplicative  factor 
that  is  proportional  to  the  dipole  term  u4. 

It  seems  worthwhile  to  stress  that  in  deriving  Eq.  (2.20)  no 
assumption  was  made  regarding  homogeneity  or  isotropy  of 
the  scatterer.  When  the  scatterer  is  statistically  homoge¬ 
neous,  as  is  often  the  case,  the  formula  (2.20)  considerably 
simplifies,  as  we  will  now  show. 


3.  SPECTRUM  OF  THE  FIELD  SCATTERED  BY 
A  STATISTICALLY  HOMOGENEOUS  MEDIUM 

When  the  medium  is  statistically  homogeneous,  the  correla¬ 
tion  function  C,(ri',  r2',  u)  will  depend  on  t\  and  r2'  only 
through  the  difference  r2'  -  r j',  and  we  will  then  write 

C,(r,',  r2',  w)  -  C,(r2'  -  r,'.  u)  when  r,'  e  V,  r2'  e  V 

« 0  otherwise. 

(3.1) 


In  this  case  Eq.  (2.18)  becomes 
C,(K„  K2,  u) 

*  j  j  C,(^ r— •/>  «)exp[-i(K,  •  r,'  +  K2  •  r2')]d3r,'d3r2',r 

(3.2) 


and  hence 

C,[-fe(J-i0).  fed  - 10).  <*>] 


V  Jv 


C,(r2'  -  r,',  «)exp[-»fed  -  *  <r2'  ~  r,')]d3r,'d3r2'. 

(3.3) 


If  we  change  the  variables  of  integration  from  ri',  r2'  to  r, 
r'  by  setting 


r  -  <*,'  +  r2')/2,  r'  -  (r2'  -  r,')  (3.4) 

and  assume  that  the  linear  dimensions  of  the  scattering 
volume  are  large  compared  both  with  the  correlation  length 
of  the  dielectric  susceptibility  [the  effective  r'  range  of  IC„(r', 
<!.')[]  and  with  the  wavelength  X  *  2-rc/u,  Eq.  (3.3)  gives 

c,l-fed  -  «0),  fed  - 10).  w]  *  VC,[fe(s  -  l0), «].  (3.5) 

where  C„(K,  u)  is  the  three-dimensional  Fourier  transform 
of  C,(r\  a;),  viz., 

C,(K,  w)  ■  |  C,(r\  «)exp(-iK  •  r')dV.  (3.6) 

On  substituting  from  formula  (3.5)  into  the  general  formula 
(2.20)  we  obtain  the  following  simpler  formula  for  the  spec¬ 
trum  of  the  scattered  field: 

S<’>(r£,  u)  *  Ifed  -  *„>,  u>]S('V).  (3-7)  [ 

This  formula  again  shows  that  the  spectrum  of  the  incident 
field  is,  in  general,  modified  by  interaction  with  the  scatter¬ 
ing  medium.  The  modification  arises  from  a  dipole  contri¬ 
bution  that  is  proportional  to  the  fourth  power  of  the  fre¬ 
quency  and  by  the  influence  of  the  dielectric  susceptibility- 
correlations,  characterized  by  the  correlation  function  C„(r', 
u). 


4.  FREQUENCY  SHIFTS  GENERATED  BY 
SCATTERING 

To  illustrate  the  effect  of  the  dielectric  susceptibility  corre¬ 
lations  of  the  scattering  medium,  let  us  consider  the  situa¬ 
tion  when  the  correlation  function  is  a  three-dimensional 
Gaussian  distribution,  i.e.,  when 
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c,(r'’ a) "  «p(“r'2/2ff2)- 

where  A  and  o  are  positive  constants,  i.e.,  they  are  indepen¬ 
dent  of  both  r'  and  “. 

The  three-dimensional  Fourier  transform  (defined  by  Eq. 
(3.6)]  of  the  expression  (4.1)  is 

C,(K,  “)  «  A  exp(— (K»)2/2].  (4.2) 

Now  with  K  «  k($  -  |0)  we  have,  since  f  and  fo  are  unit 
vectors  and  k  «  u/c, 

K 2  *  4Wc)2sin2(0/2),  (4.3) 

where  0  is  the  scattering  angle  (I  •  So  m  cos  0).  Hence  it 
follows  that,  in  this  case, 

£,[*(*  -  i0),  w]  -  A  expj-2^  a2  sin2^j  ’  (4-4> 

and  the  expression  (3.7)  for  the  spectrum  of  the  scattered 
field  becomes 

S'->(rs, “)  «  ^  (yj  *Xp[-2(c)2  “X,2(1)]S<‘I(W)‘ 


Suppose  that  the  spectrum  of  the  incident  field  is  a  Gauss¬ 
ian  function  of  rms  width  To,  centered  at  frequency  a>o,  i.e., 
that 


S("(u>)  m  B  exp 


-(“  -  fc)0)2 


where  B,  u>0  and  T0  are  positive  constants.  On  substituting 
from  Eq.  (4.6)  into  Eq.  (4.5),  we  see  that  the  right-hand  side 
contains  the  product  of  two  Gaussian  functions.  With  the 
help  of  a  product  theorem  for  Gaussian  functions  estab¬ 
lished  in  Appendix  A,  we  show  in  Appendix  B  that  the 
resulting  expression  for  the  spectrum  of  the  scattered  field 
can  be  written  in  the  form 

S(-*(r,  0;  «)  «  N(r)H(0)«*  expj- 1  ft*-*]*}  ’  <47> 

where  we  have  now  written  S("l(r,  0;  w)  rather  than  S,m>(rs, 
<•>).  The  various  quantities  that  appear  on  the  right-hand 
side  of  Eq.  (4.7)  are  defined  by  the  following  formulas: 


VAB 
‘  cV  ’ 

(4.8a) 

(4.8b) 

“0 

'o2<0)’ 

(4.8c) 

To 

'  o(0)  ’ 

(4.8d) 

H’GXSHSD" 

(4.8e) 

Xo— •  M.80 

We  see  from  Eq.  (4.8e)  that  a(6)  >  1  and  that  the  equality 
holds  only  (i.e.,  that  a(0)  ■  1]  when  either  <r  «  0  (completely 
uncorrelated  scatterer)  or  To  •  0  (monochromatic  light)  or 
when  8  ■  0  (forward  scattering).  Hence  except  in  these 
special  cases 

6(0)  <  «0,  f  (0)  <  To,  (4.9) 

implying  that  the  Gaussian  function  in  Eq.  (4.7)  is  centered 
on  a  frequency  that  is  lower  than  the  central  frequency  uo  of 
the  incident  light  and  its  rms  width  is  smaller.  The  fact  that 
6(0)  <  uo  implies  that  the  Gaussian  function  in  the  expres¬ 
sion  (4.7)  is  centered  at  a  lower  frequency  (i.e.,  is  red  shifted) 
with  respect  to  the  Gaussian  spectral  line  of  the  incident 
light,  the  magnitude  of  the  shift  depending  on  the  angle  of 
scattering, 18  0.  On  the  other  hand,  the  factor  o>4  in  the 
expression  (4.7)  is  an  increasing  function  of  the  frequency 
and  hence  will  produce  a  shift  toward  the  higher  frequencies 
(i.e.,  a  blue  shift).  Consequently  the  spectrum  of  the  scat¬ 
tered  field  will  be  either  red  shifted  or  blue  shifted  with 
respect  to  the  Bpectrum  of  the  incident  light,  depending  on 
the  magnitudes  of  these  two  contributions.  It  is  shown  in 
Appendix  C  that  the  maximum  of  the  spectrum  (Eq.  (4.7)]  of 
the  scattered  field  at  scattering  angle  0  occurs  at  frequency  v 
■  u>o'(0),  where 

"”'w,-i5s{1+[1+“’(,l(f)T}'  (4101 

It  is  customary,  especially  in  astronomy,  to  specify  frequen¬ 
cy  shifts  by  the  quantity 

V  ~  \>  “o  ~  “o'  , ,  , , , 

2  - — - - - - - — .  (4.11) 

*0  "0 

where  X©  •  2tc/u>o  is  the  original  wavelength  and  V  *  2 rd 
u»o'  is  the  corresponding  shifted  wavelength.  Evidently  r  >  0 
when  the  line  is  red  shifted  and  z  <  0  when  it  is  blue  shifted. 
On  substituting  from  Eq-  (4.10)  into  Eq.  (4.11)  we  readily 
find  that  in  the  present  case 


i+  i  +  oWi-M 


Were  the  scatterer  completely  uncorrelated  (o  ■  0),  we 
would  have,  according  to  Eq.  (4.8c),  a  ■  1  for  all  values  of  0, 
and  Eq.  (4.12)  would  reduce  to 

r  /ar_\ni« 


T 

/4r0yi> 


In  Figs.  2-4  our  main  results  are  illustrated  by  a  number  of 
computed  curves.  Figure  2  shows  spectra  of  the  scattered 
field  at  different  angles  of  scattering  for  some  selected  values 
of  the  parameters.  Figure  3  illustrates  the  behavior  of  the 
height  factor  H(0)  defined  by  Eq.  (4.8b).  The  behavior  of 
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3.3  3.4  3.5  3.6 

OJ  [1015  sec-1] 

Fig.  2.  Spectrum  [in  unit*  of  N(r)j  of  light  scattered  at  various 
angles  8,  from  a  Gaussian  correlated  medium  with  rros  width  a  ■ 
3*o.  *o  m  5500A  (u>o  m  3.427  X  1015  sec-1).  The  spectrum  of  the 
incident  light  is  a  line  with  Gaussian  profile  of  rms  width  r0  ■  10"2 
wo).  The  curve  labeled  8  «  0  also  represents  the  spectrum  (in  units 
of  N(r)]  for  the  case  when  the  scatterer  is  completely  uncorrelated  ( a 
■  0)  (see  Eq.  (4.4)1. 


In  the  caw  when 

(4iy«0)2 «  1,  «  0(1)  (4.14) 


(conditions  that  were  satisfied  in  connection  with  the  com¬ 
puted  curves),  z(8)  as  given  by  Eq.  (4.12)  can  be  approximat¬ 
ed  by  an  expression  that  clearly  indicates  the  roles  of  the 
physical  parameters.  One  finds,  after  a  straightforward  cal¬ 
culation,  that 


z(8)  *  4 


feHfl-e)) 


(4.15) 


where  Lq  ■  c/r0  is  the  coherence  length  of  the  incident  light 
(cf.  Note  18).  Further,  when  the  conditions  (4.14)  are  ful¬ 
filled,  the  fractional  shift  of  the  center  frequency  of  the 
Gaussian  factor  in  Eq.  (4.7),  defined  by  the  expression 


zc(8)  B 


u0  —  u>(6) 


(4.16) 


is  readily  found  to  be  approximately  equal  to 


H  (©) 


Fig.  3.  Height  factor  [H(«)[  [Eq.  (4.8b)),  with  iy<vo  “  lO'2.  The 
values  of  the  correlation  parameter  a/\o  are  indicated  on  the  curves. 


the  relative  frequency  shift  as  a  function  of  the  angle  of 
scattering  is  shown  in  Fig.  4.  We  see  that  for  a  certain  range 
of  directions  around  the  forward  direction  8*0,  the  line  is 
blue  shifted  (r  <  0),  but  the  magnitude  of  the  shift  decreases 
with  increasing  8,  and  eventually  red -shifted  lines  (r  >  0)  are 
produced.  The  lower  curves  in  Fig.  4  show  the  difference 
between  the  values  pertaining  to  scattering  from  a  spatially 
random  medium  of  finite  (nonzero)  correlation  length  (o  > 
0)  and  from  one  that  is  completely  uncorrelated  (<r  •  0). 
This  difference  is  seen  to  be  positive  for  all  angles  of  scatter¬ 
ing  8  ^  0,  indicating  that  when  the  scatterer  has  a  finite 
(nonzero)  correlation  length  the  spectrum  is  always  red 
shifted  with  respect  to  the  spectrum  that  is  produced  with  a 
completely  uncorrelated  scatterer. 


Z  (©) 


m 


Fig.  4.  Relative  frequency  shift  *(8),  (a),  and  the  difference  ft#)  “ 
«(9)  -  iviimvT.  (b).  of  the  spectral  line  of  light  scattered  at  various 
angle*  t  from  a  Gaussian -correlated  medium  with  rms  width  a.  The 
spectrum  of  the  incident  light  is  a  line  with  a  Gaussian  profile  of  rms 
width  f0  -  10'2  wo- 
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Hence  the  first  term  in  the  formula  (4.15)  is  a  red  shift  due 
to  the  spatial  correlations  in  the  physical  properties  of  the 
medium.  The  second  term  is  a  blue  shift  due  to  the  dipole 
factor  u4. 


5.  CONCLUDING  REMARKS 

We  have  considered  in  this  paper  the  effect  of  a  random 
medium  whose  physical  properties  do  not  change  in  time  on 
the  spectrum  of  the  light  scattered  by  it.  It  is  commonly 
stated  that  under  these  circumstances  the  frequency  of  the 
incident  light  will  not  change  by  the  process  of  scattering. 
How  is  it  possible  then  that  our  analysis  indicates  that  such 
scattering  may  produce  frequency  shifts  of  the  spectral 
lines?  The  resolution  of  this  apparent  paradox  lies  in  the 
fact  that  most  previous  work  was  concerned  with  the  scatter¬ 
ing  of  essentially  monochromatic  light,  of  well-defined  fre¬ 
quency,  t>)o,  say.  The  spectrum  of  the  scattered  light  will 
then  necessarily  also  be  essentially  monochromatic  and  of 
the  same  frequency  u>o.  However,  when  the  incident  light  is 
polychromatic ,  the  different  frequency  components  of  the 
incident  light,  which  are  scattered  in  any  particular  direc¬ 
tion,  will  be  scattered  with  different  strengths.  As  a  result, 
the  spectrum  of  the  scattered  light  will  differ  from  that  of 
the  incident  light,  even  though  the  different  frequency  com¬ 
ponents  are  uncorrelated.  As  we  showed,  the  difference 
may  manifest  itself  as  a  frequency  shift.  This  possibility  of 
generating  frequency  shifts  by  scattering  is  analogous  to  that 
discovered  not  long  ago  in  connection  with  effects  of  source 
correlations  on  the  spectrum  of  the  emitted  light.4"9  The 
only  difference  is  that  we  now  deal  with  secondary  sources, 
namely,  with  the  polarization  induced  in  the  scattering  me¬ 
dium  by  the  incident  wave.  The  induced  polarization  will, 
in  general,  be  correlated  over  finite  distances  of  the  scatter¬ 
ing  medium  and  will  thus  imitate  correlations  in  primary 
sources.  It  was  this  analogy  between  scattering  and  radia¬ 
tion  that,  in  fact,  led  to  the  present  analysis. 

Finally  we  might  mention  that  our  results  suggest  a  new 
method  for  determining  correlation  properties  of  refractive- 
index  distribution  in  spatially  random  media,  based  on  the 
analysis  of  the  spectral  changes  produced  by  a  spatially 
random  medium  when  illuminated  by  polychromatic  light  of 
known  spectral  composition. 


APPENDIX  A:  PRODUCT  THEOREM  FOR 
GAUSSIAN  FUNCTIONS 

In  this  appendix  we  will  establish  the  following  theorem. 
Theorem:  If 

G(w  -  w,\  r;)  -  exp[-(u)  -  u>,)72r/l  0  -1.2),  (Al) 

then 

g(<»3  — 1*>,;  r,)G(u>  —  <i)2!  iy 

-  G(w,  -  <4i2:  (Pi2  +  r22)'/2]G(w  -  w;  f ),  (A2) 

where 


7r«-r;  +  -rr  (A4) 

f2  r,2  r22 

To  establish  this  theorem  we  multiply  together  the  two 
formulas  (Al)  and  express  the  product  in  the  form 

G(o)  —  r,)G(o>  —  <o2;  r2)  -  exp|— g(o!)],  (A5) 


£(“>)  *  — -5- ■ -2  fr22(uj  -  w,)2  +  r,2(w  -  u>2)2] 

2*1  r2 


2r,2r,2 


(a2ta;2  —  2bu>  +  c) 


a2  -  r,2  +  r22, 
b  -  aj,r22  +  w2r,2, 
c  —  a!,2r22  +  w22r,2. 

On  completing  the  square  in  Eq.  (A6),  we  find  that 


1  o2  1  |  1 

f2  *  r,2r22  “  r,2  r22’ 


6  o.,r22  +  o>2r,2 

o2'  17  + r22 


By  using  Eqs.  (A7a)-(A7c),  one  can  readily  show  that 


iyiy 

r,2+r2 


(ill,  —  li>2)2. 


On  substituting  from  Eq.  (All)  into  Eq.  (A8)  and  from  Eq. 
(A8)  into  Eq.  (A5),  we  find  that 

G(w  —  r,)G(o>  —  «2;  r2) 

J"  (m,  ~  w2)2  1  T  (m  -  £)2~] 

“  exp[~  2(iy  +  r22)Jexp[~  2f 2  J 

-  G[ui,  -  u>2;  (r,2  +  r22)1/2]G(u,  -  u>;  f ),  (A12) 

where  the  last  step  follows  from  Eq.  (Al).  This  completes 
the  proof  of  the  theorem. 

Some  remarks  about  the  implications  of  this  theorem  are 
in  order.  Equation  (A2)  shows  that  the  product  of  two 
Gaussian  functions  is  proportional  to  a  new  Gaussian  func¬ 
tion,  G( u  -  u>;  f ),  the  constant  of  proportionality  being  G[ui, 
—  o)2;  (T,2  +  r22),/2].  The  center  frequency  and  width  of  the 
new  Gaussian  function  are  u>  and  f,  given  by  the  formulas 
(A10)  and  (A9),  respectively. 

The  relationship  between  u>  and  the  two  original  center 
frequencies  ui,  and  u>2  can  be  expressed  in  a  different  form. 
It  follows  from  Eq.  (A3)  that 

u  *  /,<• 1,  +  fjU2- 


(A3) 


(A13) 
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where 


U 


r,2 


JV  +  i? 


h' 


r,2 


P)  +r2 


(A14) 


Since  f\  +  /2  “  1,  £>  is  a  weighted  average  of  the  two  original 
center  frequencies.  In  particular,  if  0  <  o>i  <  i*>2,  and  neither 
T]  nor  f2  is  zero. 


Ol!  <  to)  <  to)2. 


(A15) 


The  relationship  between  f  and  the  original  widths,  T  i 
and  T  2,  can  be  deduced  from  Eq.(A4).  If  neither  T  i  nor  T2  is 
zero,  we  find  that  the  new  Gaussian  is  narrower  than  either 
of  the  two  original  Gaussians,  i.e.,  that 


f  <  f„  f  <  r2. 


(A16) 


Finally,  it  follows  directly  from  Eq.  (A12)  that  the  maxi¬ 
mum  value  (as  a  function  of  w)  of  the  product  G{u  —  wf, 
fi)G(to)  -  w2;  r2)  occurs  at  the  frequency  w  m  Si  and  is  equal  to 
G[(u>i  -  o>2;  (rt2  +  r22)1/2J,  whereas  the  maximum  values  of 
each  of  the  original  Gaussians  occur  at  w  *  w,  and  are  each 
equal  to  unity.  Therefore,  as  long  as  oil  *  <*>2,  the  maximum 
value  of  the  product  of  two  Gaussians  is  smaller  than  the 
maximum  value  of  either  of  the  two  original  Gaussians. 

APPENDIX  B:  DERIVATION  OF  EXPRESSION 
(4.7)  FOR  SPECTRUM  OF  THE  SCATTERED 
HELD 


We  have,  according  to  Eqs.  (4.5)  and  (4.6), 


(w  —  w0)2 

~r7~ 


where19 


If  we  set 


i  4  2  .  Je\ 

— r  =  —  a  Sin  —  • 

r2  c2  \2> 


G(w  -  v/,  Tj)  *  ex; 

w  VAB 

N  r- 


(Bl) 

(B2) 

(B3) 

(B4) 


cV 

the  expression  (Bl)  becomes 

S'-'(rs, «)  «  Nw*G{v,  r,)G(«  -  w0;  r0).  (B5) 

Now  according  to  the  product  theorem  for  Gaussian  func¬ 
tions  established  in  Appendix  A, 

G(w;  r,)G(u>  -  w0;  r0)  «  G(cu0;  O',2  +  r22),/2]G(o.-  -  £>;  f ), 

(B6) 

where 


and 


w0r,2 

r02  +  r,2 


f2‘r02  r,2' 


(B7) 


(B8) 


It  will  be  convenient  to  rewrite  Eqs.  (B7)  and  (B8)  as 
Si  1 


«o  i  +  dVr,)2 
f  i 


r0  (i  +  (r0/r1)2)1/2 

Now  from  Eq.  (B2)  it  follows  that 
Po  2 


Po  2  _  .  f 6\ 


(B9i 

(BIO) 

(Bill 


or,  when  the  relation  Xot*>o  *  c,  ( Xo  -  hJ2  z )  is  used ,  Eq .  ( B 1 1 ) 
gives 


h*)gk> 


(B12) 


On  substituting  from  Eq.  (B12)  into  Eqs.  (B9)  and  (BlO)  we 
-  find  that 


to)  ^  1 

Wo  a2 

r0  a 


where 


HOSH*’ 


(B13) 

(B14) 

(Bl  5) 


Further  we  have,  if  we  use  Eqs.  (BlO), 

r02  +  r,2  =  r,2 


r.2 


or,  if  we  also  use  Eq.  (B14), 

r02  +  r,2 «  o2r,2. 


[-(31 

m 


(B16) 


Hence  the  first  factor  on  the  right-hand  side  of  the  identity 
(B6)  has  the  explicit  form 


r  «o2 

t  2°2lV 


G[to>0;  (r02  +  r,2),/2]  -  exp[ 
or,  if  we  make  use  of  Eq.  (B12), 

G[w0;  (r02  +  r,2),/2]  -  exp^-  sin2(f  )j'  (B17) 

On  substituting  from  Eq.  (B17)  into  Eq.  (B6)  we  find  that 
G(to;;  r,)G(to)  -  to)0;  r0) 

(B,SI 

Finally,  on  substituting  from  Eq.  (B18)  into  the  expression 
(B5)  we  obtain  the  following  expression  for  the  spectrum  of 
the  scattered  field: 

S(m)(rS,  to))  -  N(r)H(9)to»4  exp[-  |  [- ^-]2}  •  (B19) 

where 
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<B20> 

The  quantities  u>(0)  and  f(0),  which  appear  on  the  right- 
hand  side  of  Eq.  (B19),  are  given  by  the  formulas  (B13)  and 
(B14),  respectively,  with  a(8)  given  by  Eq.  (B15).  The  fac¬ 
tor  N(r )  is  defined  by  the  formula  (B4). 


APPENDIX  C:  DERIVATION  OF  THE 
FORMULA  (4.10)  FOR  THE  FREQUENCY  AT 
WHICH  THE  SHIFTED  LINE  ATTAINS  ITS 
MAXIMUM 

The  formula  (4.7)  for  the  spectrum  of  the  scattered  field  may 
be  written  in  the  form 


S{m)(r,  8;  w)  -  N(r)H(8)f(u,  8),  (Cl) 

where 

<C2, 

with  u>(0)  and  f(0)  given  by  the  formulas  (4.8c)  and  (4.8d), 
respectively. 

To  determine  the  maximum  of  f(u,  0)  as  a  function  of  u 
(with  8  fixed),  we  differentiate  Eq.  (C2)  with  respect  to  oi. 
We  find  that 


dflu,  0) 

dm 


Obviously 

fc£>.  o 

do! 


(C3) 


when  <*!  *  0  or  when  «  **  oj±,  where 


4  - 


*  0, 


i.e.,  when 


<C4) 


If  we  make  use  of  Eqs.  (4.8c)  and  (4.8d),  it  follows  that 

— sH,*-'©’]T  « 

Both  of  these  roots  can  readily  be  shown  to  be  frequencies  at 
which  /( u>,  0)  and  consequently  S'“,(r,  8;  w)  attain  maximum 
values.  Only  the  one  with  the  positive  sign  is  relevant  be¬ 
cause  of  our  use  of  the  complex  analytic  signal  representa¬ 
tion  of  the  field12  (see  also  Ref.  17,  Sec.  10.2].  Denoting  this 
root  by  u>o'(0),  we  have 

""'w-  <c6’ 
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Hence,  for  fixed  0,  alt)  is  a  monotonically  increasing  function  of 
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ly  decreasing  function  of  this  ratio. 
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It  has  been  predicted  theoretically1*1  and  confirmed 
experimentally4*7  not  long  ago  that  correlations  in  the 
fluctuations  of  a  source  distribution  can  give  rise  to  fre¬ 
quency  shifts  of  lines  in  the  spectrum  of  the  emitted  radi¬ 
ation,  even  when  the  source  is  at  rest  relative  to  the  ob¬ 
server.  Because  of  the  close  analogy  that  exists  between 
the  processes  of  radiation  and  scattering,  a  similar  effect 
may  be  expected  to  occur  when  radiation  is  scattered  by 
a  “static"  medium,  i.e.,  a  medium  whose  constitutive  pa¬ 
rameters  (e  g.,  the  dielectric  susceptibility)  are  indepen¬ 
dent  of  time  but  are  random  functions  of  position,  with 
appropriate  correlation  properties.  Indeed,  a  very  recent 
theoretical  analysis1  has  shown  this  to  be  the  case. 

The  possibility  of  generating  frequency  shifts  of  spec¬ 
tral  lines  by  a  mechanism  other  than  the  motion  of  the 
source  relative  to  the  observer  or  by  gravitation  might  be 
of  particular  interest  for  astronomy,  especially  in  connec¬ 
tion  with  the  long-standing  controversy  surrounding  qua¬ 
sars  (see,  for  example,  Refs.  9-1 4).  However,  source 
correlations  and  correlations  produced  by  static  scatter¬ 
ing  do  not  imitate,  in  all  respects,  frequency  shifts  pro¬ 
duced  by  the  motion  of  a  source  relative  to  the  observer 
(i.e.,  Doppler  shifts).  The  relative  frequency  shifts15 


are,  in  general,  frequency  dependent  when  they  are  gen¬ 
erated  by  such  correlations,  whereas  they  art  frequency 
independent  when  they  are  manifestations  of  the  Doppler 
effect.  Moreover,  frequency  shifts  which  are  generated 
by  source  correlations  or  by  correlations  in  the  constitu¬ 
tive  parameters  of  static  scatterers  are  restricted  in  mag¬ 
nitude  to  values  that  are  of  the  order  of  or  smaller  than 
the  effective  widths  of  the  spectral  lines'*  Isee  the  discus¬ 
sion  following  Eq.  (10)  below]. 

In  the  present  Letter  we  investigate  whether  a  more 
general  correlation  mechanism,  involving  dynamic 
scattering,  could  give  rise  to  relative  frequency  shifts 
that  are  essentially  the  same  for  all  the  lines  in  the  spec¬ 
trum  of  radiation  originating  in  a  (stationary)  source 
and  whose  magnitude  could  be  arbitrarily  large. 

Consider  a  linearly  polarized  plane  electromagnetic 
wave  of  spectral  profile  S{,)(a)  that  propagates  in  the 
direction  specified  by  a  unit  vector  u,  incident  on  a  linear 
medium,  localized  in  space,  whose  dielectric  susceptibili¬ 
ty17  ij(r,r;»)  is  a  random  function  of  both  position  (r) 
and  time  (/).  We  assume  that  the  ensemble  which 
characterizes  the  statistical  behavior  of  the  medium  is 
homogeneous,  isotropic,  and  stationary,  at  least  in  the 
wide  sense.1*  It  has  recently  been  shown  that,  within  the 
accuracy  of  the  first  Born  approximation,  the  spectrum 
of  the  scattered  field  at  a  point  r—ra*  (|u'|  “1)  in  the 
far  zone  of  the  scatterer  is  given  by  '* 


S<“,(o>';ru\u)«Ma>“' 


to 

— U,<U  —W,(0 

c 


S(0(w)dat, 


(2) 
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correlation  function* 

<?(R,7» -<if*(M;a>)*(r+R,f +  7»> 


incident  radiation  five*  rise  to  oee  and  only  one  frequen- 
^  ey  w'  in  the  scattered  radiation  and  that,  moreover. 


of  the  dielectric  susceptibility  and  the  angular  brackets 
denote  the  average,  taken  over  the  ensemble  of  the  ran¬ 
dom  medium,  la  Eq.  (3)  V  represents  the  volume  of  the 
scatterer  (whose  linear  dimensions  are  assumed  to  be 
large  compared  to  the  correlation  distance  of  if),  w 
denotes  the  angle  between  the  direction  of  the  electric 
vector  of  the  incident  wave  and  the  direction  a'  of 
scattering,  and  c  is  the  speed  of  light  in  vacuo. 

For  our  purpose  it  is  convenient  to  rewrite  Eq.  (2)  in 
the  form 

S^Wa'.n) -Aa>'tJ~mX(o>\v»\m)SU)(»)da) , 

(6) 


where 


lf(e»',e);u',a)  J* 


a>;e> 


(7) 


We  will  refer  to  the  function  Xf(a>',a>;a',a)  as  the 
scattering  kernel. 

In  the  special  case  of  static  scattering  the  correlation 
function  C(R,T;m)  will  be  independent  of  the  temporal 
argument  T  and  we  will  then  denote  it  by  g(R;a>).  In 
this  case  Eqs.  (4)  and  (7)  imply  that  the  scattering  ker¬ 
nel  is  g'ven  by 


Let  as  now  tarn  to  the  more  general  case  of  a  medium 
whose  dielectric  susceptibility  varies  randomly  not  only 
la  space  but  also  in  time.  If  Sb>(aOBd(«-ato).  Eq. 
(6)  gives 

Sw(o'/i',»)  ■Aa»-Jlf(ai',a«B,,B) .  (11) 

This  formula  implies  that  a  frequency  component  ®  »o>o 
of  the  incident  radiation  will  give  rise  to  a  frequency 
component  a»‘  in  the  scattered  radiation  and  that,  in  gen¬ 
eral,  aVac  and,  moreover,  there  may  be  several  values 
of  a>’,  possibly  a  continuous  range  of  them,  associated 
with  any  particular  value  of  ate.  An  example  of  the 
latter  situation  is  provided  by  Briliouin  scattering  from  a 
simple  fluid  under  equilibrium  conditions.  The  scatter¬ 
ing  kernel  is  then  proportional  to  the  sum  of  three 
Lorentzian  distributions,  centered  at  frequencies  c/-c>o, 
ax>n  +  2(p/c)sin(0/2)],  and  atoll  **  2(c/c)sin(0/2», 
where  t>  is  the  speed  of  sound  in  the  fluid  and  6  is  the  an¬ 
gle  of  scattering.  Another  example  is  discussed  in  Ref. 
21. 

As  one  of  the  simplest  generalizations  of  the  formula 
(8)  to  media  which  vary  randomly  both  in  space  and  in 
time  let  us  suppose  that  the  scattering  kernel  has  the 
form 

#(<v’,e;u’.a)  ,  (12) 


Jf(o>\a);u’,u)—  g 


6(<o'  —  e>) , 


where 


where  a  is  independent  of  the  frequencies  but  may  de¬ 
pend  on  a  and  u'.  In  this  case  the  formula  (6)  gives  the 
following  expression  for  the  spectrum  of  the  scattered  ra¬ 
diation: 


(9) 


is  the  three-dimensional  spatial  Fourier  transform  of 
g(R,o>)  and  S  is  the  Dirac  delta  fuoction.  On  substitut¬ 
ing  from  Eq.  (8)  into  Eq.  (6)  we  see  that  the  spectrum  of 
the  scattered  radiation  is  given  by 


S'-Wu'.u)-/*®-* 


y  (o' -•);«’ 


S0)(  *').  (10) 


This  expression,  which  is  the  electromagnetic  analog  of 
the  main  result  of  Ref.  8  relating  to  scattering  of  scalar 
waves  by  static  media,  shows  that  spatial  correlation  of 
the  dielectric  susceptibility,  represented  by  the  two-point 
correlation  function  g(R,e>),  modifies  the  spectrum  of 
the  radiation  incident  on  the  scatterer.  However,  be¬ 
cause  the  influence  of  the  correlations  is  manifested  in 
Eq  (10)  only  through  a  multiplicative  factor,  do  oew 
frequency  components  are  generated  by  static  scattering 
li.e.,  S(“,(®Va',e)  "0  whenever  S®f,(w')“0l.  This 
conclusion  is,  of  oourse,  an  immediate  consequence  of  the 
presence  of  the  factor  d(«'-«)  in  the  expression  (8), 
which  implies  that  each  frequency  component  a>  of  the 


S(-Wd',d)  •Aa>«fl*\ow,;m\n)Su\o*>’) .  (13) 

Suppose  further  that  the  spectrum  Sb,(«>)  of  the  in¬ 
cident  radiation  consists  of  a  tingle  line  centered  on  the 
frequency  ®™£T  and  that  over  the  width  of  the  line  the 
factor  «a4/(®',e®';a,i’)  does  not  change  appreciably 
with  Then,  according  to  Eq.  (13)  the  spectrum 
S("^(o>W,e)  of  the  scattered  radiation  will  also  oonsist 
of  a  single  line,  but  this  line  will  be  centered  close  to  the 
frequency  •>'•£>',  where 

S' •S/a.  (14) 

On  substituting  from  Eq.  (14)  into  the  formula  (I)  we 
aee  that  the  spectral  line  has  been  shifted  in  frequency  by 
the  relative  amount 


_  S  —  S/a 
S/a 


•a  —  I . 


(15) 


Thus  apart  from  a  small  contribution  arising  from  the 
factor  »*/(•> ',0® the  relative  frequency  shift  is 
independent  of  the  central  frequency  S  of  the  spectral 
line  of  the  incident  radiation  and  of  the  width  of  the 
line 
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Evidently,  if  instead  of  a  single  line  the  spectrum  of 
the  incident  radiation  consisted  of  several  lines,  each 
would  be  shifted  by  essentially  the  same  amount 
z  “a  —  1,  thus  imitating  a  Doppler  shift  The  shift  will 
be  towards  the  lower  frequencies  (radshift)  if  a  >  I  and 
towards  higher  frequencies  (blueshift)  if  a  <  I  and  can, 
in  principle,  have  any  magnitude. 

Next  let  us  now  consider  a  more  general  kernel,  of  the 
form 

A’(a>',0;u<,u)  mf(w',w;u\m)e  t  (16) 


where  a  and  o  are  positive  quantities  that  are  indepen* 
dent  of  the  frequencies.  Suppose  that  the  spectrum  of 
the  incident  radiation  is  the  line 

SwU)-Jx“****,Vsr\  (17) 

where  B,  S,  and  T  arc  positive  constants.  On  substitute 
ing  from  Eqs.  (16)  and  (17)  into  the  formula  (6)  and  as¬ 
suming  that  the  factor  /(w’.srja'.a)  does  not  vary  appre¬ 
ciably  with  o>  over  the  effective  width  of  the  spectral  line 
(17),  we  obtain  for  S(,,)(a>’;ra',e)  the  expression 


S{“  W«‘.») -ABo>*/(w\S;u',*)  -u-t)'nr'du 


The  integral  in  Eq.  (18)  may  be  readily  evaluated 
(most  simply  by  making  use  of  the  to-called  product 
theorem  for  Gaussian  functions  (cf.  Appendix  A  of  Ref. 
S)]  and  one  then  finds  that 

s'-Hm's*:,  ■) -abH  2m) 

*  e-u-su)'nr*%  (I9) 

where 

r-!(ff»+rJ)l/J,  f-  -V+-rl  (20) 

«  o2  r2 J 

The  formula  (19)  shows  that  the  spectrum 
SimHai‘-.rn',u)  of  the  scattered  radiation  is  proportional 
to  the  product  of  the  factor  a>-/(a)',«;u’,u)  and  a 
Gaussian  function  centered  at  the  frequency  m'~S/a 
and  of  rms  width  T.  This  Gaussian  function  is  shifted 
with  respect  to  the  Gaussian  function  (17)  by  an  amount 
whose  relative  value  is  again  given  by  the  expression 
(IS)  and  is,  therefore,  independent  of  the  mean  frequen¬ 
cy  S  of  the  incident  radiation.  The  widths  of  the  two 
lines  are  not  the  same,  however.  Moreover,  the  factor 
a»‘4/(a>\57;u',u)  in  Eq.  (19)  will  produce  a  distortion  of 
the  line  and  possibly  an  additional  frequency  shift,  whose 
relative  value  may  or  may  not  be  small  compared  to  the 
value  x  ■*«  —  1,  depending  on  the  exact  form  of  the  func¬ 
tion  /(a>’a»;u’,u). 

It  is  evident  from  the  preceding  analysis  that  there  is  a 
possibility  that  scattering  from  some  media  whose  mac¬ 
roscopic  properties  fluctuate  randomly  both  in  space  and 
in  time  may  generate  frequency  shifts  of  spectral  lines 
which  closely  resemble  Doppler  shifts.  As  mentioned 
earlier,  this  possibility  might  be  of  particular  interest  in 
connection  with  the  long-standing  quasar  controversy. 
However,  in  order  to  apply  the  present  theory  to  quasar 
problems  it  would  be  necessary  to  have  a  reliable  model 
available  for  the  medium  through  which  radiation  origi¬ 
nating  in  these  astronomical  sources  passes,  and  no  mod¬ 
els  that  have  been  proposed  so  far  can  be  accepted  with 
full  confidence.  Moreover,  the  very  process  of  generation 
of  the  radiation  that  reaches  us  from  quasars  is  not 
currently  understood.  One  should  also  bear  in  mind  that 
throughout  the  preceding  analysis  the  incident  field  was 


assumed  to  be  a  plane  wave.  For  applications  to  quasars 
and  possibly  other  astronomical  objects,  additional 
averaging  over  a  range  of  directions  of  incidence  would 
also  have  to  be  performed.  Moreover,  extension  of  the 
analysis  beyond  the  first  Bom  approximation  and  also  to 
anisotropic  media  might  be  required.23 

There  is  a  widely  held  opinion  that  no  other  mecha¬ 
nism  except  the  Doppler  effect  and  gravitation  exist 
within  the  framework  of  present-day  physics,  which  can 
account  even  for  parts  of  the  redshifts  observed  in  the 
spectra  of  radiation  that  reaches  us  from  astronomical 
sources.  The  analysis  presented  in  this  Letter  raises  a 
question  about  the  correctness  of  this  opinion. 
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Experiments  that  illustrate  the  influence  of  source  oomelsuooi  oe  tbe  qpeetntm  of  light  detected  is  the  Ut  field  ere  reported  It 
■  demonstrated  tbet  source  eorrelsuoci  which  vioiitc  the  seal  ins  lew  [E  Wolf,  Phyv  Rev.  Lett.  S6  (19(6)  1370]  produce  ■ 
■ormtlusd  spectrum  id  the  far  field  that  is  different  from  the  normalised  spectrum  of  the  light  at  the  eouree. 


1.  Iatrodnction 

la  a  rectal  Letter  by  Wolf  [  1 J,  it  was  shown  that 
the  normalized  spectrum  of  light  is  the  same 
throughout  the  far  zone  and  is  equal  to  the  normal¬ 
ized  spectrum  of  light  at  the  source  provided  that  the 
complex  degree  of  spectral  coherence  is  a  function  of 
the  variable  k(p2-p, )  only,  where  (p2-pi )  denotes 
the  (vectorial)  distance  between  two  points  on  the 
source,  *«cu/c,  to  represents  the  angular  frequency 
and  c  it  the  speed  of  the  light.  If  the  complex  degree 
of  spectral  coherence  depends  on  frequency  only 
through  the  variable  k(p3-p, ),  the  source  is  said  to 
satisfy  the  scaling  law. 

It  is  well  known  that  the  correlation  properties  of 
light  can  be  controlled  by  application  of  the  Van  Cit- 
tert-Zemike  theorem  [2-6].  To  date,  using  conven¬ 
tional  sources  and  optical  systems,  tbe  secondary 
sources  with  controlled  correlation  that  have  been 
produced  all  satisfy  the  scaling  law. 

In  this  article  an  experimental  arrangement  that 
provides  a  secondary  source  with  a  complex  degree 
of  coherence  that  is  independent  of  the  wavenumber 
k  (Le.  violates  tbe  scaling  law)  is  described.  In  accord 
with  theoretical  predictions,  when  the  scaling  law  is 
violated  the  spectrum  in  the  far  zone  is  not  equal  to 
the  spectrum  of  light  at  the  secondary  source  and  it 
also  depends  on  tbe  observation  point.  The  experi¬ 
mental  measurements  of  the  spectrum  in  the  far  zone 
are  found  to  be  in  good  agreement  with  theory. 


2.  Coherence  la  linear  eptical  systems 

Consider  a  general  linear  system  with  an  impulse- 
response  function  given  by  h(x,  y,  (,  if,  to),  as  illus¬ 
trated  in  fig.  1.  With  a  given  transverse  component 
of  the  electric  field,  we  associate  a  complex  analytic 
signal  u(£,  it.  f)  and  define  the  spectral  amplitude 
V(i,  tr,  to)  to  be  tbe  time-truncated  Fourier  trans¬ 
form  of  w(£,  rr,  t). 

T 

UU.tr,  to) •:  J  uU.tr.t)  exp(-fiwr)  dr .  (1) 

-r 

Using  this  definition  one  can  express  tbe  spectral 
amplitude  of  the  system  output  l/n(x,  y,  to)  in  terms 
of  the  input  spectral  amplitude  t/,(£,  iy,  a>)  as  follows 

l/„(x,y,co) 

«  JJ  U,(£,t,w)  A(x,y,{,*«i0d£d»7.  (2) 


f 

1* 

A 

1 

Impulse 

\ 

nvhfwnw 

h(x.y;Ln;») 

1 

I  II 

Via  ‘  *; - 1 - iiti  f-irriimirj  it|nmhai  jwqnitei  rm  nr- 

*(x,r.L*  w). 
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In  eq.  (2)  the  spatial  coordinates  in  the  input  and 
output  planes  are  denoted  by  ({,  q)  and  (x,  y), 
respectively,  and  u»2xc/i,  where  2  is  the 
wavelength. 

When  one  considers  the  operation  of  an  optical 
system  with  broadband  illumination,  by  far  the  most 
convenient  quantity  to  calculate  is  the  cross-spectral 
density  function  defined  [7]  as 

*'(x,,y1;x2,yJ;cu) 


m.  lim 


r-ac 


<V(x,,yt\(o)  V*(X2.y»,to)'>t 
2  T 


(3) 


in  which  <...>,  denotes  an  average  over  an  ensemble 
of  sources.  Note  that  the  time-truncated  field  is  used 
since,  in  general,  the  noiselike  field  u(x,  y,  t)  is  not 
•quart  integrable. 

A  quantity  of  particular  interest  in  this  article  is 
the  spectral  intensity 

S( x,  y  (D)  •  W(x,  y  x.  y  u) ,  (4) 


which  represents  the  intensity  at  point  (x,  y )  and  fre¬ 
quency  to.  Note  that  the  total  intensity  at  point  (x, 
y)  is  obtained  simply  by  integrating  5(x,  y  w)  over 
all  temporal  frequencies. 

If  an  input  object  «,({,  7  u)  is  inserted  at  plane  I 
and  illuminated  by  a  field  with  spectral  amplitude 
7  <■>),  then  the  spectral  amplitude  {/,({,  7  w) 
kaving  plane  I  is  given  approximately  by 

U,({,  7  to)  =  7  w)  /,({,  7  &>) .  (5) 

Using  eqs.  (2),  (3)  and  (3),  the  cross-spectral  den¬ 
sity  function  H'„(x,,  y,;  x2,  y4;  to)  in  plane  II  is  writ¬ 
ten  as 

**n(  .  y. ;  .  y2 ;  w) «  J  J  J  7  w) 

X*(x, .  y, ;  {,  7  to)  h*(x2,  y2;  q';  to) 

Xd^dqctfdq'  .  (6) 

Eq.  (6)  is  a  general  expression  relating  the  cross- 
apectral  densities  in  planes  I  and  U. 

It  is  useful  to  define  a  normalized  form  of  the  cross- 
spectral  density  function  gives  as  follows: 


Au(a>) 


W',(*i,yl;x3,yi;ti>) 

lS(x,,y,;<u)  5(xJ.y1;a>))'/1  ’ 


(7) 


in  which  |^i](n>)  |  can  assume  values  in  the  range 
0<|A<ti(cv)i<l.  is  known  as  the  complex 

degree  of  spectral  coherence  (8]. 


3.  Caatrol  af  spatial  coherence 

Diagrams  of  the  optical  systems  used  in  the  exper¬ 
iments  are  shows  in  fig.  2.  In  both  system  configu¬ 
rations,  a  primary  source  illuminates  a  circular 
aperture  of  radius  a,  in  plane  I.  It  is  assumed  that 
illumination  is  provided  by  a  q ussi-homogeneous 
thermal  source  [6]  located  directly  behind  the  aper¬ 
ture  in  plane  1.  The  impulse-response  function  con¬ 
necting  planes  I  and  D  is  used  to  produce  a  secondary 
source  with  controlled  spatial  coherence  at  plane  II. 
The  spectrum  of  light  is  measured  at  the  secondary 
source,  plane  D,  and  in  the  far  field  of  the  secondary 
source,  plane  III 

A  theoretical  expression  for  the  complex  degree  of 
spectra]  coherence  in  the  secondary  source  plane  II 
can  be  calculated  using  eqs.  (6)  and  (7).  For  the  sys¬ 
tem  shown  in  fig.  2(a),  plane  I  is  located  in  the  from 
focal  plane  of  an  achromatic  telescope  objective  of 
focal  length  F.  Plane  D  is  located  in  the  back  focal 
plane  of  the  objective.  For  this  case  the  impulse- 


lb) 


Fig.  2.  Experimental  eosfigunttioas  for  reahzsuoe  af secondary 
sourer*  with  controlled  metis)  coherence  that  (a)  mtisfy  the 
acalim  law,  and  (b)  violate  the  scaling  law.  The  aperture  is  plane 
I  ■  {Humiliated  using  a  broadband,  usually  incoherent  source 
A  secondary  eource  with  a  m«ified  complex  degree  of  wcctral 
ooherence  it  formed  in  plane  0  through  application  of  the  Van 
Gmen-Zemike  theorem.  Measurement!  of  the  intensity 

are  made  is  the  mcoodary  source  plane  Q  and  ia  the  £ar  field  of 
Ac  mcoodary  source,  plane  ID. 
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response  function  connecting  planes  I  nnd  D  is  given 

by 

m(-ia>/2*cF)apH-ia>/cF)(x(+yri)] ;  (8) 

this  is  the  conventional  impulse  response  for  a  F-to- 
F  Fourier-trmnsfonn  system  [9]. 

In  the  system  arrangement  shown  in  fig.  2(b),  an 
achronutic-Fourier-transform  lens  [  10-12]  is 
inserted  between  planes  I  and  0.  With  spatially  coh¬ 
erent  light  the  effect  of  this  Fourier  achromat  is  to 
produce  an  optical  Fourier  transform  in  plane  II  in 
which  the  transform  size  is  independent  of  wave¬ 
length.  The  impulse  response  of  the  Fourier  achro¬ 
mat  is  given  by 

h*(x,  y,  {,  t.  «) 

«(-ki»o/2*cF)expl(-ktf0/cf)(J<+y>i)] ,  (9) 

in  which  m0  is  a  constant  and  corresponds  to  a  par¬ 
ticular  design  frequency. 

The  effect  of  the  impulse-response  functions  given 
byeqs.  (8)  and  (9)  is  illustrated  in  fig.  3;  these  pho¬ 
tographs  are  optica]  transform  patterns  recorded  in 
plane  II  of  the  systems  in  figs  2(a)  and  2(b),  respec¬ 
tively.  In  each  case,  a  400-/an-diameter,  circular 
aperture  is  inserted  in  plane  I  and  illuminated  with 
spatially  coherent,  broadband  light.  The  spectral 
bandwidth  of  the  light  is  200  nm.  In  fig.  3(a)  one 
notes  that  the  size  of  the  transform  pattern  scales  lin¬ 
early  with  the  wavelength  of  illumination.  The  scal¬ 
ing  of  the  pattern  with  wavelength  can  be  traced 
directly  to  the  frequency  dependence  of  A“Q"*(je,  y,  (, 
%  to)  in  eq.  (8).  The  system  in  fig.  2(a)  produces  a 
secondary  source  in  plane  II  that  obeys  the  scaling 
law.  In  fig.  3(b)  it  is  seen  that  the  scale  size  of  the 
optical  transform  pattern  produced  by  the  Fourier 
achromat  is  independent  of  the  illumination  wave¬ 
length.  The  Fourier  achromat  of  fig.  2(b)  produces  a 
secondary  source  in  plane  II  such  that  the  complex 
degree  of  spectral  coherence  is  independent  of  the 
illumination  wavelength,  and  therefore,  does  not  sat¬ 
isfy  the  scaling  law. 

Using  eqs.  (6)-(8)  gives  the  following  expres¬ 
sions  for  the  spectra]  intensity  and  complex  degree 
of  spectra]  coherence  in  plane  D  for  the  system  in  fig. 
2(a): 


SlT’(«)«ie»S<0*(«),  (10) 

Pir(»)-/,(xVz.  (11) 

where  S<0,(o>)  denotes  the  spectrum  at  the  primary 
source,  Jt(x)  is  the  Bessel  function  of  order  one, 

X«(wa,/cn{(*i-*i)J+(ya->'i)JJ,/a,  and  <c0  «  * 

constant  We  will  refer  to  sources  that  have  a  com¬ 
plex  degree  of  spectral  coherence  with  the  functional 
form  of  eq.  (1 1 )  as  Basel-correlated  sources. 

By  using  eqs.  (6),  (7)  and  (9)  one  can  calculate 
expressions  for  spectral  intensity  and  complex  degree 
of  spectral  coherence  in  plane  n  of  fig.  2(b),  and  one 
finds  that 

JTf(ft»)«r,  (<u0 tu)1  T(to)  S,e>(o) ) ,  (12) 

tf$(a>W.Cr')/X' .  (13) 

where  T(w)  is  introduced  to  account  for  transmis¬ 
sion  losses  through  the  Fourier  achromat, 
X'« (Woai/c^Kxj-x,  )*+ (y2-y,  )*) 1/J, and  jc,  is  a 
constant  Note  that  the  only  difference  in  the  expres¬ 
sions  for  the  complex  degree  of  spectral  coherence  in 
eqs.  (II)  and  (13)  is  that  to  in  eq.  (11)  has  been 
replaced  by  <o0  in  eq.  (13).  A  detailed  calculation  for 
the  spectra]  intensity  shows  that  Sf?  (r;  to)  is  inde¬ 
pendent  of  the  spatial  location  r  only  in  a  region 
whose  size  depends  on  the  spectra]  bandwidth.  In  the 
experiments  reported  below,  a  circular  aperture  is 
inserted  in  plane  D  to  ensure  that 
5ff(  r;  ci>)«$fP(ctf)forall  points  within  the  aperture. 


4.  Spectral  ef  light  la  the  hr  field 

The  spectrum  of  light  in  the  far  field  of  the  second¬ 
ary  source  located  in  plane  D  can  be  calculated  using 
the  above  linear  system  formalism  to  describe  the 
propagation  between  planes  D  and  HI  of  fig.  2.  The 
spectral  intensity  in  plane  III  is  readily  found  to  be 
given  by  the  following  formula: 

5, „(u,p,w) 

■  <Ti  (*!,*;») 

x  W„(x, ,  y, ;  x2,  y2;  co)h<m}(u,  r,  x„y,;  o>) 

X  A "•’*(«,  v.Xj.yj-.w)  dx,  6y,  dx2  d y2  .  (14) 
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F«  ).  Acy-dtU  diffraction  ptnm  recorded  ia  the  back  focal  plane,  plane  n,  of  (a)  aa  achromatic  tefctcopc  objective  km.  and  (b)  as 
achroma  tic-Fourier-transfi one  leas  Tbcae  photographs  «m  produced  by  ffiumiaatiaf  a  circular  aperture  located  ta  plane  I  with  broad- 
hud.  apaually  oobercet  h*ht.  la  (a)  the  diffrnctioo  patten  teaks  haearty  with  winleafth,  tad  m  (b)  the  tbt  of  the  diffractioa  patten 
It  independent  of  the  wavelength. 
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when  i„(x,  y.  to)  represents  a  transmission  function 
in  plane  11,  and 

h,m>(u,  v; x,  y\  to) « ( - io»/2*«) 
xexp[( +io»/2cr)(uJ  -f  p*)] 

Xexpt(— i<i>/«)(iur+iy)] .  (15) 

Taking  t„(x,  y,  o>)«  1  and  using  eqs.  (7),  (10), 
(11)  and  (14),  one  finds  that 


m 


(16) 


SSTv(r.w)«icSfrv(w).  fKtttX/F, 

*0,  r>a,*//\ 

when  r«(uJ-fv,),/2and  k2  is  a  constant. 

Using  eqs.  (7)  and  (12)— (14 )  one  obtains  the  fol¬ 
lowing  expression  for  the  spectrum  of  the  field  in  the 
far  zone  of  the  achromatized  secondary  source  of  fig. 
2(b): 

sj/!(r,  to)mKi(to/<o0)4St![((o) ,  r<w0fl|i/wF, 
«0,  r>to0a,z/(oF , 

(17) 

when  1C3  is  a  constant 

S.  Experimental  resalts 

In  the  experiments  a  100-watt  tungsten-halogen 
lamp  is  used  as  the  primary  source  and  is  located  in 
dose  proximity  behind  a  circular  apertun  of  radius 
Oi*2.5  mm  in  plane  1  (see  fig.  2].  The  focal  length 
Fit  180  mm.  Using  these  values  the  cohennce  inter¬ 
val  at  the  secondary  source  plane  11  is  24  ym  at 
2*20*330  nm. 

A  circular  apertun  of  radius  o2*0.S  mm  is  located 
in  plane  II.  The  far-field  condition  is  satisfied  when 
x»ojj/2*0.S  m  at  2*  500  nm.  To  obtain  an  accept¬ 
able  SNR  in  our  detection  system,  the  distance  z  is 
taken  to  be  one  meter.  The  apertun  in  plane  II  is  also 
used  to  ensun  that  the  spectrum  of  the  secondary 
aource  in  the  experimental  arrangement  of  fig.  2(b) 
is  independent  of  the  spatial  location  within  the 
apertun  (see  eq.  (12)). 

The  spectral  intensity  of  the  light  at  the  various 
observation  points  is  measured  using  a  grating 
monochromator. 

Fig.  4(a)  shows  curves  for  the  spectral  intensity 


Fig  4  Normalized  spectral  intensity  at  the  secondary  source  plant 
U  and  in  the  far  field  of  the  secondary  source,  plant  III.  for  (a  1  a 
source  that  obeys  the  scahnt  la*,  and  (b)  a  source  in  which  the 
complex  decree  of  spectral  coherence  is  independent  of  the  illu- 
mination  frequency  (i.e.,  violates  the  acalini  law).  In  (a)  the 
spectral  intensity  is  normalized  to  the  peai  value  of  the  spectral 
intensity  in  the  respective  planes.  In  (b)  the  spectral  intensity  in 
planes  II  and  HI  art  scaled  so  that  5Ti"  wSf/l  (0,  u0)  ■  1 


measured  at  the  source  plane  11  and  at  various  loca¬ 
tions  in  the  far  field,  plane  III,  for  the  system  config¬ 
uration  of  fig.  2(a).  As  expected,  Bince  this  source 
obeys  the  scaling  law,  the  normalized  spectrum  is 
indeed  invariant  on  propagation. 

In  the  experiments  involving  the  system  configu¬ 
ration  of  fig  2(b),  an  all-glass  Fourier  achromat  [11] 
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Fig  5.  Matured  values  of  Ihe  aontulued  spectral  iaientiiy  in 
the  far  field,  plane  III,  for  cate  (b)  (i.e.  when  the  source  does  not 
obey  the  scaling  law).  The  (pectraJ  intensity  is  measured  at  off- 
axis  angles 6m 0‘,  10*.  IS*.  20*.  Note  thai  the  spectral  intensity 
versus  w  narrows  as  the  off-axis  angle  *  increases.  The  intensity 
dau  arc  normalized  with  respect  to  SmfO.oi,).  where  u,« 
3.2x10'*  s'1. 


is  used.  The  lens  consists  of  three  widely  spaced  lens 
groups.  The  lens  is  capable  of  forming  250  adjacent, 
achromatized  Airy  disks  across  a  7-mm  diameter  in 
the  optical  transform  plane.  Achromatization  is  well 
corrected  for  wavelengths  ranging  from  approxi¬ 
mately  470  did  to  590  nm  (4.0xl0IJ>fc)>3.2x  1015 

The  spectral  intensity  measured  on-axis  at  planes 
II  and  III  for  the  Fourier-achromat-system  configu¬ 
ration  is  shown  in  fig.  4(b)  -  the  solid  curves  are 
plotted  from  laboratory  measurements  and  the 
dashed  curve  is  the  theoretical  prediction  for  the 
spectral  intensity  as  given  in  eq.  (17). 

In  fig.  5,  measurements  of  the  spectral  intensity  at 
various  ofT-axis  points  in  plane  Ill  of  fig.  2(b)  are 
given.  Note  that  the  width  of  the  spectral  intensity 
versus  to  narrows  as  the  angle  of  observation 
increases.  For  the  off-axis  points,  the  measured  val¬ 
ues  of  the  spectral  intensity  are  somewhat  lower  than 
that  predicted  by  eq.  (17).  This  can  be  explained  by 
the  fact  that  the  far-fieid  condition  is  only  approxi¬ 
mately  satisfied.  For  the  calculation  of  the  spectral 
intensity  at  off-axis  points  in  plane  III,  one  should 
actually  use  the  impulse  response  that  corresponds 
to  a  Fresnel-diff faction  geometry.  However,  for  the 


sake  of  brevity,  this  calculation  is  omitted  in  the 
present  treatment. 


i  Discassion 

As  observed  in  the  experiments  and  consistent  with 
Wolf's  theoretical  predictions  1 1),  the  spectrum  of 
light  in  the  far  field  depends  on  the  correlation  prop¬ 
erties  of  the  light  at  the  source.  If  the  optical  field  at 
the  source  obeys  the  scaling  law,  the  spectrum  of  light 
on  propagation  is  invariant.  Departures  from  the 
scaling  law  produce  changes  in  the  spectrum  that 
depend  on  the  correlation  properties  of  the  source  and 
on  the  location  of  the  observation  point. 

In  our  experiments  a  circular  aperture  is  inserted 
in  plane  I  of  fig.  2  to  produce  a  secondary  source  at 
plane  II  that  is  Bessel-correlated  (see  eqs.  (11)  and 
(13)].  A  Fourier  achromat  (see  fig.  2(b))  is  used  to 
produce  a  field  correlation  in  plane  II  that  is  inde¬ 
pendent  of  the  illumination  frequency  (a  departure 
from  the  scaling  law).  In  this  case  the  spectrum  in 
the  far  field  is  different  from  the  spectrum  measured 
at  the  source  and  changes  for  different  observation 
points  in  the  far  field  (see  eq.  (17)  and  figs.  4  and 
5] .  While  a  Bessel-correlated  source  was  used  in  these 
experiments,  it  is  noted  that  virtually  any  allowed 
functional  form  for  the  complex  degree  of  spectral 
coherence  can  be  produced  through  application  of  the 
Van  Cmert-Zernike  theorem.  For  example,  a  gaus- 
Sian -correlated  source  can  be  produced  by  inserting 
a  gaussian  transmission  function  in  plane  1.  A  Four¬ 
ier  achromat  inserted  between  planes  I  and  II  elimi¬ 
nates  the  wavelength  dependence  in  the  complex 
degree  of  spectral  coherence  in  plane  II.  A  different 
frequency  dependence  for  the  complex  degree  of 
spectral  coherence  oould  be  produced  by  a  modifi¬ 
cation  of  the  design  rules  used  for  the  Fourier 
achromat. 

The  effects  of  source  correlation  on  the  spectrum 
is  to  be  distinguished  from  changes  in  the  spectrum 
due  to  diffraction  of  fully  or  partially  coherent  light 
by  an  aperture  located  in  the  source  plane.  Although 
both  effects  can  modify  the  spectrum  of  light  on 
propagation,  their  physical  origin  is  different.  For 
example,  suppose  that  one  placed  a  diffraction  grat¬ 
ing  in  the  secondary-source  plane  D.  Using  eq.  (14) 
and  the  appropriate  expression  for  the  transmission 
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ftmction  tD(xy,a>)  of  s  diffraction  grtting,  one  can 
write  aa  expression  for  the  spectrum  Sm(u,  v;  to)  in 
the  far  field.  The  transmission  function  ln(x,  y,  to) 
will  affect  Sm(  u,  v\  to)  only  if  the  coherence  interval, 
associated  with  the  cross-spectr*]  density  Wu(.x,,  y,; 
jrJ(  y,;  to),  is  on  the  order  of  or  greater  than  the  grat¬ 
ing  period,  i.e.  the  light  is  spatially  coherent  over  at 
least  a  portion  of  the  grating  structure.  Similiar  rea¬ 
soning  can  be  used  for  the  case  in  which  a  refractive 
dement,  such  as  a  prism,  is  placed  in  the  secondary- 
source  plane. 

In  our  experiments,  a  1 -mm -diameter  circular 
aperture  is  inserted  in  plane  II.  The  purpose  of  this 
aperture  is  to  ensure  that  the  far-field  condition  is 
satisfied  (approximately)  in  the  measurements  of 
$,„(«,  v,  to).  As  stated  above,  the  coherence  interval 
at  plane  11  is  24  /an.  Since  the  coherence  interval  is 
much  less  than  the  object  structure  (aperture  diam¬ 
eter),  the  illumination  in  plane  II  is  essentially 
incoherent  and  the  aperture  does  not  modify  the 
spectrum  of  light  in  the  far  field.  In  our  experiments, 
changes  in  the  spectrum  on  propagation  are  due  to 
the  correlation  properties  of  the  secondary  source. 
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Wolf  has  recently  shown  that  the  spectrum  of  radiation  from  an  extended  source  changes  on  propaga¬ 
tion  unless  a  certain  scaling  condition  is  obeyed  by  the  degree  of  spectral  coherence  across  the  source 
For  a  large  class  of  source-coherence  functions,  the  change  may  be  such  as  to  produce  red  shifts  or  blue 
shifts  of  spectral  lines.  We  have  performed  an  acoustic  experiment  with  two  small  partially  correlated 
sources  and  demonstrated  Wolfs  prediction  of  frequency  shifts  of  spectral  lines  by  this  mechanism. 

FACS  numbers  42.68  Hf.  42.IO.Mg,  43.20.+g 


Wolf'  has  shown  theoretically  that  the  spectral  density 
of  radiation  in  the  far  zone  of  an  extended  source  may  be 
different  from  the  spectral  density  of  the  individual  radi¬ 
ating  source  points  if  there  is  some  spatial  coherence  be¬ 
tween  the  source  fluctuations  at  different  points  in  the 
source  region.  Morris  and  Faklis  subsequently  confirmed 
experimentally  the  noninvariance  of  the  spectrum  for  ra¬ 
diation  from  a  particular  optical  source.2  In  later  publi¬ 
cations  Wolf3  4  has  shown  that  the  center  frequency  of  a 
spectral  line  could  be  shifted  from  the  center  frequency 
of  the  same  line  as  seen  at  the  source.  This  shift  could 
be  either  a  red  shift  or  a  blue  shift  depending  on  the 
specific  form  of  the  source-correlation  function. 

Wolf s  main  result  is  the  following  expression  for  the 
spectral  density  of  the  field  at  frequency  to,  at  a  point  R u 
<u2  1 )  in  the  far  zone,  emitted  by  sources  of  a  wide 

class: 

St/(J?u)«Sg(<u)>Jp(ku,<u)  (1) 

Here  Sq((o)  is  the  spectral  density  at  each  source  point, 
pQ(ku,<i>)  is  the  spatial  Fourier  transform  of  the  degree 
of  spectral  coherence  pq(i\  —  r2,<u)  of  the  source  fluctua¬ 
tions  at  two  typical  source  points,  and  k—to/c,  c  being 
the  velocity  of  propagation  of  the  radiation.  Formula 
(l)  shows  that,  in  general,  the  spectrum  of  the  far  field 
will  not  be  equal  to  the  source  spectrum  Sq((o )  but  will 
be  modified  by  the  source  correlations. 

Wolfs  analysis  was  based  on  correlation  theory  of  sca¬ 
lar  fields  and  hence  must  also  apply  to  acoustic  fields. 
To  test  experimentally  the  dependence  of  the  spectrum 
of  the  emitted  field  on  the  correlation  properties  of  the 
source,  we  considered  a  very  simple  radiating  system, 
consisting  of  two  small  correlated  acoustical  sources  with 
the  same  source  spectrum  Sq(o>)  and  we  measured  the 
spectrum  Su(to)  of  the  emitted  field  in  the  far  zone.  The 
correlation  between  the  two  sources  can  be  characterized 
by  a  correlation  coefficient  pq((o)  which  is  a  function  of 


the  frequency.  For  our  case  one  finds  that  Eq  (1) 
reduces  to  the  following  expression  for  the  spectral  inten¬ 
sity  Sc/ (to)  at  the  observation  point: 

St/(a>)*5g(fl>)[l  +Re/ip(w)].  (2) 

This  formula  is  derived  and  pq(w)  is  defined  in  the  ac¬ 
companying  Letter  [Wolf,5  Eq.  (8)1,  where  it  is  also 
shown  that  the  factor  Il  +  Re^g(w)l  in  Eq.  (2)  can  pro¬ 
duce  a  field  spectrum  Su(to )  whose  central  frequency  is 
shifted  either  towards  the  lower  frequencies  (producing  a 
red  shift)  or  towards  the  higher  frequencies  (producing  a 
blue  shift),  compared  with  the  central  frequency  of  the 
source  spectrum  Sg(co). 

We  have  demonstrated  this  “Wolf  shift”  with  a  two- 
point  acoustic  source  which  has  a  frequency-dependent 
cross  correlation  of  the  two  source  points.  We  have  ob¬ 
served  both  red  shifts  and  blue  shifts. 

Our  experiment  is  illustrated  in  Fig.  1.  Two  loud¬ 
speakers  of  diameter  </«  9  cm  are  separated  by  a  dis¬ 
tance  D™40  cm  and  located  at  distances  R\  and  /?: 
from  a  microphone  with  an  aperture  of  1.5  cm.  We 
chose  m.  None  of  these  dimensions  ap¬ 

peared  to  be  critical  and  were  not  expected  to  be.  An 
acoustic  spectral  line,  i.e.,  band-limited  random  noise, 
with  center  frequency  n>o  (wavelength  ko)  and  width  Aco 
was  broadcast  from  each  speaker. 

In  Fig.  2  we  show  the  detected  spectral  density  of  the 
field  from  source  1,  labeled  A,  and  source  2,  labeled  B 
These  spectra  were  measured  at  the  microphone  one  at  a 
time,  with  the  other  source  switched  off,  and  had  the 
same  shape  as  the  spectra  at  the  source  which  we  verified 
by  also  measuring  the  spectrum  of  the  voltage  wave  form 
which  drove  each  speaker.  The  intensities  of  the  two 
sources  at  each  frequency  were  equal  to  within  5%;  the 
center  frequencies  were  both  at  1180  Hz  and  the 
linewidth  at  half  maximum  amplitude  was  —300  Hz. 
The  spectral  density  of  the  sound  at  the  microphone  with 
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Source  1  Source  2 


Detector 


FIG.  I.  The  configuration  of  our  experiment.  Two  acoustic 
sources  each  with  diameter  d  are  located  distance  D  apart. 
The  detector  is  a  microphone  at  distances  and  Ri  from  the 
two  sources. 

both  sources  switched  on  is  shown  by  the  curve  labeled  C 
in  Fig  2.  It  has  a  peak  at  1020  Hz  exhibiting  a  red  shift 
of  160  Hz.  We  will  now  explain  how  this  was  achieved. 

Although  both  sources  radiated  band-limited  random 
noise  with  the  same  spectral  density,  they  were,  in  fact, 
partially  correlated.  Each  spectrum  had  two  parts,  an 
anticorrelated  component,  a  broad  line  centered  at  1180 
Hz,  and  a  correlated  component  centered  at  1020  Hz. 
The  anticorrelated  component  of  each  line  was  prepared 
by  modulating  the  phase  of  a  Hewlett-Packard  3325A 
frequency  synthesizer  with  low-pass  filtered  white  noise 
which  created  a  300-Hz-wide  spectral  line  centered  at 
the  synthesizer’s  set  frequency.  This  signal  was  passed 
through  a  L-C  notch  filter  set  at  1020  Hz  then  split  into 
two  equal  parts,  and  one  part  was  inverted  to  yield  two 
anticorrelated  random  signals  xO)  and  — x(r).  The 
correlated  components  were  generated  by  passing  the 
output  of  a  second  Hewlett-Packard  3325A  synthesizer, 
which  was  set  to  the  same  frequency  as  the  first  and 
phase  modulated  by  an  independent  source  of  random 
noise,  through  a  L-C  band-pass  filter  tuned  to  the  same 
frequency  as  the  notch  filter.  This  signal  was  divided 
into  two  equal  parts,  each  denoted  by  yO),  and  added  to 
the  uncorrelated  signals  to  give  the  two  signals  x(f) 
+y(t)  and  -  xU)+y(,i)  which  drove  the  speakers. 

Assuming  x(t)  and  y(t)  to  be  independent  random 
processes,  one  can  calculate  the  degree  of  spectral  coher¬ 
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FIG.  2.  Square  root  of  the  spectral  density  of  the  acoustic 
spectral  lines  in  three  cases.  The  curves  labeled  A  and  B  are 
the  results  of  measurements  at  the  microphone  when  one 
source  at  a  time  is  switched  on.  The  curve  labeled  C  represents 
the  red-shifted  acoustic  spectral  line  when  both  sources  are 
switched  on  to  act  as  a  correlated  extended  source. 

ence  (jq(iu)  and  we  find  that 

Mq(o>)  *"l2S>(£u)/Sp(£u)]  —  I,  (3) 

where  is  the  spectrum  of  y(t).  Using  Eq.  (2)  we 

find  the  spectral  density  at  the  observation  point  is 

Su(to)*Sy(a>)%  (4) 

i.e.,  it  is  proportional  to  the  spectrum  of  the  correlated 
components  of  the  signals  which  drove  the  two  speakers. 

We  also  produced  a  blue  shift  of  a  spectral  line  by 
moving  the  center  frequency  of  the  anticorrelated  com¬ 
ponent  of  the  source  spectra  to  840  Hz  keeping  the 
correlated  components  the  same.  With  both  sources 


FIG  3.  The  curves  labeled  A  and  B  have  the  same  meaning 
as  in  Fig  2.  The  curve  labeled  C  now  represents  a  blue-shifted 
line 
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switched  on,  the  peak  in  the  spectrum  at  the  detector  was 
at  1000  Hz,  thus  exhibiting  a  I60-Hz  blue  shift,  as  seen 
in  Fig.  3. 

There  were  a  few  other  features  in  the  spectra  which 
we  now  wish  to  explain.  The  weak  line  present  at  300 
Hz  and  the  feature  at  very  low  frequency  were  both 
direct  electrical  pickup  from  the  power  lines  and  had  no 
effect  on  the  experiment.  The  jagged  features  in  the 
spectra  in  Fig.  2  seem  to  be  an  effect  of  the  room  in 
which  the  experiments  were  performed.  The  rectangular 
room  was  not  totally  anechoic  and  the  room  resonances 
made  it  impossible  to  simultaneously  obtain  smooth 
profiles  for  all  the  lines. 

In  summary,  we  have  confirmed  experimentally  the 
theoretical  prediction  that  lines  in  the  spectrum  of  radia¬ 
tion  from  a  partially  coherent  source  may  indeed  be 
shifted  with  respect  to  their  locations  in  the  emitted 
spectra  of  the  individual  source  points.  The  significance 
of  this  effect  goes,  however,  beyond  the  present  illustra¬ 
tion  of  principle.  It  demonstrates  a  new  mechanism  by 
which  spectral  lines  from  extended  sources  may  be  fre¬ 
quency  shifted.  This  effect  may  be  important  in  astro¬ 
physics.  The  effect  of  spatial  correlations  between  emit¬ 
ting  elements  of  an  extended  astronomical  source  have 
not  previously  been  considered  (with  the  possible  excep¬ 
tion  of  astronomical  masers).  Specifically,  the  cosmolog¬ 
ical  interpretation  of  the  quasar  red  shifts  has  been  ques¬ 


tioned4,7  and  Wolf  has  suggested  that  this  new  mecha¬ 
nism  may  provide  contributions  to  the  red  shifts  of 
quasers.3  However,  there  has  not  yet  been  offered  any 
firm  suggestion  for  the  origin  of  possible  correlations 
among  radiating  atoms  in  astronomical  sources  but  the 
correlation  lengths  required  for  producing  observable  red 
shifts  need  only  be  of  the  order  of  a  wavelength,3  which 
seems  a  relatively  modest  requirement. 
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Frequency  shifts  of  the  optical  spectra  emitted  by  a  pair  of  suitably  correlated  small  sources  have  been  recently  predicted  by 
Wolf  and  demonstrated  experimentally  in  the  acoustic  domain  by  Bocko.  Douglass  and  Knox  We  sho*  that  a  small  modification 
in  the  form  of  the  correlation  functions  used  by  Wolf  leads  to  a  model  of  the  correlated  pair  that  can  be  easily  synthesiied  in  the 
optical  domain.  The  results  of  experimental  tests  exhibiting  redshtfts  and  blueshifts  are  presented. 


The  existence  of  spatial  correlations  within  a  source 
afTects  the  way  in  which  the  spectral  properties  of  the 
light  emitted  by  the  source  propagate  (1-4}.  Fur¬ 
thermore,  it  has  been  proved  by  Wolf  that  a  suitably 
correlated  source  can  give  rise  to  redshifts  or  blue- 
shifts  of  the  emitted  light  with  respect  to  the  case  of 
an  uncotTelaled  source  [5,6].  In  view  of  the  impor¬ 
tance  of  this  prediction,  in  particular  for  astronom¬ 
ical  problems,  experimental  proofs  of  it  are  very 
significant.  In  order  to  facilitate  experiments,  Wolf 
offered  w  hat  is  perhaps  the  simplest  model  of  a  source 
exhibiting  redshifts  and  blueshifts  of  the  above  type 
(7).  It  consists  of  a  pair  of  small  sources  suitably 
correlated.  Bocko,  Douglass  and  Knox  (8}  suc¬ 
ceeded  in  s>  nthesizing  an  acoustic  source  of  this  type 
and  found  a  complete  agreement  with  the  theoretical 
predictions. 

In  spile  of  its  simplicity,  the  form  of  the  correla¬ 
tion  functions  used  in  Wolfs  paper  is  not  very  easy 
to  synthesize  in  the  optical  domain.  To  overcome 
this,  we  propose  a  small  modification  of  the  form  of 
the  correlation  functions  to  be  used  in  the  source 
model.  We  will  show  that  in  this  way  redshifts  and 
blueshifts  can  still  be  exhibited  while  making  the  ex¬ 


perimental  synthesis  of  the  source  quite  easy.  We  shall 
begin  by  recalling  a  few  essential  results  from  the  pa¬ 
per  of  Wolf  (7J.  Two  small  fluctuating  sources,  par¬ 
tially  correlated,  are  located  at  certain  points,  say  P, 
and  P2  The  space-frequency  description  of  the  cor¬ 
relation  between  the  source  fluctuations  is  given  by 
the  cross-spectral  density  [9]  H’0(P,,  P2,  w).  where 
ai  is  the  angular  frequency.  The  normalized  version 
of  the  cross-spectral  density  is  the  degree  of  spectral 
coherence  (9] 

JtQ (Pi <  Pj.  w)  =  H q (P,, P2,  &y)/SQ (a>) ,  (1) 

where  SQ( w)  is  the  spectrum  (assumed  to  be  the 
same)  of  each  of  the  two  source  distributions  Sup¬ 
pose  that  the  spectrum  of  the  field  produced  by  the 
two  sources  is  observed  at  any  point  P  on  the  per¬ 
pendicular  bisector  of  the  line  joining  P,  and  Pj.  Wolf 
showed  that  such  a  spectrum,  say  Su(P,  cu),  is  given 
by 

5u(P,tu)S!(2/RJ)SQ(tu)  [l  +  Re/rQ(o>)]  .  (2) 

where  Re  denotes  the  real  pari  and  R  is  the  common 
distance  of  the  point  P  from  the  two  sources.  Here¬ 
after,  the  explicit  dependence  of  both  fjQ  and  H'0  on 
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P,  and  P:  is  dropped.  Eq.(2)  clearly  shows  that  the 
field  spectrum  Sv  at  the  point  P  is  different  from  the 
source  spectrum  5Q  unless  the  degree  of  spectral  co¬ 
herence  is  independent  from  u  (in  particular,  for 
either  uncorrelated  sources,  i.e.  for  nQ{(o)  =  0,  or 
perfectly  correlated  sources,  i.e.  for  jiQ(cd)=:  1 ).  In 
certain  cases,  the  field  spectrum  can  be  redshifted  or 
biueshifted  with  respect  to  the  spectrum  of  the  source 
fluctuations.  This  can  be  seen  through  a  suitable 
choice  of  H'Q(o>)  and  SQ(o>).  We  shall  first  consider 
the  form  chosen  by  W'olf  in  ref.  [7],  In  order  to  fa¬ 
cilitate  the  comparison  with  the  modified  form  to  be 
seen  shortly,  we  shall  quote  the  Wolfs  formulas  in 
the  following  manner,  slightly  different  (but  of  course 
equivalent)  to  the  original  one 

If  q  (tu) =B  exp[  -  (aj-u'o  )3/2Sq  ] 

—A  exp[-  (cu-cj0)j/2<5o]  .  (3) 

Sq(co)=A  exp[  —  (cu— too  )j/2<5q],  (4) 

where  A ,  £u0,  cuo,  So  and  S'0  are  all  positive  constants 
with  the  same  meaning  as  in  ref.  [7],  B  is  also  a  po¬ 
sitive  constant  and  must  satisfy  the  inequality 
B/A<2  (see  ref.  [ 7 J  for  further  details).  On  in¬ 
serting  from  eqs.  (3)  and  (4)  into  eqs.(l )  and  (2), 
the  following  expression  is  obtained  for  the  spectrum 
of  the  field  at  the  point  P 

Sv(?,oj)  =  (2B/R3)txv\-(u-w'0)3/262).  (5) 

Hence,  a  redshift  is  produced  if  w'0  <  co0  and  a  blue- 
shift  in  the  opposite  case. 

We  now  modify  If ’Q  and  SQ  in  the  following  way 


Wq(o ))=Aa  exp[  -  (lD-(Da  )3/263a  ] 

- Ag  t\v[  -  {u-w,)3 /26]]  , 

(6) 

SQ(iu)-Aa  cxp[ -  (w-cja  )3/263a  ] 

+A,  exp[  -  (lo-cjg)3 /26f  ]  , 

(7) 

where  Aa,  At,  u>a ,  cofi,  ia  and  St  are  positive  con¬ 
stants.  For  the  moment,  no  constraint  is  imposed  on 
A0IAp  On  comparing  eqs.  (6),  (7)  with  eqs. (3),  (4), 
we  see  that  the  present  case  differs  from  the  one  con¬ 
sidered  by  Wolf  mainly  for  the  form  of  Sq(cu).  How¬ 
ever,  see  eq.(l),  the  normalized  version  of  the 
correlation  function  /cQ(o;)  is  also  modified.  Pro¬ 
ceeding  as  before,  we  easily  obtain 

Su(o>)  =  (4/l0/KJ)exp(-(tu-a>0)J/2<5£]  ■  (8) 


In  order  to  see  that  eq.(8)  can  entail  a  redshift  or  a 
blueshift,  let  us  refer  to  the  case  Aa=At,  6a-6g  and 
\ti)a-w,\  <26a.  In  this  case,  the  overlapping  of  the 
two  gaussian  profiles  given  by  eq.(7)  gives  rise  to  a 
curve  in  which  the  two  component  lines  are  not  re¬ 
solved  and  the  source  spectrum  has  a  maximum  at 
Wm=  {ioa+U)t)l2.  Accordingly,  the  field  spectrum, 
see  eq.(8),  is  redshifted  if  ua<(Ou  and  biueshifted 
in  the  opposite  case.  An  example  is  shown  in  fig.  1, 
where  the  numerical  values  are  approximately  the 
same  as  in  the  experiment  to  follow.  Furthermore,  it 
can  be  easily  seen  that  the  shift  phenomena  appear 
even  if  the  conditions  on  the  parameters  of  eqs.  (6) 
and  ( 7 )  are  somewhat  relaxed. 

The  main  virtue  of  the  present  source  model,  as 
compared  to  the  previous  one,  is  that  an  experimen¬ 
tal  synthesis  of  it  is  easier  to  produce. 

In  order  to  see  this,  we  begin  by  noting  that  both 
the  cross-spectral  density  and  the  power  spectrum 
described  by  eqs. (6),  (7)  are  the  sum  of  two  terms 
with  weighting  factors  Aa  and  Af.  Were  Ag- 0,  the 
remaining  terms  in  eqs. (6)  and  (7)  would  represent 
perfectly  correlated  source  fluctuations  at  P,  and  P2 
with  a  gaussian  power  spectrum  centered  at  coa.  In 
fact  on  evaluating  the  degree  of  spectral  coherence 
we  would  find  Mq(oj)  *  I.  On  the  other  hand,  were 


Fig  1 .  Redshift  *nd  blueshift  produced  by  two  sources  correlated 
as  in  eqs  (6)  and  (7).  Curve  (a):  spectrum  obtained  when  the 
sources  ndiate  simultaneously  with  the  follow in|  parameters: 
cu.w3.0l  x  101’  rad/s;  tu,«3.04x  10”  rad/s;  d.«<5,w2.5x  I0” 
rad/s ,AamA,=  I  (arbitrary  units)  Curve  (b ):  same  as  curve  (a) 
except  that  cu.w3.04x  10”  rad/s  and  tUjwJOIxlO1’  rad/s. 
Curve  (c):  field  spectrum  (multiplied  by  a  factor  2)  produced  at 
the  observation  point  by  a  single  source. 
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AgsO,  we  would  find  j/q (a/)*  - 1,  evidentiating  that 
the  second  term  alone  in  eqs.(6)  and  (7)  accounts 
for  perfectly  anticorrelated  fluctuations  of  the 
sources.  For  Aa*  0  and  Ae+  0,  the  source  correlation 
has  both  a  correlated  and  an  anticorrelated  compo¬ 
nent.  These  two  components  are  mutually  indepen¬ 
dent  as  implied  by  the  fact  that  they  are  summed 
together  in  cqs.(6)  and  (7). 

Starting  from  the  previous  interpretation  it  is  not 
difficult  to  devise  possible  schemes  for  synthesizing 
two  small  sources  describable  by  eqs.(6)  and  (7). 
As  an  example,  let  us  refer  to  the  scheme  of  fig.  2. 
The  two  sources  are  realized  as  small  holes,  centered 
at  P,  and  P2,  pierced  in  an  opaque  mask  M.  Each 
hole  is  illuminated  by  two  different  beams  of  light  by 
means  of  a  beam-splitter  BS.  One  beam,  say  Ba,  pro¬ 
vides  a  spatially  coherent  and  cophasal  illumination 
at  P,  and  P2.  Its  power  spectrum  is  described  by 
eq.(7)  letting^/=0.  The  other  beam,  say  B*.  whose 
power  spectrum  is  obtained  by  eq.  (7 )  letting  Ao=0, 
gives  a  spatially  coherent  but  anticorrelated  (i.e.  with 
a  n  dephasing)  illumina.  on  at  P,  and  P2.  The  an¬ 
ticorrelation  condition  can  be  easily  satisfied  acting 
on  the  angular  position  of  BS  in  order  to  produce  a 
suitable  tilling  of  the  wavefronts  of  B^  with  respect 
to  M.  The  spectrum  of  the  field  produced  by  the  two 
sources  at  a  point  P  located  at  the  same  distance  2? 
from  both  sources  is  then  examined  by  a  spectrum 
analyzer  SA. 

Of  course,  the  two  sources  synthesized  in  this  way 
are  secondary  instead  of  primary  sources,  but  this 
does  not  change  the  basic  mechanism  leading  to  the 
redshift  or  blueshift  phenomenon. 

We  performed  an  experimental  test  of  the  synthe¬ 
sis  procedure  using  a  scheme  derived  from  that  of 
fig.  2,  with  some  unessential  changes  (such  as  re¬ 
placing  holes  by  rectilinear  slits  to  improve  the  lu- 


Fig  2  Experimental  letup 


minosity  of  the  system).  Each  of  the  two  beams  Bn 
and  B;,  spatially  coherent  across  the  distance  from 
P,  to  P2>  was  prepared  in  a  standard  way  by  colli¬ 
mating  the  light  emerging  from  a  sufficiently  narrow 
slit  illuminated  by  a  primary  halogen  lamp.  Filters 
were  inserted  across  each  beam  in  order  to  produce 
the  required  spectral  distribution.  In  a  first  experi¬ 
ment,  aimed  at  evidentiating  the  redshift  phenom¬ 
enon,  the  spectrum  of  Ba  was  centered  at  a 
wavelength  of  526  nm  (giving  tu„= 3.01  x  1  O'5  rad/s) 
with  a  linewidth  of  8  nm,  whereas  the  spectrum  of 
B/  was  centered  at  520  nm  (u^=3.04x  1 0’ 5  rad/s) 
with  a  linewidth  of  10  nm.  In  fig.  3,  curve  (c)  gives 
the  spectrum  produced  at  P  by  a  single  source  (at 
either  P,  or  P2).  Curve  (a)  gives  the  spectrum  ob¬ 
tained  at  P  when  both  sources  are  turned  on.  A  shift 
of  the  spectral  maximum  tow  ards  lower  frequencies 
of  about  I  x  I013  rad/s  is  revealed.  In  order  to  pro¬ 
duce  a  blueshift,  the  spectra  of  B0  and  Be  were  in¬ 
terchanged  {cja- 3  04 x  10'5  rad/s;  cu,=  3.01  X  I0,J 
rad/s).  The  system  was  adjusted  in  such  a  way  that 
the  spectrum  produced  at  P  by  a  single  source  is  still 
represented  by  the  curve  (c).  The  spectrum  ob¬ 
tained  in  this  new  situation  when  both  sources  are 
radiating  is  given  by  curve  (b).  A  blueshift  of  about 
2x  10' 3  rad/s  is  apparent.  The  asymmetry  between 
redshift  and  blueshift  in  the  present  experiment  is  to 


Fig  3  Experiment*!  multi  (heavy  dots )  of  the  measurement  of 
field  spectrum  for  the  redshiftmg  (curve  (a ) )  and  the  blueshift- 
ing  (curve  (b) ).  Curve  (c)  represents  the  spectrum  (multiplied 
by  a  factor  2 )  produced  by  a  single  source.  For  the  sake  of  clarity , 
the  experimental  points  have  been  connected  by  lines. 
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be  escribed  to  differences  of  the  two  spectre!  com¬ 
ponents  (both  in  linewidth  end  in  weight). 

While  this  paper  wes  managed  for  publication,  the 
results  of  a  similar  experiment  were  presented  by 
W.H.  Knox  and  R.S.  Knox  at  the  1 987  Annual  Meet¬ 
ing  of  the  Optical  Society  of  America. 

This  work  was  supported  by  Minister©  Pubblica 
Istruzione  (Italy). 
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Shifts  in  spectral  lines,  due  to  source  correlation  and  violation  of  scaling  law.  have  been  discussed  with  reference  to  spectrora- 
diometnc  measurements.  The  large  scatter  found  in  the  tmercompanson  of  spectroradiometrtc  scales  maintained  b>  standards 
laboratories  is  attributed  to  the  spectral  shifts. 


1.  Introduction 

Recently  Wolf  (1-3]  has  shown  that  statistical 
correlation  of  light  emitted  by  a  partially  coherent 
source  can  produce  a  frequency  shift  in  the  spectrum 
observed  in  the  far  field.  This  shift  is  produced  only 
when  the  correlation  function  of  the  emitted  radia¬ 
tion  does  not  satisfy  a  certain  scaling  law.  A  few  ex¬ 
periments  [4-8]  have  been  performed  to  prove  this 
theoret  cal  prediction  which  has  far  reaching  con¬ 
sequences  in  the  field  of  optical  measurements.  The 
effect  can  be  directly  observed  in  spectroradiometric 
measurements  where  the  spectral  power  of  a  source 
is  measured  as  power  emanating/incident  per  unit 
area  per  unit  wavelength.  To  measure  the  spectral 
power  distribution,  different  apertures  are  placed  in 
the  optical  path  between  the  source  and  the  mono- 
chrorr  tor.  These  apertures  make  the  incoherent  light 
partially  coherent  and  a  source  correlation  is  intro¬ 
duced  which  along  with  the  optics  involved  violates 
certain  scaling  law  thereby  modifying  the  spectral 
distribution  of  the  source  in  the  far  zone.  While  dis¬ 
cussing  the  sources  of  errors  in  speciroradiometry, 
Moore  [9]  has  stated  emphatically  that  the  spectral 
transmission  characteristics  of  monochromators  are 
notoriously  non-uniform.  Nc  suitable  explanation  of 
this  non-uniform  behaviour  of  monochromator 
transmission  was  given.  On  the  basis  of  our  experi¬ 
mental  observations,  it  is  felt  that  shifts  in  the  spec¬ 
tral  lines  due  to  source  correlation  (2,3]  may  have 
some  bearing  with  the  unusual  behaviour  of  the 
monochromator  transmission. 


In  spectroradiometric  measurements  it  is  usual  to 
place  an  integrating  sphere  before  the  entrance  slit  of 
the  monochromator  so  that  radiation  enters  the 
monochromator  only  after  undergoing  multiple  re¬ 
flections  within  the  sphere.  Accurately  measured 
small  apertures  of  mm  order  dimensions  are  put  on 
the  surface  of  the  integrating  sphere.  These  apertures 
and  the  optics  involved  modify  the  degree  of  spectral 
coherence  of  the  light  incident  on  the  monochrom¬ 
ator  and  the  monochromator  behaves  in  an  erratic 
manner. 

This  communication  describes  the  effects  of  ap¬ 
ertures  placed  in  the  optical  path  and  the  observed 
shifts  in  the  source  spectrum.  We  have  also  endea¬ 
voured  to  provide  an  explanation  for  the  unexpected 
nonuniform  behaviour  of  the  monochromator. 


2.  Theory 

Consider  a  quasi-homogencous  source  whose 
spectrum  is  a  line  having  gaussian  profile 

SQ(tu)  =  (l/iyjir) 

X  exp[  —  (a>— o>o):/2<5:)  ,  <5«*ca>0.  (1) 

and  whose  degree  of  spectral  coherence  is  also  gaus¬ 
sian  and  independent  of  ui  viz. 

ik>(r)s=exp(-r:/2<J:)  ,  (2a) 

and 

AjfAJMSv/^O’expt-*1^).  (2b) 
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where  £  is  called  effective  coherence  length  and 
/2q(A)  is  a  three-dimensional  Fourier  transform  of 
fiQ(r).  Wolf  has  shown  that  such  a  source  will  emit 
light  whose  spectrum  in  the  far  zone  will  be  red 
shifted  if  the  degree  of  spectral  coherence  and  co¬ 
herence  length  is  independent  of  frequency  io.  In 
general,  the  spectral  profile  and  degree  of  coherence 
may  not  be  true  gaussian.  In  such  cases  a  shift  may 
occur  on  either  side  of  the  source  spectrum. 

The  spectrum  of  light  in  the  far  zone  generated  by 
a  source  represented  by  eqs.  ( 1 )  and  (2 )  is  given  by 

5T(w)=(aAV'^r) 

X  exp(  -  (cu-tua/cr  )/2(<5/a):]  ,  (3) 

where  ar=  1  +  (<5£/c):,  c  being  the  velocity  of  light 

Eq.  ( 3 )  shows  that  the  emitted  light  in  the  far  zone 
has  a  gaussian  profile.  The  frequency  is  not  centred 
at  u0  but  is  shifted. 

It  is  thus  evident  that  if  some  correlation  exists  in 
the  light  emitted  by  a  source,  the  spectrum  in  the  far 
zone  is  shifted  with  respect  to  the  source  spectrum 
and  the  amount  of  shift  depends  on  the  effective  cor¬ 
relation  length. 


3.  Experimental  results 

To  study  the  effect  of  apertures  placed  in  between 
the  source  and  monochromator  vis-a-vis  the  corre¬ 
lation  properties  of  the  source,  the  experimental  set 
up  (fig  1 )  used  was  an  integrating  sphere  of  0.5  m 
diameter  coated  with  BaS04  paint.  A  450  W  tung¬ 


RLANE-t  PLANE -II  PLANE  Ht 

INPUT  SECONDARY  EAR  FIELD 

Fi|  I .  Experimental  arrangement  for  realizing  secondary  source 
with  varying  degree  of  spectral  coherence 


sten  halogen  lamp  was  placed  at  the  centre  of  the 
sphere.  An  aperture  of  0.24  cm  diameter  was  put  on 
the  surface  of  the  sphere.  After  multiple  reflections, 
the  light  falls  on  this  aperture.  (The  aperture  plane 
is  assigned  as  source  plane  1. )  Different  filter-lens 
and  filter-mirror  combinations  were  used  to  pro¬ 
duce  a  secondary  source  at  the  focus  of  the  combi¬ 
nation  (considered  as  secondary  plane  II).  Aper¬ 
tures  of  varying  sizes  from  0.24  cm  to  1 .0  cm  were 
put  at  the  secondary  plane  II.  The  spectrum  of  light 
in  far  zone  i.e.  in  plane  III,  was  recorded  with  the 
help  of  a  0.5  m  Jarrell- Ash  grating  monochromator 
without  and  with  the  apertures.  It  was  observed  that 
without  the  apertures,  the  spectrum  of  light  in  plane 
III  was  independent  of  the  location  of  the  mono¬ 
chromator.  Insertion  of  an  aperture  in  plane  11  made 
the  spectrum  shift.  Observations  with  varying  ap¬ 
erture  sizes  were  made  and  it  was  observed  that  after 
a  critical  size  of  the  aperture  the  shift  disappeared. 
The  observations  were  repetitive  and  reproducible 
and  are  presented  in  table  1 . 

The  circular  aperture  placed  at  plane  I  produced 
a  secondary  source  at  plane  II.  Various  combinations 
of  filter-lcns/mirror  were  used  to  eliminate  the 
wavelength  dependence  of  the  complex  degree  of 
spectral  coherence  which  is  produced  through  ap¬ 
plication  of  Van  Cittert-Zemike  theorem.  The  illu¬ 
mination  thus  produced  in  plane  II  violated  the  scal¬ 
ing  law.  Introduction  of  a  circular  aperture  in  plane 
II  helped  in  modifying  the  spectrum  in  the  far  zone. 


T»ble  I 

Spectral  shift  variation  with  different  aperture 


Peak 

transmission 

wavelength 

Half 

band¬ 

width 

(nm) 

Shifted  peak  transmission 
wavelength  ( nm )  with 
secondary  aperture 

dia  2  4  mm  dial  0  mm 

422  0 

9 

421.0 

422.0 

484.1 

9 

483.6 

484  1 

512.4 

5 

314.1 

5124 

566.1 

13 

564  1 

566  1 

609.1 

8 

610.3 

609  1 

652.0 

8 

653.2 

652  0 
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4.  Discussion 

The  above  observations  reveal  that  if  the  size  of 
the  aperture  is  very  small  the  shift  in  the  spectral  line 
is  large  and  when  the  aperture  size  is  increased  the 
shift  reduces.  In  all  spectroradiometric  measure¬ 
ments  it  is  a  common  practice  that  small  apertures 
of  the  size  of  mm  are  used  on  the  integrating  sphere. 
To  the  best  of  our  knowledge  no  work  has  been  re¬ 
ported  where  the  partial  coherent  nature  introduced 
by  propagation  through  the  optical  system  and  ren¬ 
dered  to  the  light  emitted  from  the  aperture  of  an 
integrating  sphere  is  taken  into  consideration.  The 
experimental  evidence  [8]  as  well  as  the  theoretical 
stimulations  have,  however,  indicated  that  due  weight 
has  to  bt  given  to  the  correlation  properties  of  light 
emitted  by  the  apertures  and  the  effect  of  correlation 
property  on  the  spectrum  of  the  emitted  light.  If  these 
above  factors  are  considered,  it  is  expected  that  par¬ 
tially  coherent  light  when  viewed  through  interposed 
optics  and  monochromator  will  show  a  spectral  shift. 
The  nature  and  magnitude  of  the  spectral  shift  will 
depend  on  the  optical  path  and  the  violation  of  scal¬ 
ing  law  for  the  spectral  degree  of  coherence.  This 
qualitatively  explains  the  unexpected  behaviour  in 
the  measurement  through  a  monochromator  which 
w  as  earlier  ignored  or  was  taken  as  measurement  un¬ 
certainly.  Far  reaching  consequences  of  this  effect 
can  be  expected  throughout  the  entire  electromag¬ 
netic  spectrum.  For  example  the  large  scatter  in  the 
1974  intercomparison  [9]  of  the  spectral  irradiance 
scale  of  eight  standard  laboratories  could,  to  some 
extent,  be  attributed  to  the  correlation  property  of 
the  light  emitted  from  the  integrating  sphere  and  also 
due  to  different  optical  set  ups  used  by  different  na¬ 
tional  laboratories. 

In  a  senes  of  papers.  Wolf  has  shown  that  if  the 
spectral  degree  of  coherence  follows  scaling  law,  the 
far  field  has  the  same  spectrum  as  that  of  the  non 
coherent  light  It  means  that  for  a  source  which  is 
effectively  spatially  incoherent  i.e.  <*—  0  and 

( * )  —  0.  no  shift  will  occur.  The  above  prediction 
has  been  verified  by  us.  It  is  also  evident  from  the 
theory  that  for  a  source  with  |/r0(r)  |  =  1,  this  shift 
will  be  very  small  but  for  a  partially  coherent  light 
which  violates  scaling  law,  the  shift  is  appreciable.  In 
our  observation  we  have  seen  that  if  the  size  of  the 
secondary  apertures  is  less  than  the  size  of  the  co¬ 


herence  length  produced  by  the  primary  aperture,  the 
shift  is  less  and  if  a  third  aperture  with  a  size  less  than 
the  coherence  length  produced  by  the  secondary  ap¬ 
erture  is  placed  at  plane  III,  no  shift  is  observed.  With 
increase  in  the  size  of  the  aperture  at  the  secondary 
plane  we  have  observed  that  the  shift  increases  and 
that  after  a  certain  size  of  the  aperture  the  shift  dis¬ 
appears.  The  observations  are  consistent  with  the 
theoretical  findings  of  Wolf. 

In  order  to  establish  that  the  size  of  an  aperture 
placed  before  a  radiometer  plays  a  significant  role  in 
spectral  power  measurements,  we  have  taken  obser¬ 
vations  by  placing  a  filter  having  a  !'(/.)  response 
before  the  aperture.  Fig.  2  shows  the  shift  in  the 
spectral  distribution  of  a  lamp- !'(/.)  filter  combi¬ 
nation  in  presence  and  in  absence  of  a  secondary  ap¬ 
erture.  The  observed  spectral  shift  was  10  nm  In  the 
case  of  interference  filters,  the  shifts  observed  were 
of  the  order  of  their  half  bandw  idth.  For  1(A)  filter, 
however,  the  shift  did  not  correspond  to  its  large  half 
bandwidth  but  was  relatively  very  small.  We  atirib- 
ute  the  smaller  shift  with  the  !'(/.)  filter  to  a  partial 
violation  of  the  scaling  law.  It  was  also  found  that, 
if  the  size  of  secondary  aperture  is  more  than  some 
critical  size,  1(A)  spectral  distribution  does  not 
change  during  propagation. 


WAVELENGTH  (nm) 

Fig  2.  Spectral  distribution  of  !'(*)  filler  using  lens  optics  t  X  i 
without  secondary  aperture  and  (B I  with  secondary  aperture 
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It  is.  therefore,  concluded  that  mere  placing  an  ap¬ 
erture  of  any  arbitrary  value  makes  the  spectrora- 
diometric  and  radiometric  measurements  erroneous. 
It  should  be  ensured  before  taking  any  measurement 
that  the  spectral  distribution  of  the  source  and  the 
filters  or  monochromator  combination  should  not 
change  by  the  optics  involved  in  the  path  of  the  light 
beam.  Otherwise  the  measurements  will  not  be  true 
measurements  and  may  result  in  large  deviations 
from  the  true  values. 
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